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ABSTRACT 

The  operational  use  of  military  aircraft  in  forward  area  situations 
has  necessitated  a  comprehensive  look  at  the  critical  factors  which  define 
the  aircraft  flotation  performance  and  operations  capability  on  semi-  and 
unprepared  soil  runways.  This  investigation,  which  is  a  part  of  a  con¬ 
tinuing  program,  wa<*  directed  primarily  at  defining  the  drag-sinkage 
response  of  multiple  wheel  landing  gear  on  soil  and  the  development  of 
multiwheel  flotation  criteria  to  permit  comparative  evaluation  of  the 
relative  merits  of  various  landing  gear  configurations. 

The  total  multiwheel/soil  interaction  study  consisted  of  four  parts. 
Part  1  was  an  evaluation  of  existing  full  scale  field  test  data  from  air¬ 
craft  or  test  carts  operating  with  multiple  tire  configurations  (twin  and 
tandem).  The  second  part  was  a  Twin  Plate  Vertical  Load  (dynamic) 

Test  Program  which  yielded  sinkage  interaction  effects  from  single  and 
twin  plate  load  tests  in  sand  and  clay.  An  analytical  study  was  made  in 
Part  3  using  a  lumped  parameter  iteration  technique  together  with  an 
elastic-plastic  soil  model  for  investigating  adjacent  load  (twin)  interaction 
effects.  Part  4  was  the  Rolling  Multiple  Wheel  Verification  Test  Program. 
The  results  of  the  study  show  that  certain  multiwheel  configurations  and 
spacing  are  beneficial  in  terms  of  sinkage-drag  performance.  The  results 
of  the  multiwheel/ soil  interaction  study  were  used  to  develop  multiwheel 
flotation  criteria  (guidelines  for  evaluating  performance). 

The  flotation  variable  of  braking  was  also  studied  on  a  preliminary 
basis.  Braked  tire/soil  interaction  equations  suitable  for  defining  the 
braking  coefficient  on  soil  were  developed,  and  the  results  of  a  compara¬ 
tive  study  using  these  braked  tire  equations  were  favorable. 
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SECTION  I 

INTRODUCTION  AND  SUMMARY 


Comprehensive  efforts  at  studying  the  problems  associated  with  the 

operation  of  aircraft  on  forward  area  soil  runways  have  been  conducted 

/ «  2  3  4  5) 

over  the  past  several  years'  ’  *  ’  *  .  Those  factors  which  influence  the 

flotation/operation  performance  of  aircraft  on  soil  runways  are  aircraft 
landing  gear  configuration,  tire  inflation  pressure,  aircraft  surface  speed, 
soil  strength  and  texture,  surface  roughness,  turning  effects,  landing 
impact,  braking  action,  multipass  operations  of  the  aircraft.  The  results 
of  past  research^  have  identified  what  have  been  termed  the  primary  and 
secondary  variables  as  related  to  aircraft  flotation/operation  performance. 
The  primary  variables  are  drag,  sinkage,  multiple  wheel  effects,  and 
braking  effects.  Secondary  variables  include  roughness,  speed,  multipass, 
landing  impact,  and  turning. 

In  order  to  develop  the  necessary  criteria  to  permit  the  efficient 
and  effective  utilization  of  soil  runways  while  including  adequate  safety 
provisions,  it  is  necessary  to  define  the  interaction  effects  between  the  tire 
and  soil  during  rolling,  turning,  and  braking  modes.  The  current  research 
effort  described  in  this  report  is  a  part  of  a  comprehensive  aircraft  flotation 
research  program.  The  objective  of  the  program  is  to  analytically  define 
landing  gear-soil  interaction  and  to  develop  a  system  for  rating  the  relative 
flotation  capacity  of  landing  gear  contact  elements  and  landing  gear  systems 
during  aircraft  operation  on  semi-  and  unprepared  soil  runways. 

Flotation  in  this  context  is  defined  as  the  tire -soil  interaction  and  the 
ability  of  the  soil  to  support  the  tire  load  system. 

Phase  1^  of  this  program  included  a  survey  of  the  flotation  problem, 
establishment  of  the  critical  flotation  parameters,  and  an  investigation  of 
available  flotation  data  leading  to  the  development  of  a  flotation  analysis 
equation.  Phase  II  ^  included  the  development  of  an  empirical  sinkage 
prediction  equation,  development  of  a  lumped  parameter  simulation  sinkage 
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prediction  technique  for  moving  tires,  conducting  the  Rolling  Single  Wheel 
Verification  Tests,  and  the  development  of  the  single  wheel  Relative  Merit 
Index  (RMI)  system  for  defining  comparative  flotation  capacity.  Phase  III 
of  the  research  program  consists  of  two  segments.  The  first  segment, 
which  is  detailed  in  this  report  includes: 

-  Multiwheel  sinkage-drag  analysis  equations 

-  Twin  Plate  Vertical  Load  Test  Program 

-  Lumped  parameter  rnultiwheel  simulation  computer 
program  for  twin  and  tandem  configurations 

-  Rolling  Multiwheel  Verification  Tests 

-  Preliminary  studies  of  braking  effects. 

The  second  segment  of  Phase  III  will  concentrate  on  more  detailed  studies 
of  braking  as  well  as  preliminary  evaluations  of  the  effects  of  speed, 
turning,  multipass,  and  roughness  on  flotation  performance. 

The  results  of  the  current  research  program,  as  well  as  the  results 
of  mobility  and  flotation  studies  conducted  by  others,  were  used  in  developing 
the  Quick  Reference  -  Tires  On  Soil  Flotation  Guide  shown  in  Table  1. 
Reference  to  Figure  1  indicates  considerable  progress  has  been  made  in 
recent  years  in  defining  tire/soil  flotation  performance. 
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Quick  Reference-Tires  on  Soil  Flotation  Guide 
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Wolfing  Drag  Z  =  Instantaneous  Soil  Sinkage  R  =  Braked  Drag 

B 

Vertical  u<»ad  D  =■  Tire  Outside  Diameter  b  -  Tire  Se  ction  Width 


SECTION  II 

MULTIWHEEL  SINKAGE  AND  DRAG  ANALYSIS 


1.  General  Considerations  in  Maitiwheel  Performance 

The  growth  in  size  and  weight  of  military  aircraft  with  designed 
operational  capability  on  soil  runways  has  led  to  an  increasing  use  of 
multiple  wheel  landing  gear  configurations.  Possible  multiple  tire  con¬ 
figurations  include  twin  (dual),  tandem-tracking  and  tandem-nontracking 
as  shown  in  Figure  2.  Previous  studies^’  ^  had  shown  that  adjacent  loaded 
areas  influence  the  performance  of  aircraft  tires  operating  on  soil  run¬ 
ways.  The  nature  and  magnitude  of  this  influence,  however,  had  not  been 
established  in  these  studies. 

A  comprehensive  four-part  program  was  undertaken  to  study  the 
tire-soil  sinkage  and  drag  interaction  effects  for  multiple  tire  landing 
gears.  Part  1  consisted  of  an  evaluation  of  existing  full  scale  field  test 
data  from  aircraft,  or  test  carts  operating  with  multiple  tire  configurations. 
The  second  part  was  a  Twin  Plate  Vertical  Load  Test  Program  in  which 
both  single  and  twin  plates  were  dynamically  penetrated  into  cohesive 
(clay)  and  cohesionless  (sand)  type  soils.  An  analytical  study  was  made 
in  Part  3  using  a  lumped  parameter  iteration  technique  together  with  an 
elastic-plastic  soil  model  (soil  simulation)  for  studying  adjacent  load 
interaction  effects.  Part  4  was  the  full  scale  P„olling  Multiwheel  Veri¬ 
fication  Tests. 

2.  PART  1  -  Review  of  Existing  Multiwheel  Test  Data 

(4  6  7  8  9) 

Availability  of  existing  multiple  tire  test  data'  ’  ’  ’  ’  1  was 
restricted  to  cohesive  soils.  A  summary  of  available  data  for  the  drag 
performance  of  twin  tires  on  clay  at  various  tire  spacings  is  shown  in 
Figure  3.  The  comparison  is  made  between  the  drag  ratio,  R/P  (R  =  total 
rolling  drag  for  both  tires,  P  =  total  vertical  load  on  both  tires),  for  the 
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Twin  Wheel  Configuration 
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Tandem  Wheel  Configuration 


b  =  Tire  Section  Width 

D  =  Tire  Outside  Diameter 

m  =  Fraction  of  Tire  Outside  Diameter 

n  =  Fraction  of  Tire  Section  Width 

r'  =  Fraction  of  Tandem-Nontracking 
Tire  Spacing 

r'  =  ■  —  1  where  L  =  "^(nb)2+(mD)2 

V2  Z 
vb  +D 


Figure  Z.  Multiple  Tire  Configurations 
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Figure  3.  Twin  Tire  Performance  on  Clay 


1.0 
.9 
.8 
.7 
6 

per  Tire  Tandem 

7 -  .5 

~  Single  Tire  ^ 

.3 
.2 
.1 
0 

Ratio  of  Tandem-Tracking  Spacing  to  Tire  Diameters, 


i S  56X16  TIRES,  BUCKSHOT  CLAY 

a  20.00-20  TIRES,  BUCKSHOT  CLAY 

0  10.00-20  TIRES,  CLAY 

V  33.75X9.75  TIRES,  LONDON  CLAY 

•  40X14  TIRES,  HARPERS  LAKE  CLAY 

R  =  ROLLING  DRAG 
P  »  VERTICAL  LOAD 

J _ 1  ‘  1 _ I _ I _ I _ L_ 

10  2.0  30  4.0 


Figure  4.  Tandem-Tracking  Performance  on  Clay 


m 


7 


twin  tire  case,  and  the  drag  ratio  for  the  single  (isolated)  tire  case.  The 
test  data  shown  in  Figure  3  was  primarily  for  low  sinkages  (<3/4").  A 
similar  drag  ratio  comparison  was  made  using  available  tandem-tracking 
test  data  as  shown  in  Figure  4.  Again  the  test  data  was  primarily  for  low 
sinkages. 

Reference  to  Figures  3  and  4  shows  considerable  scatter  in  the  test 
data  but  generally  indicates  that  within  certain  spacing  limitations,  some 
reduction  in  rolling  drag  can  be  expected  for  twin  and/or  tandem-tracking 
tire  configurations  over  that  for  single  tires  when  operated  on  cohesive 
soils  at  low  sinkages.  Since  drag  has  been  shown  to  be  directly  related 
to  sinkage  through  the  relationship 

|  =  0.018  +  3.23  (|)  (1) 


where 

Z  =  instantaneous  sinkage 
D  =  tire  diameter, 

corresponding  reductions  in  sinkage  can  be  inferred  from  the  use  of 
Equation  (1).  No  usable  test  data  was  found  for  multiple  tires  operated 
on  cohesionless  soils. 

3.  PART  2  -  Twin  Plate-Vertical  Load  Test  Program 
Objective 

The  purpose  of  the  dual  vertical  plate  tests  was  to  obtain  initial 
information  and  data  concerning  the  sinkage  interaction  effects  produced 
by  adjacent  dynamic  plate  loads  as  compared  to  a  single  isolated  plate 
load  under  similar  test  conditions.  Plate  spacing  was  also  varied  to 
permit  a  detailed  look  at  the  influence  of  this  variable.  While  dynamic 
vertical  plate  load  tests  do  not  simulate  the  action  of  a  rolling  tire,  such 
tests  can  be  used  to  study  the  influence  of  sinkage  interaction.  As  indi¬ 
cated  by  Equation  (1),  since  sinkage  is  directly  related  to  drag  for  rolling 
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tires,  the  sinkage  results  of  the  plate  tests  can  be  used  to  interpret  drag 
effects  by  inference. 

Test  Program 

To  accomplish  the  objectives,  a  series  of  single  and  twin  vertical 
load  plate  tests  were  conducted  in  two  different  soils:  buckshot  clay  and  a  mix¬ 
ture  of  mortar  sand  and  Yuma  sand.  Dynamic  loads  were  applied  through  the 
plates  at  a  constant  rate  of  penetration.  Sufficient  cone  penetrometer  tests 
were  made  to  fully  define  the  soil  properties  of  each  test  specimen.  A 
summary  of  the  test  program  is  given  in  Table  2.  The  full  details  of  the 
soil  properties,  soil  preparation,  test  procedure,  and  uniformity  of 
placement  are  given  in  Appendix  I. 

Test  Setup 

A.  Instrumentation 

Dynamic  vertical  loads  were  applied  to  the  test  specimens  by 
an  MTS  5000  pound  dynamic  load  system  containing  two  service  modules 
and  a  load  frame.  The  soil  was  contained  in  a  test  box  of  the  dimensions 
16  inches  by  13  inches  by  9  inches  deep.  The  complete  test  system, 
including  the  soil  test  box,  is  shown  in  Figure  5.  Loads  were  monitored 
by  a  SR-4  Baldwin-Lima-Hamilton  0.  5KIP  load  cell  attached  to  the  loading 
head.  Both  loads  and  sinkages  were  recorded  on  a  Honeywell  Visicorder. 

All  plate  tests  in  sand  were  conducted  at  a  rate  of  12.  5  inches 
per  second,  all  plate  tests  in  clay  at  a  4.  2  inches  per  second  velocity, 
and  cone  penetrometer  tests  were  run  at  a  rate  of  1.  25  inches  per  second. 

B.  Test  Plates 

The  application  of  the  load  to  the  soil  was  accomplished  through 
the  use  of  circular  plates,  1-1/2  inches  in  diameter  and  5/16  inch  thick. 

The  twin  plate  tests  used  two  such  plates  mounted  in  an  aluminum  block 
with  a  shaft  which  permitted  varying  the  plate  spacing.  Figure  6  shows 
the  single  and  dual  plate  test  equipment. 
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TABLE  2 


DUAL  VERTICAL  PLATE  TEST  PROGRAM  (DESIGN) 


Test 

Soil 

Soil  Stiffness 

Plate  Spacing,  nD 

Plate  Re; 

No. 

Type 

Cl  (lb) 

avg 

p 

Tests 

IS 

Sand* 

9.0 

Cone  Only 

2S 

II 

9.0 

Single 

3S 

II 

9.0  . 

Single 

4S 

tl 

9.0 

Cone  Only 

5S 

II 

9.0 

Single 

6S 

II 

9.0 

Cone  Only 

7S 

II 

9.0 

1-1/2  D 

8S 

II 

9.0 

1-1/2  D 

Repeat 

9S 

It 

9.0 

Cone  Only 

10S 

II 

9.0 

2-1/2  D 

Repeat 

1  IS 

II 

9.0 

2-1/2  D 

Repeat 

12S 

II 

9.0 

4  D 

1C 

Buckshot  Clay  30-35 

Single 

2C 

II 

20-25 

1-1/2  D 

Repeat 

3C 

It 

30-35 

2-1/2D 

Repeat 

4C 

II 

20-25 

Single 

5C 

II 

30-35 

4  D 

Repeat 

6C 

II 

30-35 

2-1/2  D 

Repeat 

7  C 

?l 

30-35 

Single 

8C 

1  1 

20-25 

1-1/2  D 

Repeat 

9C 

II 

20-25 

4  D 

Repeat 

*  Mixture  of  Mortar  Sand  and  Yuma  Sand 

n  =  Plate  Spacing/ Plate  Diameter 

Dp=  Plate  Diameter  =  1.50" 

Cl  =  Average  Cone  Index 
avg 
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Figure  5.  Plate  Test  System 
not  reproducible 


C.  Test  Penetrometer 


A  cone  penetrometer  was  used  to  measure  the  soil  Jtrength 
and  uniformity  after  placement.  The  cone  was  made  of  stainless  steel, 
ground  to  a  smooth  finish,  and  was  of  the  same  dimensions  as  the  Army 
cone  penetrometer  used  in  mobility  studies.  The  cone  penetrometer  is 
shown  in  Figure  6. 

Sand  T ests 


A  typical  cone  and  single  plate  penetration  resistance  recording 

is  shown  in  Figures  7  and  8.  The  results  of  the  cone  penetration 

tests  which  reflect  the  uniformity  of  the  sand  test  sections  are  given  in 

Appendix  1  along  with  the  formula  used  in  calculating  the  average  soil 

cone  penetration  resistance  (CPR  )  for  each  test  section.  The  results 

avg 

of  the  single  and  twin  plate  tests  conducted  in  the  sand  sections  ai  e  given 

TABLE  3 

PLATE  FORCE  RESISTANCE  DATA,  SAND 
(Initial  Load) 


Test  Force  Resistance  per  Plate 

Sequence  Comments  in  lb  at  Penetration  of 

No.  1/4’  1/2"  3/4"  1" 


2S 

1st  Single  Plate  Test 

14.  3 

13.  5 

17.  6 

20.7 

3S 

2nd  Single  Plate  Test 

16. 

18.  1 

21.  1 

25.  0 

SS 

3rd  Single  Plate  Test 

6.  o 

9.  2 

12.  5 

16.0 

7S 

1st  Twin  at  1-1/2  Dp 

8.  5 

10.  2 

11.  8 

14.  5 

8S 

2nd  Tv'in  at  1-1/2  Dp 

12.6 

15.  6 

17.  6 

ro 

o 

rv 

10S 

1st  Twin  at  2-1/2  Dp 

6.  9 

9.  2 

11.  9 

15.  1 

1 1  s 

2nd  Twin  at  2-1/2  Dp 

10.0 

12.  7 

14.  8 

17.0 

12S 

1st  Twin  at  4  Dp 

8.  9 

10.  9 

13.  1 

15.  8 

in  Table  3.  Figure  9  shows  this  same  data  presented  as  the  average 

penetration  resistance  (CPR  )  versus  sinkagc  for  each  plate  spacing. 

a  v  c: 

F  was  determined  by  the  results  of  three  tests  (2S,  3S,  and  5S)  for 
avg 

the  single  plate,  and  two  tests  for  the  twin  plate  tests,  at  spacings  of 


Sinkage 


01  SEC 


Figure  7.  Typical  Cone  Penetration  Test  Result,  Sand 


Figure  8.  Typical  Single  Plate  Test  Result,  Sand 


m  Resistence  Per  Plate 


nd  Tests 


1-1/2  D  and  2-1/2  D  .  Only  one  twin  plate  test  was  conducted  at  a 
P  P 

spacing  of  4  D  .  Analysis  of  Figure  9  shows  that  more  force  (per  plate) 

P 

is  required  in  a  twin  plate  situation  for  a  given  sinkage  than  for  a  cor¬ 
responding  single  piate  penetrated  to  the  same  sinkage,  with  the  exception 
of  the  closest  plate  spacing,  1-1/2  Dp. 

Clay  Tests 

Figures  10  and  11  show  the  typical  recordings  for  a  cone  penetration 

test  and  a  twin  plate  test,  respectively,  in  clay.  The  results  of  the  cone 

penetration  tests  which  reflect  the  uniformity  of  the  clay  test  sections  are 

given  in  Appendix  I  along  with  means  of  determining  the  average  cone  penetration 

resistence  (CPR  for  each  test  section.  The  results  of  the  single  and  twin 
umf 

plate  penetration  tests  are  summarized  in  Table  4.  These  results  along  with 


TABLE  4 

PLATE  FORCE  RESISTANCE  DATA,  CLAY 
(Initial  Loading) 


Test  Force  Resistance  per  Plate 

Sequence  Comments  in  lb.  at  Penetration  of 


No. 

1/4" 

1/2" 

3/4" 

1" 

1C 

1st  Single  Plate  Test 

84.0 

95.  2 

97.  4 

101.  5 

2C 

1st  Twin  at  1-1/2  Dp 

55.0 

61.  5 

63.  5 

64.  3 

3C 

1st  Twin  at  2-1/2  Dp 

74.  3 

80.  5 

82.  0 

81.  5 

4C 

2nd  Single  Plate  Test 

51.  9 

57.  8 

61.  5 

63.6 

5C 

1st  Twin  at  4  Dp 

74.  8 

83.  0 

81.6 

80.  6 

6C 

2nd  Twin  at  2-1/2  Dp 

77.  8 

85.  8 

86.  5 

86.  5 

7C 

3rd  Single  Plate  Test 

87.  7 

97.  9 

98.  9 

98.  6 

8C 

2nd  Twin  at  1-1/2  Dp 

56.  5 

59.  5 

64.  0 

67.0 

9C 

2nd  Twin  at  4  Dp 

74.  0 

82.  2 

80.  9 

80.0 

the  soil  strength  data  were  used  to  analyze  the  results  in  a  manner  similar 

to  that  used  for  sand  as  given  in  Figure  12  which  shows  the  variation  of 

/CPR^.j  for  each  plate  spacing  including  the  single  plate  condition. 

Each  Fa  represents  the  average  of  the  results  presented  in  Table  4, 

Avg  r  r 

and  CPR  is  defined  in  Appendix  I.  Reference  tc  Figure  12  for  clay  shows 


15 


6 


NOTE:  Intersecting  lines  were  caused  by  nonuniformities 
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Figure  12.  F^^/CPR^^  vs.  Sinkage,  Clay  Tests 


a  trend  similar  to  that  for  sand,  that  is,  within  certain  spacing  limitations, 
more  force  is  required  per  plate  for  dual  plate  configurations  at  a  particular 
sinkage  than  is  required  for  a  single  isolated  plate  at  the  same  sinkage. 

The  tendency  of  the  plate  spacing  lines  (see  Figure  12)  to  intersect  was 
caused  by  variations  in  strength  with  depth  in  the  buckshot  clay  test  speci¬ 
mens.  Tests  involving  the  single  plate  and  the  plates  at  a  spacing  of  1-1/2  Dp 
had  test  sections  exhibiting  a  slight  increase  in  strength  with  depth,  while 
tests  involving  plate  spacings  of  2-1/2  Dp  and  4  Dp  had  test  sections  showing 
a  slight  decrease  in  strength  with  depth. 

Analysis  and  Interpretation  of  Test  Results 

The  results  of  the  sand  and  clay  tests  were  previously  interpreted 

in  terms  of  plate  resistance  to  penetration  at  different  sinkages.  A 

clearer  picture  of  sinkage  interaction  effects  is  gained  by  plotting  sinkages 

at  a  common  dimensionless  plate  force  resistance  (F  /CPR  .,).  Figure  13, 

r  avg  umf 

which  was  developed  from  Figures  9  and  12  for  F  /CPR  =1.5  for  sand 
r  °  avg  avg 

and  F  /CPR  . ,  =  2.9  for  clay,  shows  this  interaction.  Reference  to 
avg  umf 

Figure  13  and  use  of  Equation  (1)  would  suggest  approximately  a  15%  to 
20%  reduction  in  rolling  drag  for  optimum  spacing  of  twin  tires  on  sand 
while  optimum  spacing  of  twin  tires  in  clay  would  lead  to  an  approximate 
20%  to  25%  drag  reduction  for  relatively  low  sinkages  (<3/4").  It  should 
be  recognized  that  such  conclusions  are  only  approximate  since  plate 
size  (diameter)  and  magnitude  of  sinkage  are  very  likely  influencing  factors. 

The  reduction  in  sinkage  caused  by  the  proximity  of  an  adjacent 
loaded  area  is  caused  by  an  interference  between  the  two  plastic  zones 
associated  with  each  surface  load.  Figure  14a  shows  qualitatively  the 
elements  of  a  bearing  capacity  theory  as  developed  by  Terzaghi* 1 2). 

Reference  to  Figure  14a  indicates  a  triangular  wedge  of  soil  undergoing 
upward  displacement  as  indicated  by  the  dashed  line.  As  the  second 
loaded  area  P^  moves  toward  the  first  loaded  region  Pj,  the  plastic  zones 
begin  to  overlap  creating  interference  in  the  slip  fields.  As  indicated  by 
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Influence  of  Sinkage  Interaction,  Sand  and  Clay 


(a)  Terzaghi  Elements  of  Bearing  Capacity 
Depicting  Development  of  Slip  Fields 


*—  Spacing  <1.5  Dp—* 


(b)  Adjacent  Load  Aieas  Acting  as  a  Single  Load 


Figure  14.  Adjacent  Loading  Plastic  Zone  Interaction 
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the  results  of  the  plate  tests,  within  certain  spacing  limitation  this  inter¬ 
ference  of  slip  fields  is  beneficial  in  reducing  sinkage  (and  drag  reduction 
if  inference  of  rolling  tire  performance  is  valid).  For  spacings  less  than 
approximately  1-1/2  times  the  plate  diameter,  the  adjacent  load  effects 
are  detrimental  since  the  two  loaded  areas  tend  to  act  as  a  single  load  as 
indicated  in  Figure  14b. 

In  general,  then,  the  results  of  the  vertical  twin-plate  tests  con¬ 
ducted  in  sand  and  clay  have  shown: 

1.  Within  certain  spacing  limitations,  the  presence  of  an  adjacent 
loaded  area  will  cause  a  reduction  in  sinkage  (over  that  which 
occurs  for  a  single  plate)  for  plates  under  similar  vertical 
load  conditions. 

2.  Optimum  spacing  very  likely  exists  leading  to  the  maximum 
reduction  in  sinkage.  This  phenomenon  is  suggested  by  both 
the  plate  tests  results  and  an  analysis  of  bearing  capacity 
theory. 

If  the  performance  of  rolling  twin  tires  can  be  inferred  from  the 
results  of  plate  tests  based  on  sinkage  interaction,  then  reductions  in 
rolling  tire  drag  (for  relatively  low  sinkages)  can  be  expected  through  use 
of  twin  type  tire  configurations  with  proper  spacing. 

4.  PART  3  -  Analytical  Approach 

A.  Twin  Wheel 

The  analytical  approach  employed  in  studying  the  multiwheel 
effects  on  sinkage  and  drag  is  similar  to  that  employed  previously  in 
the  single  wheel  analysis^  .  Multiwheel  effects  on  the  sinkage  of  the  tire 
into  the  soil  were  first  studied  and  then  related  back  to  drag  by  use  of 
previously  established  drag-sinkage  relationships. 

The  sinkage  of  a  multiwheel  landing  gear  configuration  is 
different  from  that  of  a  single  wheel  primarily  due  to  the  difference 
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in  stress  distribution  and  yielding  patterns  within  the  soil  medium  resulting 
from  the  interactions  of  the  other  tires.  For  a  twin  wheel  configuration, 
since  both  wheels  are  moving  about  the  same  wheel  axis,  the  interaction 
effect  is  instantaneous.  For  tandem-tracking  wheel  configuration,  since 
one  wheel  is  moving  in  front  of  the  other,  the  interaction  effect  is  over  a 
certain  time  period;  the  stress  and  yielding  pattern  in  a  region  of  soil 
rolled  over  by  the  front  wheel  changes  over  the  time  period  it  takes  for 
the  trailing  wheel  to  arrive  in  the  region  for  interaction. 

As  an  initial  approach  to  studying  the  multiwheel  interaction 
effects,  the  sinkage  prediction  technique  for  the  stationary  vertical  pulse 
load  problem  was  first  developed.  This  technique  for  a  vertical  pulse  load 
problem  is  similar  to  the  one  considered  for  the  single  wheel  analyzed  in 
Phase  II  ^  of  this  research  program,  except  that  two  surface  contact 
pressure  distributions  instead  of  one  are  applied  on  the  surface.  Since 
the  interaction  effect  is  simultaneous,  this  problem  simulates  the  twin  wheel 
configuration.  Tandem  wheel  configuration  will  be  considered  in  a  later 
section  in  which  rolling  multiwheels  are  studied. 

In  the  following  paragraphs,  the  two-dimensional  plane  strain 
approximation,  the  multiwheel  effects  evaluation,  and  the  method  of 
solution  of  the  vertical  pulse  load  problem  are  discussed. 

Two-Dimensional  Plane  Strain  Approximations 

The  vertical  pulse  load  problem  with  two  surface  contact 
pressure  distributions  is  a  three-dimensional  problem  and  its  solution 
by  the  lumped  parameter  iteration  method  or  the  finite  element  method 
would  require  the  use  of  a  three-dimensional  space  mesh.  Comparison 
of  this  problem  with  the  three-dimensional  moving  single  wheel  problem 
considered  in  Phase  indicated  that  the  computer  run  time  required 

for  this  analysis  would  be  quite  excessive  in  relation  to  the  benefit  in 
results  that  might  be  obtained,  and  the  problem  with  computer  core 
capacity  would  be  difficult  to  overcome.  Therefore,  the  two-dimensional 
plane  strain  approximated  problem  was  used  in  the  analytical  development. 
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The  two-dimensional  problem  is  much  less  difficult  in  com¬ 
parison  to  the  three-dimensional  problem,  and  yet,  can  give  a  reliable 
indication  of  the  influence  of  the  stress  distribution  within  the  soil  on 
the  sinkage.  The  reason  is  that  the  stress  distribution  in  the  region 
immediately  between  the  tires  is  the  major  influence  on  the  sinkage  due 
to  the  twin  wheel  configuration,  and  the  stress  distribution  in  this  region 
for  the  three-dimensional  problem  is  quite  similar  to  that  for  the  two- 
dimensional  problem. 

The  Two-Dimensional  Multiwheel  Vertical  Pulse  Load  Problem 

The  two-dimensional  problem  considered  is  a  semi-infinite 
half-space  soil  medium  with  two  infinitely  long  uniformly  distributed 
pressure  strips  applied  on  the  surface.  The  magnitude  of  the  uniform 
pressure  varies  with  time  and  was  taken  to  be  equal  to  the  average  pres¬ 
sure  experienced  by  a  point  near  the  surface  of  the  soil  when  a  rolling 
tire  moves  across  the  point.  The  pressure-time  dependence  has  been 
measured  experimentally  for  land  vehicles^’ ^  ^  and  has  the  form  of 
a  pressure  pulse  as  shown  in  Figure  15.  This  pressure-time  curve  is 
the  same  as  the  one  used  in  the  single  wheel  analysis^.  A  sectional 
view  of  the  soil  medium  with  the  applied  pressure  is  shown  in  Figure  16. 

In  this  figure,  the  width  of  each  pressure  distribution,  b,  is  equal  to  the 
tire  footprint  width  of  the  twin  wheel  configuration.  It  can  be  seen  that 
the  x=0  plane  is  a  plane  of  symmetry;  therefore,  it  is  necessary  only  to 
consider  half  of  the  soil  medium  and  use  symmetry  boundary  conditions 
at  the  x=0  plane. 

The  Elastic-Plastic  Soil  Medium 

A  review  of  the  existing  soil  models  indicated  that  the  most 
suitable  soil  model  for  the  present  analysis  is  an  elastic-perfectly  plastic 
material  with  the  elastic  deformations  governed  by  Hooke's  law,  the  plastic 


Figure  15.  Typical  Vertical  J.m~'tlse  Pressure  Curve 


Soil  Medium  With  Applied  Pressure  Strips  for  Vertical  Load  Problem 


(1  ^formations  governed  by  an  incremental  stress-strain  relation  which 
is  based  on  the  normality  flow  rule,  and  the  plastic  yielding  governed  by 
the  Drucker-Prager  yield  criterion^  The  soil  parameters  of 

this  model  consist  of  the  elastic  Young’s  modulus,  the  cohesion,  and  the 
friction  angle.  These  parameters  can  be  determined  experimentally 
from  soil  samples  so  that  analytical  results  may  be  compared  with  ex¬ 
perimental  results.  Strain-hardening  was  not  incorporated  in  this  model 
because  no  satisfactory  theory  with  experimentally  determinable  parameters 
is  currently  available. 

Method  of  Evaluation  of  the  Multiwheel  Effects 

The  multiwheel  effects  were  evaluated  from  the  results  of 
the  two-dimensional  vertical  pulse  load  problem  in  the  following  manner. 
First,  the  sinkages  due  to  a  specific  vertical  load  for  various  values  of 
tire  spacing,  nb  (see  Figure  16),  were  obtained.  The  spacing  was  ex¬ 
pressed  as  nb,  where  b  is  the  tire  width  and  n  is  a  numerical  factor 
greater  than  or  equal  to  one.  For  example,  spacings  for  twin  wheel  of 
n=  1.6,  2.3,  3.1,  3.9,  4.7,  5.4,  and  infinity  were  analyzed,  where  the 
case  for  n=°°  corresponds  to  the  single  wheel  problem  and  is  the  same 
as  the  case  for  n=0.  The  multiwheel  effects  were  then  evaluated  by 
studying  the  ratio  of  the  sinkages  for  n?!®  (twin  wheel)  and  that  for 
n=®  (single  wheel). 

Method  of  Solution  -  Lumped  Parameter  Iteration  Method 

In  the  lumped  parameter  iteration  method,  a  lumped  parameter 
model  is  used  for  developing  the  governing  equations.  The  model  chosen 
for  the  two-dimensional  plane  strain  problem  is  quite  similar  to  that  for 
the  axisymmetric  problem  used  in  the  single  wheel  analysis.  It  also  is 
composed  of  alternately  interwoven  mass  and  stress  points,  except  that 
the  points  are  oriented  at  a  45  degree  angle  from  the  soil  surface  as  shown 
in  Figure  17.  A  (r;,  ()  coordinate  system  is  used  for  designating  the 
diagonal  directions,  and  the  ^-coordinate  is  used  for  indicating  the 
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direction  normal  to  the  n  and  '  coordinates  (that  is,  normal  to  the  plane 
of  the  paper).  The  displacements  in  the  n  and  Q  directions  are  IT  and  V, 
respectively.  Then,  the  displacements,  u  and  w,  in  the  x  and  z  directions 
are  related  to  U  and  V  as  follows: 


VT  , 

u  =  U  cos45  0  -  V  sin45  *  =  (U-V) 

(2) 

VT 

w  =  U  sin45  °  +  V  cos45  5  =  (U+V) 

(3) 

and  the  stresses  o  ,  a  ,  and  t  are  related  to  o  ,  o„,  and  t  ,  by  the 
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transformation  equations, 
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and 

op=v(oTi  +  oc)  . 


(4) 

(5) 

(6) 


(7) 


The  pressure  boundary  conditions  are  applied  on  the  soil  surface,  z=0, 

through  a  row  of  fictitious  stress  points  (Row  j=l  in  Figure  17).  Outside 

the  loaded  surface,  the  normal  and  shear  stresses  of  the  fictitious  stress 

points  are  zero.  At  the  fictitious  stress  points  which  are  above  the  loaded 

area,  o  =-p,  0  =0,  and  t  =0,  where  the  prescribed  surface  shear  stress 
z  x  xz 

was  assumed  to  be  zero  for  the  present  problem.  With  the  inversion  of 
the  transformation  equations,  Equations  (4),  (5),  and  (6),  the  stresses  at 
the  loaded  fictitious  stress  points  in  the  45  degree  directions  were  obtained 
as  follows: 


o 

h 


a  +  c 

X  z 


2 


+  T 

XZ 


£ 

2 


(8) 


28 


0  +0 

X  z 


-  T 


£ 

2 


(9) 


o 


C 


2 


xz 


T 


•nC 


£ 

2 


(10) 


where  p  is  the  prescribed  surface  pressure  according  to  the  pressure 
pulse  curve. 

An  advantage  of  this  plane  strain  model  over  the  axisymmetric 
model  is  that  there  are  no  stress  points  at  the  surface  plane  whose  vertical 
stresses  must  be  approximated. 

Governing  Equations,  Numerical  Procedures,  and  Computer 

Programs 

The  governing  equations  for  this  problem  are  the  dynamic 
equations  of  motion,  the  quadrature  equations,  the  yield  criterion,  the 
incremental  strain -displacement  relations,  the  incremental  stress-strain 
relations,  and  the  equations  for  stress  correction.  The  derivations  of 
these  equations  and  the  development  of  the  numerical  procedures  for  the 
plane  strain  lumped  parameter  model  follow  the  same  basic  approach  as 
those  for  the  axisymmetric  model  used  in  the  single  wheel  analysis,  so 
they  are  not  repeated  here.  The  reader  is  referred  to  Reference  13  for 
the  approach  used  in  the  derivations  and  the  development  of  the  numerical 
procedures.  For  convenience,  the  governing  equations  and  the  numerical 
procedures  are  presented  without  derivation  in  Appendix  II.  A  computer 
program  was  written  based  on  these  governing  equations  and  numerical 
procedures.  A  general  flow  chart  of  the  computer  program  is  shown  in 
Figure  18.  A  Fortran  IV  source  program  listing  of  the  computer  program, 
a  sample  set  of  input  data,  a  list  of  definition  of  symbols,  and  some 
remarks  about  running  the  computer  program  are  given  in  Appendix  III. 

Test  Cases  and  Results 

Seven  cases  were  run  with  the  computer  program  developed 
for  twin  wheel  simulation.  All  the  cases  have  the  same  soil,  load,  and 
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computational  parameters,  which  are  listed  below.  This  set  of  parameters 
is  that  of  a  typical  multiwheel  aircraft  tire-soil  interaction. 


Soil  Parameters: 

Density 

Poisson's  ratio 
Young's  Modulus 
Cohesion 
Friction  angle 
Yield  stress  in  shear 


p  =  130  lb.  /cu.  ft. 
v  =  0.  45 
E  =  8950  psi 
c  =  2000  psf 
CO  =  15  0 
k  =  2440  psf 


(These  soil  parameters  correspond  approximately  to  a  clay  soil  with 
CBR  =  8  to  10. ) 


Load  Parameters: 

Tire  Footprint  Width 
Peak  Surface  Pressure 
Time  Duration  of  Pulse 


b  =  1 2.  0  inches 

p  =  24600  psf 
max  r 

t  =  0.  05  sec. 
d 


Computational  Parameters: 
Space  Mesh  Size 
Time  Increments 
Finite  Boundary  Size 


h  =  3.  0  inches 
At  =  6.  25  x  10  ^ 

Depth  =  60" 

Width  =  69"  to  111"  from  the  line  of 
symmetry 


The  only  difference  between  each  of  the  cases  is  wheel 
spacing.  One  of  the  cases  corresponds  to  the  single  wheel  case,  ri=°° 

(or  ri=0),  and  the  rest  of  the  cases  correspond  to  twin  wheel  configurations 
with  spacing  ratios  of  rf=1.6,  2,  3,  3.  1,‘  3.9,  4.7,  and  5.4.  In  order  to 
minimize  the  influence  of  the  finite  boundary  on  the  sinkage  when  the 
wheel  spacings  were  changed,  the  distance  between  the  finite  boundary  and 
the  edge  of  the  applied  pressure  strip  was  maintained  constant  by  changing 
the  distance  of  the  finite  boundary  from  the  line  of  symmetry.  The  peak 
instantaneous  sinkages  for  each  case  and  their  comparison  with  that  for 
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the  single  wheel  case  are  shown  in  Table  5.  Sinkage  in  this  analysis 
was  taken  as  the  average  of  the  vertical  displacements  of  the  mass  points 
which  are  immediately  under  the  applied  surface  pressure  distribution. 


TABLE  5 

TWIN  WHEEL  SINKAGES  FOR  DIFFERENT  WHEEL  SPACINGS 


Case 

Spacing 

Ratio 

n' 

Instantaneous 
Sinkage,  Z 
(inches) 

z  /  z 

Twin  Single 

Dimensionless 
Sinkage,  Z/b 

1 

CD 

0.48 

1.000 

0.0403 

2 

1.6 

0.  51 

1.  050 

0.0423 

3 

2.  3 

0.45 

0.942 

0.0379 

4 

3.  1 

0.  43 

0.  886 

0.0357 

5 

3.  9 

0.42 

0.  871 

0.0351 

6 

4.  7 

0.43 

0.  882 

0.0355 

7 

5.  4 

0.43 

0.  900 

0.  0362 

In  Figure  19.  the  sinkage  for  twin  wheel,  normalized  with 
respect  to  the  sinkage  for  single  wheel  (n^oo),  is  plotted  against  the  wheel 
spacing  ratio,  n1.  This  curve  shows  that  as  wheel  spacing  decreases 
from  infinity,  the  sinkage  gradually  decreases  due  to  the  interaction 
between  the  two  tires,  and  reaches  a  minimum  at  a  particular  wheel 
spacing.  The  wheel  spacing  corresponding  to  the  minimum  sinkage 
probably  depends  cn  many  factors,  for  example,  intensity  of  wheel  load, 
size  of  contact  area,  soil  strength,  etc.  For  the  particular  loads  and 
soil  conditions  of  these  test  cases,  the  minimum  falls  between  n'=3. 0 
andn'=4.5.  Further  decrease  of  the  wheel  spacing  from  the  minimum 
sinkage  point  leads  to  higher  sinkage  because  the  soil  medium  between 
and  under  the  twin  contact  area  starts  to  interact  strongly  and  begins  to 
act  as  a  whole  mass  similar  to  that  underneath  a  single  wheel. 

The  trend  of  the  curve  shown  in  Figure  19  was  also  exhibited 
in  the  experimental  results  obtained  from  the  dual  circular  plate  vertical 
load  test  given  in  Section II  -  Part  2  and  shown  in  Figure  13.  The  plate 
tests  show'ed  the  minimum  sinkage  as  occurring  between  n’=2  and  n'  =  3. 
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Consideration  of  the  ideal  slip  line  field  using  Terzahgi  bearing 
capacity  theory^  ^  also  indicated  that  the  minimum  si.nkage  probably  falls 
between  n'=2.  5  and  n'=3.  5  for  clays. 

The  vertical  displacements  of  the  surface  mass  points  were 
plotted  against  the  horizontal  distance  from  the  axis  of  symmetry  for 
several  times  during  the  downward  movement  of  the  soil  surface.  These 
graphs  are  shown  in  Figure  20,  21,  22,  and  23  for  n'a®,  nU1.6,  n'=3.  1, 
and  n'=4.7.  The  comparison  of  these  graphs  shows  the  difference  of  the 
patterns  of  soil  surface  deflection,  ws,  for  the  various  wheel  spacings.  In 
particular,  the  deflection  pattern  for  n'  =  l.6  is  more  similar  to  the  pattern 
for  a  single  large  contact  area  than  that  for  n'=3.  1.  This  shows  the  reason 
why  the  sinkages  for  n  between  1.  0  and  2.  0  tend  „o  be  higher  than  that 
for  n'=3. 1. 

Some  preliminary  test  cases^^  run  with  a  smaller  finite 
boundary  size  and  a  smaller  grid  size  show  the  same  results  as  those 
reported  here,  but  because  of  the  limited  boundary  size,  some  influence 
of  the  finite  boundary  was  evident. 

With  the  present  computer  program,  soil  parameters,  load 
parameters,  and  computational  parameters,  the  length  of  computer  run 
time  required  for  each  case  is  approximately  40  minutes.  With  the 
completion  of  the  CDC  6600  computing  facility  at  WPAFB,  this  computer 
run  time  can  be  reduced  considerably  thus  permitting  the  evaluation  of 
additional  cases  to  determine  the  dependency  of  optimum  spacing  on  the 
load  and  soil  parameters. 

B.  Rolling-Tandem  Wheels 

The  stationary  vertical  pulse  load  problem  discusseo  n  the 
preceding  sections  docs  not  simulate  the  tandem  wheel  configuration 
because  of  the  time-lapse  effect  between  leading  and  trailing  wheels. 
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In  progress,  however,  is  the  development  of  a  computer  program  to 
simulate  multiwheel  effects  in  rolling  tandem  wheel  configurations.  For 
the  same  reasons  given  in  developing  the  twin  wheel  simulation  format, 
the  actual  three-dimensional  problem  for  rolling  tandem  wheels  was 
approximated  by  the  corresponding  two-dimensional  plane  strain  problem. 

The  Two-Dimensional  Plane  Strain  Approximation 

The  two-dimensional  plane  strain  approximation  to  the  problem 
for  rolling  tandem  wheels  has  the  same  pressure  boundary  condition  and 
soil  medium  assumption  as  for  the  stationary  vertical  load  pulse  problem. 
The  differences  are  as  follows: 

a)  The  two  infinite  pressure  strips  move  in  the  x-direction  at  the 
aircraft  ground  velocity,  as  shown  in  Figure  24a. 

b)  The  pressure-time  dependence  is  different.  As  the  pressure 
strip  moves,  the  uniform  pressure  is  increased  gradually  to  a  pressure 
equal  to  the  landing  gear  vertical  load  divided  by  the  total  contact  area, 
as  shown  in  Figure  24b.  The  manner  and  rates  at  which  the  pressure  is 
increased  initially  depends  on  the  stability  of  the  numerical  calculation 
and  the  rate  with  which  steady  state  is  achieved. 

c)  There  is  no  longer  a  plane  of  symmetry  between  the  pressure 
strip,  so  two  applied  pressure  strips  must  be  considered. 

d)  The  width  of  each  applied  pressure  strip  is  equal  to  the  tire 
footprint  length,  1,  and  the  wheel  spacing  is  expressed  as  mD,  where 

m  is  a  fraction  of  tire  outside  diameter  and  D  is  the  tire  outside  diameter 
see  Figure  24). 

Method  of  Solution- Lumped  Parameter  Iteration  Method 

The  method  of  solution  of  the  two-dimensional  problem  for 
rolling  tandem  wheels  is  similar  to  the  method  of  solution  for  the  vertical 
load  (multiwhecl)  problem.  The  method  of  prescribing  the  boundary  con¬ 
ditions,  governing  equations,  and  numerical  procedures  is  the  same. 

The  governing  equations,  presented  in  Appendix  II,  are  also  applicable 


=  Tire  Footprint  Length 
(a)  Soil  Medium  with  Moving  Strip  Pressure 


(b)  Pressure-Time  Curve 

Figure  24.  Rolling  Tandem-Wheel  Interface  Boundary  Condition 


here.  The  main  difference  was  the  boundary  condition  defined  by  the  two 
moving  pressure  strips.  The  pressure  strip  distributions  were  moved  in 
the  x-direction  at  the  aircraft  horizontal  ground  velocity,  V,  in  the  following 
manner.  The  prescribed  equivalent  pressure  exerted  on  the  mass  point 
located  at  the  extreme  left  of  each  loaded  area  was  reduced  to  zero  in  incre¬ 
mental  steps  through  the  time  h/V;  at  the  same  time,  the  equivalent 
pressure  exerted  on  the  mass  point  located  at  the  extreme  right  of  each 
loaded  area  was  increased  incrementally  from  zero  to  the  prescribed  value 
in  the  same  manner. 

Method  of  Evaluation  of  the  Multiwheel  Effects 
The  method  for  evaluating  the  multiwheel  effects  of  the  rolling 
tandem  wheel  configuration  from  the  results  of  the  two-dimensional  moving 
pressure  strips  problem  is  as  follows:  First,  the  case  with  m=0  corres¬ 
ponding  to  the  single  moving  wheel  is  run  and  the  drag,  R_.  ,  ,  is  obtained 

Single 

from  its  sinkage  (see  Equation  1).  Then,  the  cases  with  other  values  of  m 

are  run  and,  for  each  case,  the  total  drag,  R_  ,  is  obtained  from  sum- 

Tandem 

ming  the  drags  corresponding  to  the  sinkages  of  the  leading  and  trailing  pres¬ 
sure  strips.  The  multiwheel  effects  are  then  evaluated  by  studying  the  ratio, 

R 

Tandem 
P 

_ 7_ 

(  ^Single 

1  p  j 

for  various  wheel  spacing  factor,  m. 

Computer  Program  and  Test  Cases 

A  computer  program  has  been  written  based  on  the  method  of 
solutions  just  outlined.  This  program  is  now  being  used  to  evaluate  the 
tandem-tracking  cases  with  wheel  spacing  factors  of  m=®,  1.05,  1.  5,  and 
2.0.  All  these  cases  have  the  same  soil,  load,  and  computational  parameters, 
and  these  parameters  are  the  same  as  those  used  in  the  twin  wheel  cases 
conducted.  The  complete  results  of  all  the  cases,  together  with  the  listing 
of  the  final  computer  program,  will  be  presented  at  a  later  date  as  a 
Research  and  Development  computer  program  interim  report. 
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5. 


PART  4  -  Multiwheel  Verification  Tests 


Pur  pos  e 

The  purpose  of  this  test  program  was  to  study  the  variables 
that  control  the  performance  of  aircraft  tire  flotation  while  operating  in 
rnultiwheel  configurations.  The  program  provided  data  to: 

1)  confirm  the  preliminary  multiwheel  flotation  performance 
criteria  with  tires  operating  in  twin,  tandem,  and  tandem¬ 
nontracking  configurations; 

2)  use  in  the  development  of  semi- empirical  relationships  for  the 

j 

extension  of  the  single  wheeL  prediction  equations  to  multi¬ 
wheel  data;  and 

3)  aid  in  preliminary  verification  of  the  finite  element  based 
analytical  rolling  wheel  sinkage  prediction  techniques  (both 
twin  and  tandem)  currently  under  development. 

During  the  testing,  the  effects  of  twin  wheel  spacing,  tandem  wheel 
spacing  in  tracking  configuration,  and  wheel  offset  combined  with  tandem 
spacing  in  the  nontracking  configuration  were  investigated  with  the  ob¬ 
jective  of  evaluating  the  drag  modification  effects  of  the  various  spacings. 

T est  Program 

In  order  to  accomplish  the  objectives,  the  following  test 
program  was  designed  such  that  the  sinkages  and  drag  loads  would  be  in 
the  range  of  application  to  aircraft  flotation  analysis.  The  test  program, 
which  is  shown  in  Table  6  for  clay  type  soil  and  in  Table  7  for  sand  type 
soil,  was  designed  for  a  Type  III  7.  00-6,  6  PR  tire  with  a  35%  operating 
deflection.  Each  multiple  wheel  configuration  had  two  7:00-6  tires 
equally  loaded  and  all  tests  were  run  at  10  fps.  The  following  parameters 
were  measured  for  each  test: 

Vertical  Load 

Drag  Load 

Tire  Deflection 
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TABLE  6 

iVlulti'wheel  Test  Program  for  Clay  Soil 


Vertical 
As  sembly 
Load  (*) 

P 

Lateral 

Spacing 

nb 

Longitudinal 

Spacing 

mD 

1200 

1.  5b 

. 

•  1 

2.  Ob 

. 

1  • 

3.  0b 

- 

" 

4.  0b 

- 

1400 

- 

2.  on 

" 

- 

3.  on 

" 

- 

4.  on 

•• 

2.  0b 

2.  on 

1 1 

3.0b 

2.  on 

" 

3.  0b 

3.  on 

" 

2.  0b 

2.  on 

1000 

1.  5b 

. 

" 

1.  5h 

- 

" 

2.  0b 

1 1 

2.  0b 

- 

3.  0b 

- 

4.  0b 

- 

1200 

- 

1.  513 

,, 

- 

2.  on 

" 

- 

3,  on 

•« 

4.  on 

800 

1.  5b 

. 

" 

1.  5b 

II 

2.  Ob 

_ 

1000 

- 

1.  513 

TABLE  7 

Multiwheel  Test  Program  for  Sand  Soil 


Bogie 

Configuration 


Twin 


Tandem  -  T racking 


Tandem-  Nontracking 


I  Twin 

i 


Tandem  Tracking 


Single 


Soil 

Strength 

Cl  (psi) 
avg'r 


Vertical 

Assembly  Lateral 
Load  (#)  Spacing 
P  nb 


Longitudinal 

Spacing 

mD 


1600 

1.  5b 

1.  75b 

- 

2.  Cb 

- 

3.  0b 

- 

3.  0b 

- 

4.  0b 

- 

1400 

- 

1.  5D 

" 

- 

2.  on 

" 

- 

3.  0D 

- 

4.  0D 

1400 

2b 

1.  5D 

2b 

2.  0D 

2b 

3.  0D 

3b 

1.  5D 

,r 

3b 

2.0D 

’* 

3b 

3.  0D 

800 

1.  5b 

- 

2.0b 

3.  0b 

3.  0b 

. 

’  * 

4.  0b 

600 

- 

1.  5D  1 

- 

2.  0D 

- 

3.  on 

4.  on 

O 

o 

- 

l.  5n 

600 

- 

3.  on 

700 

- 

- 

400 

- 

40-45 


25-30 


40-45 
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Soil  Strength 

Soil  Sinkage 

Wheel  Velocity  (horizontal) 

The  lateral  and  longitudinal  spacings  listed  in  Tables  6  and  7  were 
selected  to  permit  a  thorough  evaluation  of  the  multiwheel  phenomena 
based  on  the  preliminary  multiwheel  flotation  criteria. 

Test  Equipment 

All  multiple  rolling  wheel  tests  in  the  program  were  conducted 
attheU.S.  Army  Waterways  Experiment  Station,  Vicksburg,  Mississippi, 
at  the  model  wheel  facility  of  the  Mobility  and  Environmental  Division 
between  the  dates  of  August  25,  1970  and  November  12,  1970.  WES 
modified  their  dynamometer  to  conduct  the  multiwheel  tests,  and  provided 
adequate  professional  and  technician  support  for  conducting  the  tests. 

The  modified  dynamometer  shown  in  Figure  25  used  the  existing  single 
wheel  dynamometer  as  a  frame  to  mount  a  lever  arm  to  hold  the  multi¬ 
wheel  carriage.  The  lever  arm  pivot  is  located  at  the  front  of  the  existing 
single  wheel  carriage  and  applies  the  load  to  the  geometric  center  of  the 
multiwheel  carriage  where  the  carriage  pivots  freely.  The  test  data  is 
monitored  using  some  existing  instrumentation  on  the  single  wheel  dyna¬ 
mometer  and  some  new  instrumentation  mounted  on  the  multiwhecl 
carriage.  Figure  25  shows  the  carriage  in  a  twin  wheel  configuration 
with  the  multiwheel  cari;age  fastened  to  the  lever  arm  so  that  it  cannot 
rotate  relative  to  the  lever  arm.  Figure  26  shows  the  carriage  in  a 
tendem  configuration  with  the  carriage  free  to  rotate  so  as  to  distribute 
the  load  equally  to  each  iir>.  A  complete  description  of  the  basic  Test 
Facility  is  available  in  a  previous  WES  report' 

Test  Tire 

In  order  to  improve  the  correlation  between  the  multiwhecl 
results  and  previous  single  wheel  testing,  the  tire  size  was  chosen  to  be 
the  7:00-6,  6  PR  Type  III  tire  which  was  used  in  the  single  wheel  testing 
program.  Table  8  gives  the  tire  geometry  data  for  the  two  7:00-6  tires 
used  in  this  program. 
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C  onfiguration 


REPRODUCIBLE 

... '  €t ' 


Soil  Tests  and  Preparation 

The  two  soil  types  chosen  for  these  multiwheel  tests  were  buckshot 

clay  and  mortar  sand,  both  of  which  were  used  in  the  previous  single  wheel 

testing.  Two  purposes  were  fulfilled  by  the  soil  tests  conducted,  first  to 

insure  an  accurate  description  of  the  test  soil  and  its  uniformity;  second  to 

allow  possible  correlation  to  other  tire-soil  interaction  theories  by  collecting 

as  much  information  concerning  the  soil  as  possible.  The  soil  tests  that 

were  conducted  are  moisture  and  density  determination,  mobility  cone 

penetration  resistance  (Cl),  and  the  California  Bearing  Ratio  (CBR).  A 

complete  description  of  each  soil  test  is  given  in  Appendix  IV.  The  summary 

of  the  moisture-density  determinations  are  in  Table  9.  The  summary  of  the 

correlation  data  taken  to  relate  CRB  and  Cl  is  presented  in  Table  10. 

avg 

TAT_E 9 

MOISTURE-DENSITY  DAI  A 
Design  Average  Conditions 


Soil  Strength 

Dr}  Density 

Moisture  Content 

Soil  Type 

Clavg  (0  to  6") 

?d(pcF) 

W,  (%) 

Buckshot  Clay 

20  _  .. 

75.  6 

42 

Mortar  Sand 

40 

101.6 

less  than  1.  0 

TABLE  10 

TEST  BED  SOIL  CONDITIONS 


Scil  Type 

Design 

Soil  Strength 
Clavg  (0"  to  6") 

Cl 

avg 

0"  to  6" 

CBR 

CBR  0.  1” 

Buckshot  C’ay 

20 

22.  3 

0.  70 

Mortar  Sand 

40 

40.  4 

2.  36 
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Undrained  triaxial  tests  were  also  run  on  undisturbed  samples  of  clay  and 
bulk  samples  of  sand  taken  from  the  soil  carts  to  determine  the  cohesion 
and  friction  components  of  the  soil  strength.  These  test  results,  which  are 
described  in  detail  in  Appendix  IV,  are  summarised  below  in  Table  11. 


TABLE  11 

TRIAXIAL  TEST  RESULTS 


Design 

Triaxial  Results 

Moisture 

Soil  Strength 

c 

4> 

Content 

Soil  Type 

Clavg  (0-6") 

(psf) 

deg 

W% 

Buckshot  Clay 

20 

1.79 

- 

42 

Mortar  Sand 

40 

36.  5 

less  than  1.  0 

Test  Results  -  Buckshot  Clay 

The  finalized  test  results  for  the  28  tests  run  in  buckshot  clav 
are  presented  in  Table  12.  The  data  presented  represents  the  parameters 
that  were  recorded  for  permanent  records  in  both  computer  tape  format 
and  oscillograph  charts.  The  following  comment?  clarify  the  column 
headings  and  notations.  Three  tests  (numbers  4,  6,  and  7)  are  considered 
to  be  void  due  to  excessive  sinkage  and  basic  instability  of  the  test  carriage 
during  the  test  run.  The  spacings  listed  in  the  third  column  represent 
actual  test  spacing®  based  on  the  tire's  nominal  diameter  (D)  of  18"  and 
its  nominal  width  (b)  of  6.8".  The  soil  strength  (Cl^^j  represents  the 
average  of  at  least  five  before  traffic  tests,  and  is  given  in  terms  of  the 
average  penetration  resistance  over  the  first  six  inches  of  soil  profile  in 
psi.  Soil  strength  measurements  taken  in  the  tire  ruts  alter  one  pass  are 
given  in  Appendix  IV.  Note  that  the  east-west,  and  front-rear  designations 
have  been  used  consistently  to  designate  the  relative  position  of  each  tire 
to  the  lest  section. 
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The  rut  depth  measurements  given  represent  the  difference 
between  the  original  soil  profile  and  the  soil  profile  after  the  passage  of 
the  test  carriage.  The  sinkage  values  presented  were  calculated  as  a 
given  percentage  of  the  rut  depth  based  on  evaluation  of  the  rebound 
measured  by  deflection  pins  buried  in  the  tire  paths.  A  description  of  the 
deflection  pins  is  given  in  Appendix  IV.  Soil  uniformity  measurements 
and  an  outline  of  the  basic  test  procedure  is  also  presented  in  Appendix  IV. 

Using  the  analysis  techniques  established  during  development 
of  preliminary  multiwheel  theories,  the  verification  data  was  reduced 
and  plotted  as  shown  in  Figures  27,  28,  and  29,  to  indicate  the  multiwheel 
performance  relative  to  spacing  for  the  twin,  tandem-tracking  and  tandem¬ 
nontracking  configurations.  This  multiple  wheel  performance  is  evaluated 
by  a  comparison  of  drag  ratios  as  shown  in  Equation  11. 

mm  =  (R/P)s  '  (r/p)m  <n> 

where: 

(R/P)  =  single  wheel  drag  ratio 

(R/P).  =  multiple  wheel  drag  ratio 
M 

Mw  =  multiple  wheel  drag  modifier 

M 

In  order  to  plot  this  data,  a  value  of  the  single  wheel  drag  ratio  for  each 
specific  condition  was  required.  The  single  wheel  drag  ratio,  (R/P)  , 
used  in  this  comparison,  was  arrived  at  by  using  the  data  collected  and 
verified  during  the  single  wheel  test  program^  ^  and  not  the  value  that 
would  be  obtained  by  use  of  the  single  wheel  nomograph^  The  reason 
for  this  selection  is  the  close  control  and  accuracy  associated  with  the 
single  wheel  tests  previously  conducted  which  used  the  same  soil  and 
test  tires,  whereas  the  nomograph  represents  an  averaged  result  of 
numerous  tests  with  and  without  test  controls. 

In  general,  the  multiwheel  data  for  twin  and  tandem-tracking 
tests  grouped  by  the  magnitude  of  the  tire  sinkage.  Thus,  it  appears 
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re  load 


Figure  27.  Twin  Tire  Drag  Performance  on  Clay-Multiwheel  Tests 
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Figure  28.  Tandem-Tracking  Tire  Drag  Performance  on  Clay-Multiwheel  Tests 
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that  a  series  of  sinkage  ranges  exist  in  which  the  value  will  be  constant 
for  a  given  sinkage  at  a  given  spacing.  The  sinkage  ranges  are  approximately 
from  0  to  0.5",  0.5"  to  1.0",  and  1.0"  to  4.0",  called  low,  intermediate 
and  high  sinkage  ranges,  respectively.  Reference  to  Figure  27  shows  two 
data  points  plotted  far  below  the  rest  of  the  data.  No  explanation  for  this 
phenomena  exists  although  those  two  tests  were  run  in  a  clay  section  different 
from  all  other  c’  ly  tests.  Note  also  that  in  Figure  29,  the  tandem-nontracking 
data  for  both  sand  and  clay  have  been  plotted  together,  as  the  trends  in  this 
data  indicate  very  little  if  any  effect  of  tandem-nontracking  configurations 
for  the  spacings  used.  It  should  also  be  noted  that  due  to  limited  funding  the 
low  sinkage  tandem-nontracking  test  condition  was  not  run,  and  therefore, 
could  not  be  evaluated.  The  curves  labeled  "high"  and  "low"  sinkage  in  the 
multiple  wheel  test  data  plots  are  an  attempt  to  bracket  the  test  data  while 
reflecting  the  effects  of  different  magnitudes  of  sinkage. 

Test  Results  -  Mortar  Sand 

The  test  results  for  the  29  multiwheel  tests  run  in  the  mortar 
sand  are  presented  in  Table  13.  The  following  comments  explain  the  column 
headings  used  in  Table  13.  The  spacings  listed  in  the  third  column  repre¬ 
sent  actual  test  spacings  based  on  the  tire's  nominal  diameter  (D)  of  18" 

and  its  nominal  width  (b)  of  6.8".  The  soil  strength  (Cl  )  represents  the 

avg 

average  of  at  least  five  tests  before  traffic,  and  is  defined  as  the  average 
penetration  resistance  over  the  first  six  inches  of  soil  profile.  Soil 
strength  measurements  taken  in  the  tire  ruts  after  one  pass  are  given  in 
Appendix  IV.  The  east-west,  and  front-rear  designations  have  again  been 
used  to  designate  the  relative  position  of  each  tire  to  the  test  section.  The 
sinkage  listed  is  actually  the  measured  rut  depth,  yet  is  considered  to  be 
the  total  sinkage  composed  of  both  elastic  deformation  (very  small  in 
mortar  sand)  and  the  permanent  deformation  since  sands  exhibit  only  minimal 
rebound.  The  soil  uniformity  measurements  and  an  outline  of  the  basic  test 
procedure  is  presented  in  Appendix  IV. 
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Muitiwheel  Drag-Sinkage  Data-Sand  Soil 
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Figure  30.  Twin  Tire  Drag  Performance  on  Sand-Multiwheel  Tests 


The  sand  data  was  evaluated  based  on  drag  performance  and 
plotted  as  shown  in  Figures  30  and  3 1  to  indicate  the  multiwheel  performance 
relative  to  wheel  spacing.  The  value  of  the  single  wheel  drag  ratio  was 
obtained  from  the  data  taken  during  the  single  wheel  testing  program  as 
explained  previously.  As  for  clay,  the  results  show  that  the  data  tends  to 
group  by  the  magnitude  of  sinkage  and  the  same  three  ranges  generally  apply. 
Two  of  the  tandem-tracking  tests  were  run  on  a  much  weaker  strength  soil 
as  shown  in  Figure  31.  Some  additional  twin  wheel  performance  data  on  % 

sand  was  obtained  from  unpublished  work  by  WES.  This  data  although  some¬ 
what  scattered,  tended  to  substantiate  the  trends  shown  in  the  twin  wheel 
results  of  Figure  30. 

Summa  ry 

Although  the  multiwheel  tests  presented  here  are  not  in  them¬ 
selves  conclusive  due  to  the  limited  test  program,  the  basic  trends  shown 
by  this  data  when  combined  with  the  plate  test  results  and  the  analytical 
results  will  permit  the  development  of  design  and  operational  criteria  for 
multiwheel  configurations.  The  test  results  generally  indicate  that  the 
spacinp  of  twin  wheels  is  critical  in  minimizing  drag  for  multiwheel  con¬ 
figurations  but  that  tandem-tracking  spacing  and  tandem  nontracking  spacing 
are  less  critical  in  influencing  multiwheel  drag.  The  results  of  the  multi¬ 
wheel  analysis  and  its  application  to  multiwheel  flotation  criteria  are  pre¬ 
sented  in  Section  IV. 

6.  Summary  of  Multiwheel  Drag  and  Sinkage  Performance 

The  results  of  the  review  of  existing  twin  tire  test  data  (Figure  3), 
the  Twin  Plate  Tests  (Figure  13),  the  analytical  results  (Figure  19),  and 
the  Multiwheel  (Twin)  Verification  Tests  (Figures  27  and  30)  generally  show 
that  for  the  proper  twin  wheel  spacing,  a  reduction  in  the  drag  ratio  (and 
sinkage  ratio)  over  that  for  single  wheels  will  occur  for  low  to  moderate 
sinkages  of  twin  tires  operating  on  both  sand  and  clay.  For  high  sinkages, 
the  twin  tire  drag  ratio  will  generally  exceed  the  single  wheel  drag  ratio  for  twin 
tires  in  clay  soil  but  remain  less  than  the  single  wheel  drag  ratio  for  sand  soils. 
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Although  the  analytical  results  for  tandem-tracking  tires  arc  not  yet 
available,  the  results  of  the  review  of  existing  data  (Figure  4),  and  the 
Multiwheel  (Tandem)  Verification  Tests  (Figures  28  and  31)  show  a  slightly 
improved  performance  (-15%)  over  single  wheel  performance  in  sand  soil 
but  an  increase  in  drag  ratio  (  =  10%)  on  clay  soil  for  most  tandem-tracking 
tire  spacings. 

The  results  of  the  analytical  and  experimental  studies  conducted  to 
date  on  multiple  wheel  tirc/soil  interaction  effects  permit  the  following 
preliminary  conclusions. 

TWIN  WHEEL  -  CLAY  SOIL  -  Slight  reduction  (*10%)  in  the  rolling  drag 
ratio  (R/P)  in  clay  soils  for  twin  tires  is  gained  over  that  for  isolated  single 
wheels  by  proper  spacing  of  the  tires  for  low  sinkage  situations.  For  high 
sinkage  situations,  increases  occur  (*15%)  in  rolling  drag  ratio  over  that 
for  single  isolated  wheels.  The  proper  twin  wheel  spacing  lies  between 
2.0  b  and  3.  0  b  where  b  is  the  tire  width. 

TWIN  WHEEL  -  SAND  SOIL  -  Reductions  in  the  rolling  drag  ratio  (R/P) 
for  the  proper  spacing  are  possible  in  sand  type  soil  for  twin  wheel  con¬ 
figurations  in  comparison  to  single  wheels  for  all  sinkages.  Potential 
reductions  in  the  drag  ratio  of  10%  to  20%  are  evident.  The  proper  twin 
wheel  spacing  lies  between  1.5  b  and  2.5  b. 

TANDEM- TRACKING  -  CLAY  SOIL  -  The  operation  of  tandem-tracking 
tires  in  clay  soil  will  yield  some  reduction  (*10%)  in  the  drag  ratio  (R/P) 
for  the  proper  tire  spacing  in  low  sinkage  situation  when  compared  to 
single  wheel  performance  but  higher  drag  ratios  (  =  15%)  for  high  sinkage 
conditions.  Available  evidence  suggests  a  desirable  tandem-tracking  spacing 
of  1.5  D  to  2. 0  D. 

TANDEM -TRAC KING  -  SAND  SOIL  -  The  operation  of  tandem-tracking  tires 
in  sand  type  soils  will  yield  slight  benefits  (up  to  10%)  in  reducing  the  drag 
ratio  over  that  of  a  single  wheel  for  proper  twin  wheel  spacing.  The  spacing 
of  the  tandem-tracking  tires  does  not  appear  to  influence  the  magnitude  of 
the  drag  ratio. 
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TANDEM-NONTRACKING  -  CLAY  AND  SAND  SOIL  -  Although  the  number 
of  tandem-nontracking  tests  conducted  in  clay  and  sand  type  soil  was  limited, 
this  preliminary  data  suggests  that  tandem-nontracking  configurations  have 
little  influence  on  sinkage  and  drag  interaction  between  aircraft  tires  for 
the  normally  encountered  tire  spacings. 

As  indicated  above,  the  results  indicate  that  multiple  wheel  geo¬ 
metric  configuration  (dual,  tandem,  spacing,  etc.)  is  an  important  parameter 
in  determining  the  rolling  drag  performance  of  aircraft  tires  on  soil  and 
should  receive  increased  consideration  in  aircraft  design  when  soil 
operational  performance  is  a  requirement. 


62 


SECTION  III 

BRAKING  SINKAGE  AND  DRAG  ANALYSIS 


Although  the  flotation  requirements  for  aircraft  landing  gear  design 
and  aircraft  tire  selection  should  be  based  primarily  on  the  minimization  of 
rolling  wheel  drag  (minimization  of  sinkage)  during  takeoff  operations,  the 
braked  tires  on  soil  effects  are  an  important  factor  in  defining  the  maximum 
drag  loads  on  landing  gears  and  in  determining  aircraft  stopping  distances 
on  unprepared  runways.  The  braked  tire  on  soil  problem  is  considerably 
different  from  rolling  tire  action  since  a  resultant  shear  force  will  exist  at 
the  tire/soil  interface  for  braked  tires.  The  major  variables  present  for 
braked  tire  conditions,  in  addition  to  those  considered  for  rolling  tires,  are 
tire  slip  and  the  shearing  stresses  at  the  interface.  Large  increases  in 
drag  and  in  some  cases  sinkage  occur  for  braked  tires  over  that  for  rolling 
tires. 

The  major  portion  of  research  efforts  on  studying  the  tire/soil 
response  for  slip  conditions  has  been  conducted  on  powered  wheels.  Very 
little  theoretical  or  exoerimental  work  has  been  accomplished  for  braked 
tires  at  speeds  and  conditions  peculiar  to  aircraft  operations.  Those 
studies'  ’  which  have  been  accomplished  to  date  have  been  restricted 
to  field  efforts  directed  at  determining  experimental  braking  coefficients 
for  limited  conditions.  The  present  state  of  the  art  of  soil  mechanics 
precludes  the  development  of  a  theoretical  braking  analysis  due  to  the  lack 
of  an  adequate  equation  of  state  describing  the  stress-strain-time  charac¬ 
teristics  of  soil.  For  these  reasons,  a  dual  approach  is  being  taken  to 
developing  a  more  complete  understanding  of  braked  tire/soil  interaction 
phenomena.  The  semi-analytical  braking  theory  described  below  is  an 
attempt  to  provide  a  useable  theory,  partially  verified  by  existing  experi¬ 
mental  data,  for  defining  braked  tire/soil  interaction.  The  finite  element 
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braking  analysis  currently  in  progress  is  also  summarized  in  this  section 
and  is  an  attempt  to  take  a  long  range  theoretical  look  at  the  braking  problem. 

1.  Semi-Analytical  Braking  Theory 

v  The  braked  tire  condition  creates  a  force-slip  condition  where  the 

degree  of  .braking  is  defined  by  the  tire  slip,  which  is  given  as 


V 

S  =  <3—  -  i)  100 
a 

where 


(12) 


S  =  percent  tire  slip 

V  =  peripheral  speed  of  wheel,  and 
w 

V  =  horizontal  velocity  of  the  axle. 

The’  slip  is  considered  to  be  positive  when  the  wheel  is  in  a  powered  con¬ 
dition  and  negative  under  braked  conditions.  A  fully  braked  (locked) 
wheel  is  a  condition  of  100  percent  negative  slip. 

a  u  r  ,  (2,7,18-23)  .  .  .  , 

A  number  of  references  were  used  in  developing 

the  braking  drag  expression.  The  significant  variables  included  in  the 

analysis  were: 

Loads 

P  =  vertical  tire  load 
Rg  =  braked  tire  drag  force 


Tire  Geometry 

b  =  tire  section  width 

D=  tire  outside  diameter 

l  =  tire  footprint  length  on  rigid  surface 


Soil  Strength 

c  =  soil  cohesion 

4)  =  soil  friction  angle 

Cl  =  average  cone  index  over  0"  to  6" 
avg 

=  slope  of  cone  index  versus  depth  profile  in  sand 
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Other 


S  =  percent  tire  slip 
V  =  aircraft  horizontal  velocity 

The  braking  drag  resistance  as  shown  in  Equation  (13)  is  based  on 
the  premise  that  the  total  drag  resistance  can  be  treated  as  the  sum  of  two 


components: 


rb  =  rr  +  rt 


the  horizontal  soil  resistance  to  forward  motion  exclusive  of  soil  shear 

resistence  during  braking  (R  ) ,  plus  the  horizontal  component  (R  )  of  the  net 

shearing  force  resistance  (T)  between  the  tire  and  the  soil  as  a  function  of 

slip.  Figure  32  shows  the  loading  and  interface  conditions  for  a  braked  tire. 

(19) 

Justification  exists  for  the  use  of  Equation  (13)  in  that  form.  Schuring 

has  stated  that  the  rolling  resistance  (R  )  is  independent  of  slip  and  that  the 

K  (21) 

coefficient  of  friction  (related  to  R^,)  is  a  function  of  slip  alone.  McCrae 

has  also  assumed  that  for  driven  wheels  at  any  slip,  the  term  R^  can  be 

✓ 

determined  from  a  rolling  resistance  formula.  Both  R  and  R  are 

R  I 

influenced  by  the  magnitude  of  sinkage.  R  is  influenced  through  use  of  the 

R 

rolling  drag  equation  (also  see  Equation  1) 

-~=f(Z/D)  (14) 

while  R^,  varies  with  sinkage  since  the  size  of  the  shear  stress  contact  area 

at  the  tire/soil  interface  is  influenced  by  the  magnitude  of  the  sinkage.  A 

( 18 ) 

review  of  some  existing  information'  '  for  braked  tires  on  soil  indicates 
that  sinkages  increase  only  slightly  for  all  percents  slip  on  clay  soil  while 
sinkages  increase  markedly  for  braked  tires  on  cohesionless  soils.  This 
phenomena  is  shown  very  clearly  in  Figures  33  and  34. 

R  Term,  Clay  and  Sand  Type  Soil 
R 


As  indicated  previously,  the  horizontal  .resistance  to  forward  motion 
exclusive  of  shear  resistence  is  independent  of  slip  and  can  be  approximated 


p 


R  =  forward  motion  soil  drag 
R 

P  =  vertical  soil  reaction 
s 

T  =  tangential  shear  force 

-  equivalent  plane  of  contact 


Figure  32.  Braked  Tire/Soil  Interface  Conditions 
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%  Negative  Slip  (S) 


Figure  33.  Influence  of  Negative  Slip  on  Sinkage,  Clay  Soil 


Braked  Sinkage 
Free  Rolling  Sinkage 


Figure  34.  Influence  of  Negative  Slip  on  Sinkage,  Sand  Soil 


from  the  rolling  drag  equation  for  cohesive  soils^  as 


-jT  =3.85  (Z/D) 


(15) 


Since  sinkages  in  clay  soil  increase  only  slightly  throughout  the  negative 

slip  range,  the  R  term  is  assumed  to  remain  essentially  constant.  This 

conclusion  is  partially  verified  as  shown  in  Figure  35  by  the  results  of 

an  experimental  effort^^  in  which  R  was  measured  separately  from  R  , 

R  T 

In  sand  type  soils,  the  R  term  will  increase  in  magnitude  since 

K 

sinkages  increase  with  increasing  negative  slip.  A  change  in  R  caused 

by  braking,  however,  should  correspond  to  the  same  change  in  drag  (R) 

for  rolling  tires  for  an  equivalent  change  in  sinkage.  A  comparative 

study^^  of  this  phenomena  is  shown  in  Figure  36  and  tends  to  substantiate 

the  use  of  the  rolling  drag  equation  for  sand.  The  rolling  drag  equation  for 
,(1)  .  . 

sand  is  given  by 


R 


R 


=  0.048  +  2.  77  (Z/D) 


(16) 


It  should  be  noted  that  R  for  sand  will  increase  throughout  the  negative 

R 

slip  range  thus 


<V 


S/0 


=  <v 


+  A'R. 


s=o 


(17) 


where 

ARr  =  f  (AZ) 

R^  Term,  Clay  and  Sand  Type  Soil 

The  tangential  force  (T)  at  the  tire/soil  interface  is  the  integral  of 

the  shear  stresses  over  an  equivalent  plane  of  contact  as  shown  in  Figure  32. 

The  horizontal  component  of  T  is  the  term  R^  and  contributes  significantly 

to  the  total  braked  drag  resistance  R  .  If  Coulombs  law  for  determining 

B 
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R  =  Horizontal  Resistance  to  Forward 

Motion  Exclusive  of  Shear  Resistance 


(VS£Q 

(Rr)s=o 


ro  Negative  Slip 


Figure  3$.  Influence  of  Negative  Slip  on  R^,  Clay  Soil 


AR  =  Change  i-  Rolling  Drag 
ARp  =  Change  in  Forward  Motion 

Resistence  Exclusive  of  Shear 
Resistence 


braked 


rolling 


Figure  36.  Comparison  of  Braked  Rp  With  Rolling  Tire 
Drag  for  Increasing  Sinkage,  Sand  Soil 


shear  force  at  the  interface  is  valid,  the  shearing  resistance  at  the  inter¬ 
face  could  be  computed  as  a  function  of  slip  as 

T  =  (c  +  o  tan  <)>)  n(S)  (18) 

where 

t  =  shearing  stress 
o  =  effective  normal  stress 
n(S)  =  nonlinear  function  which  varies  with  slip 

For  the  rolling  tire  case  for  which  the  braking  torque,  B',  is  zero, 
experimental  evidence  indicates  that  although  there  are  shearing  stresses 
at  the  tire/soil  interface,  the  net  shearing  force,  T,  is  zero,  and  the 
resultant  of  the  vertical  and  horizontal  soil  reaction  passes  through  the 
center  of  the  tire.  As  a  braking  torque  is  applied  continuously  to  the  tire, 
the  negative  slip  increases  creating  additional  shear  stresses  at  the  inter¬ 
face.  These  shearing  stresses  result  in  a  net  shearing  force  as  given  by 

T  =  Az(c  +  otan«J>)u(S)  (19) 

where 

A  =  area  of  equivalent  plane  of  contact  at  sinkage  Z. 

The  resultant  force  at  the  interface  no  longer  passes  through  the  tire 
center,  but  acts  through  an  eccentric  point  creating  a  moment  equivalent 
to  the  braking  torque,  B',  for  a  steady  state  braking  condition.  The 
horizontal  component  of  T  is  given  by 

Rt  =  A^(c  +  otan$)  n(S)  cos  8  (20) 

Reference  to  Figure  32  indicates  that  the  angle  6  defining  the  equivalent 
plane  of  contact  is  given  as 

e  =  90*  -  [sir'Ml  -~)  +  1/2 cos-1  (1  -  ^y)]  (21) 
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The  use  of  Equation  (20)  necessitates  the  accurate  description  of 
the  tire/soil  contact  area  as  a  function  of  sinkage.  For  certain  clay  type 

soils  this  A  term  would  have  to  include  sidewall  contact  areas.  Using 

^  (3  21) 

simplifying  assumptions,  expressions  have  been  developed  by  others'  ’ 

which  analytically  define  the  footprint  length  (i  )  and  footprint  width  (b  ) 

i  Z  Z 

from  which  could  be  determined  if  sidewall  effects  are  ignored. 

Analysis  of  these  expressions  indicates  that  for  tires  operating  in  their 

normal  deflection  range  (±l07o  of  their  normal  deflection),  that  the  contact 

area  can  be  written  as  a  functional  expression  in  the  following  manner 


f(Az,  Z,  D) 


"gfr  =  f  (Z/D)  (23) 

(18) 

The  results  of  a  controlled  braking  study  conducted  by  WES'  for 
the  Air  Force  were  used  to  examine  Equations  (20)  and  (23)  leading  to  the 
development  of  an  R^,  expression  for  cohesive  and  cohesionless  type  soils. 
These  tests  were  conducted  using  two  different  tires  (D  =  14"  and  28-1/2"), 
one  tire  deflection  (D  =  257o),  various  soil  strengths  and  vertical  loads  in 
both  sand  and  clay,  and  negative  slips  ranging  from  0  to  1007o. 

R^.,  Cohesive  Soil 

For  cohesive  soils,  the  tan  $  is  zero  in  Equation  (20)  and  the 

cohesion,  c,  can  be  replaced  by  the  average  cone  index,  Cl  (or  Cl  ). 

avg 

Equations  (20)  and  (23)  can  then  be  combined  to  give  the  dimensionless 
expression 


Cl  •  Ds  (Z/D)' 


-  s{S)  cos  ? 


The  clay  test  results'  were  plotted  based  on  Equation  (24)  and  Figure  37 
shows  that  a  reasonable  convergence  of  the  data  exists  for  the  case  n  =  1/2. 
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R  values  were  determined  by  subtracting  a  calculated  R  value  based  on 
the  magnitude  of  the  instantaneous  sinkage,  Z  (see  Equation  15).  The  data 
used  in  developing  Figure  37  is  summarized  in  Appendix  V.  Based  on 
Figure  37  then,  the  braked  tire  horizontal  shearing  force  at  the  interface 
can  be  approximated  by 

Rt  =  0.  11  CI-D2(Z/D)1/2u(S)  (25) 

R^,  Cohesionless  Soil 


In  cohesionless  soils,  the  cohesive  term,  c,  in  Equation  (20)  is  zero. 
Replacing  ~  in  Equation  (20)  by  a  =  P/A,  where  A  is  the  rigid  surface  tire 
contact  area,  Equation  (20)  and  (23)  can  be  combined  as 


R. 


D2(Z/D) 


n 


=  a  tan0  •  u(S)  cos  p 


(26) 


(181 

Initial  analysis  of  the  sand  data1  '  indicated  that  the  R_  term  in  sand  was 

T 

only  a  weak  function  of  Z/D  (n  is  very  small)  and  consequently  within 
acceptable  accuracy  limits,  Equation  (26)  can  be  written  as 


-j~=-  =  a  tan</»-u(S)cos  p  (27) 

Observations  of  braked  tire  performance  in  soil  shows  that  considerable 
sand  flow  takes  place  at  the  tire  soil  interface.  This  sand  disturbance  and 
flow  very  likely  causes  the  shearing  strength  to  be  determined  by  some 
large  deformation  equilibrium  void  ratio  condition  in  sand  rather  than  the 
initial  soil  shear  strength.  If  this  phenomena  is  true,  analysis  of  Equation  (27) 
should  show  R^,  to  be  approximately  constant  for  all  values  of  initial  shear 
strength  for  sand  i'or  all  other  factors  remaining  constant.  Figure  38  which 
includes  soil  strength  changes  up  to  an  order  of  magnitude  of  three,  shows 
that  the  initial  soil  strength  does  not  have  a  significant  influence  on  the 
magnitude  of  R^  for  braked  tires  on  -and.  Analysis  of  Figure  38  in  con¬ 
junction  with  Equation  (27)  leads  to  the  expression 


R 


T 


(28) 


f  «- 


CI-D*(Z/D)  '  cosO 


o 


Cl 


avg 


Figure  38.  Horizontal  Shear  Component  Evaluation,  Sand  Soil 
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As  for  the  clay  data  analysis,  R  values  were  determined  by  subtracting  a 
calculated  R  value  based  on  the  magnitude  of  the  instantaneous  sinkage,  Z 
(see  Equation  16).  The  data  used  in  developing  Figure  38  is  summarized  in 
Appendix  V. 

Variance  of  R  With  Slip 

As  the  negative  slip  increases,  the  R^,  value  in  both  sand  and  clay 

increases  in  magnitude.  This  increase  in  R^,  with  increasing  slip  occurs  in 

a  nonlinear  manner  as  shown  in  Figure  39.  Although  Figure  39  shows  some 

differences  in  the  growth  of  R  /R  with  slip  for  sand  and  clay,  as  a  first 

T  Tmax 

approximation,  the  quantity  f_i(S)  in  Equation  (25)  and  (28)  can  be  written  as 

S  1  /3 

M(s)  *  rI  “  <Too)  <*’> 

T 

max 


Variance  of  Sinkage  With  Slip,  Cohesionless  Soil 

Due  to  the  large  increases  in  sinkage  for  braked  tires  on  sand  type 

soils,  these  sinkages  must  be  known  since  it  significantly  influences  the 

resulting  drag,  Rfi.  By  comparing  the  sinkage  at  100%  slip  (zs^q)  to  the 

sinkage  at  zero  slip  (Zg_Q)  for  the  test  data  from  the  WES  study(*®),  Figure  40 

was  developed  which  shows  that  the  initial  load-strength  ratio  (o /Cl  )  is 

avg 

critical  in  determining  this  ratio.  Although  considerable  research  remains 
to  be  accomplished  in  order  to  accurately  define  the  growth  of  sinkage  for 
braked  tires  on  sand,  Figure  40  can  be  used  to  estimate  the  sinkage  at  S  =  100% 
in  the  following  manner: 

1.  Calculate  the  sinkage  for  a  rolling  tire  (S=0)  using  the  previously 
developed  sinkage  prediction  equations  for  sand^.  This  sinkage 
is  Z  in  Figure  40. 

o  —  U 

2.  Determine  a/CI  for  the  braked  tire  and  use  Figure  40  to  find 

avg 

(Z„,n)  /Z_  From  this  value  the  maximum  braked  tire 

SfO  max  S=0 

sinkage  [Z_,n)  in  sand  can  be  determined. 

5*0  max 
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R  =  Horizontal  Component  of 
^  Net  Shear  Force 
Resistence 

Rt  =  Maximum  Rt 
•••max  0  A  <3 


o  *  ■ 


max  o.  < 


A  Q 


WES(I<,)  Data 

0  Cl  =  17,  clay 

□  Cl  =  iO,  clay 

O  CI=  30,  clay 

A  Cl  =  i<>,  cla\ 

9  G,  =  1  i.  6.  sand 
b 

■  G  =  6.  6,  sand 
b 

A  G  =  18.  0,  sand 


A  G  =  11.5,  sand 
b 


%  Negative  Slip  (S) 

Figure  39.  Increase  in  Braked  Shear  Force  at  the 
Tire/Soil  Interface,  Sand  and  Clay 


^S^O  max 


WES,10)  Data 
O  9.00-14 
A  4.00-7 


,  )  =  Maximum  Braked  Sinkage 

S^0  max 

Z  .  =  Free  Rolling  Sinkage 


Load  -  Strength  Ratio, 


Figure  40.  Variation  of  Braked  Sinkage  Ratio  With 
Load-Strength  Ratio,  Sand  Soil 


Braking  Equations  -  Summary 


Cohesive  Soil 


R 


_B 

P 


R 


R 


R. 


„  Z  0.  11  CI*DS  Z  1/2 
3*  85  ^  - p - <d  )  U(S) 


(30) 


Cohesionless  Soil 


0.048  +  2.  77  (“ )  u(S) 


(31) 


where  |j(S)  is  given  by  Equation  (29). 

2.  Comparative  Study  Using  Semi-Analytical  Braking  Theory 

Limited  experimental  efforts  have  been  sponsored  by  the  Air  Force 
to  conduct  braked  tire  on  soil  tests  on  a  full  scale  basis.  Douglas^, 
Boeing^,  and  Lockheed^  have  conducted  these  test  programs.  The 
results  of  these  tests  provided  an  opportunity  to  compare  the  braking 
drag  analysis  equations  developed  from  controlled  laboratory  testing  to 
the  results  of  these  field  operations.  Equations  (30)  and  (31)  together 
with  other  information  was  used  in  developing  a  computer  program  to 
conduct  this  comparative  study.  The  results  of  the  study  are  summarized 
below. 


Lockheed  Tests 

(2) 

The  Lockheed  Corporation  in  1968,  conducted  a  series  of 
rolling  and  braked  tire  tests  at  the  NASA/Langley  test  track  using  a  29  inch 
diameter,  Type  III  tire.  Braked  tire  tests  were  conducted  on  both  buckshot 
clay  and  a  beach  sand  common  in  the  Langley  area.  The  soil  strength, 
vertical  load,  and  tire  deflection  for  each  test  are  summarized  in  Table  14 
and  Table  15  for  clay  and  sand  respectively.  Braked  drag  predictions  in 
the  clay  soil  (Buckshot  clay)  were  based  on  sinkage  by  adding  to  the 
measured  rut  depth  a  rebound  deflection.  This  rebound  deflection  was 
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Braking  Test  Variables,  Lockheed  Clay  Tests 
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estimated  based  on  data  from  previous  tests  conducted  at  WES^  (also  in 
Buckshot  clay)  where  both  rebound  and  rut  depth  were  measured.  Table  16 
presents  the  measured  rut  depth  and  the  corrected  rut  depth  (sinkage) 
using  the  results  of  the  Phase  II  study^  while  Table  17  gives  the  measured 
rut  depths  for  the  sand  tests.  Very  little  rebound  occurs  in  sand  type  soils 
and  consequently  the  rut  depths  were  used  as  sinkages. 

The  results  of  the  braked  drag  predictions  are  summarized  in 
Table  18  for  the  clay  soil  tests  and  in  Table  19  for  the  sand  soil  tests. 
Comparisons  between  the  predicted  braking  drag  using  Equations  (30)  and 
(31)  and  the  measured  braking  drag  (results  of  Table  18  and  19)  are  shown 
in  Figures  41  and  42  for  clay  and  sand  respectively.  In  both  Figures, 

±20%  lines  are  shown  in  order  to  provide  an  estimate  of  the  accuracy  of 
the  prediction  technique.  Considering  the  limitations  which  exist  in 
preparing  uniform  soil  test  beds  and  the  accuracy  with  which  tests  measure¬ 
ments  on  soil  can  be  taken,  the  results  as  shown  in  Figures  41  and  42  are 
quite  favorable. 

Boeing  Tests 

(7) 

The  Boeing  Company  in  1964,  carried  out  a  series  of  aircraft 

field  tests  including  both  roiling  and  braking  tire  conditions  using  a  Boeing 

367-80  (KC-135  prototype/  aircraft.  The  tests  were  conducted  on  a  lean 

clay  lake  bed  located  at  Harpers  Lake,  California.  The  aircraft  was 

equipped  with  a  high  Potation  landing  gear  (basic  707  type  bogie)  with 

46  x  16,  Type  VII  tires.  Only  one  braking  test  was  conducted  on  soil 

and  a  summary  of  the  test  variables  and  the  resulting  measured  braking 

coefficient  (R  /P)  as  well  as  the  predicted  braking  coefficient  are  shown 
B 

in  Table  20. 

Douglas  Tests 

Douglas  Aircraft  Company  also  ran  a  series  of  aircraft  on 
soil  field  tests  in  1964-65  using  a  C-5A  prototype  landing  gear  configuration. 


78 


TABLE  16 

Braking  Rut  and  Sinkage  Data, 
Lockheed  Clay  Tests 
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Braking  Drag  Comparison,  Clay  Soil 
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Measured  Drag  vs.  Predicted  Drag*  Sand  Soil 


Table  20.  Comparative  Braking  Study  Results--Boeing  Data 
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Uning  clay  prediction  equations. 


The  test  program  which  was  also  conducted  at  Harpers  Lake,  California, 
included  rolling  aircraft  tests  and  a  limited  number  of  braking  tests  on  a 
soil  runway.  Although  tire  data  was  not  specified,  the  tires  were  very 
likely  40  x  14,  Type  VII  tires.  The  braking  tests  were  conducted  on  a  lean 
clay  with  a  CBR  rating  of  16.  Analysis  of  soil  strength  data,  however, 
indicates  that  the  rated  CBR  16  runway  had  CBR  values  ranging  from  8 
to  34  which  certainly  influenced  the  resulting  measured  braking  coefficients. 
A  summary  of  the  test  variables  and  the  resulting  comparisons  using  the 
braking  prediction  equations  is  given  in  Table  21. 

The  results  as  given  in  Table  21  show  the  predicted  braked  drag 
ratio  to  be  significantly  higher  than  the  measured.  These  differences  are 
likely  related  to  the  magnitude  of  the  sinkages  used  in  the  prediction 
equations.  Due  to  the  variability  in  soil  strength  at  the  test  site,  the 
estimated  sinkages  shown  in  Table  21  could  be  considerably  in  error. 

Based  on  the  results  of  the  comparative  study,  the  braked  drag 
prediction  equations  developed  herein  can  be  used  as  a  first  approximation 
for  determining  aircraft  braking  perfoimance  on  soil  runways.  These 
braking  equations  which  were  developed  for  sand  and  clay  type  soil  can 
also  be  used  to  bracket  the  braked  drag  ratios  for  intermediate  type  soils. 

3.  Finite  Element  Braking  Analysis 

As  indicated  previously,  the  interaction  of  aircraft  tires  and  soil 
runways  while  an  aircraft  is  being  braked  is  a  complex  phenomenon.  In 
actuality,  braking  does  not  occur  independently  of  other  equally  compli¬ 
cated  effects  such  as  the  rolling  action  of  the  tires,  the  speed  of  the  air¬ 
craft,  and  turning.  Additionally,  the  proximate  tires  interact,  the 
behavior  is  three-dimensional,  the  composition  and  material  properties 
of  soils  are  usually  uncertain,  etc. 

The  purpose  of  this  section  is  to  outline  the  numerical  approach 
currently  under  development  for  the  prediction  of  the  linkage  of  aircraft 
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tires  during  braking.  In  order  to  study  the  effects  of  braking  on  sinkage, 
braking  is  isolated  from  other  effects  such  as  those  mentioned  above. 

It  may  be  appropriate  at  a  later  time  to  attempt  to  include  as  many  factors 
as  possible  in  a  single  analysis;  however,  the  approach  outlined  below  has 
been  adopted  to  obtain  an  initial,  "first-cut",  grasp  of  the  braking 
phenomenon. 

The  remainder  of  this  section  is  presented  in  four  principal  parts. 

(1)  The  first  part  contains  a  description  of  the  specific  problem  which  will 
be  considered.  The  problem  definition  includes;  a  listing  of  the  physical, 
geometric,  and  material  assumptions  which  are  made  (e.  g.  ,  isolation  of 
braking  effects,  plane  strain,  etc.);  a  description  of  the  simulated  loading; 
a  description  of  the  region  of  half-space  in  which  a  solution  will  be 
obtained;  and  a  discussion  of  the  assumed  boundary  conditions  at  the 
extremities  of  the  region  of  solution.  (?.)  The  second  part  contains  a  brief 
description  of  the  mathematical  model  which  is  intended  to  represent  the 
behavior  of  an  idealized  sample  of  soil  subjected  to  simulated  weight  and 
brake  loadings.  This  discussion  contains:  the  proposed  finite  element; 
the  assumed  displacement  state  for  the  finite  element;  the  resulting  stiffness 
matrix;  and  a  possible  finite  element  modeling  of  the  soil  sample.  (3)  The 
third  part,  contains  a  description  of  the  solution  procedure  which  is  pro¬ 
posed  (load-incrementation  technique).  (4)  The  last  part  contains  a  brief 
summary  of  rhe  current  progress  of  the  development  of  the  finite  element 
braking  analysis. 

PART  I  -  Problem  Definition 

The  idealized  problem  considered  is  defined  below  by  listing  the 
assumptions  and  discussing  the  loading,  region  of  solution,  and  boundary 
conditions. 

A  s  s  umptions 

-  A  single  wheel  is  in  contact  with  the  sample  of  soil  under 

8^ 


consideration. 


-  Only  the  effects  of  vertical  loading  and  braking  are  con¬ 
sidered  (the  aircraft  is  not  turning,  the  soil  surface  is  smooth,  etc.). 

-  The  deformation  of  the  soil  material  under  the  loading 
considered  results  in  a  state  of  plane  strain.  This  is  equivalent  to 
assuming  that  the  tire  is  infinite  in  width.  Sidewall  shear  is,  therefore, 
neglected;  however,  this  is  most  likely  a  secondary  effect.  Consequently, 
the  problem  is  reduced  to  two  dimensions. 

-  The  loading  is  assumed  to  be  applied  slowly  so  that  accel¬ 
eration  effects  of  the  soil  can  be  neglected. 

Loading 

Figure  43  shows  the  portion  of  the  half-plane  which  is  loaded 

by  a  uniform  vertical  pressure,  p  ,  and  a  uniform  shear  distribution,  p  . 

n  rs 

The  indicated  loading  is  intended  to  represent  a  very  simplified  simulation 
of  the  effect  of  an  aircraft  tire  during  braking.  The  shear  distribution, 
p  ,  is  a  function  of  the  contact  area  between  the  tire  and  the  soil,  tjie 

s  r 

cohesion,  the  friction  angle,  and  the  percent  slip  between  the  tire  and  the 
soil.  The  loading  is  assumed  to  remain  constant  during  the  deformation. 

Region  of  Solution 

Thus  far  the  problem  has  been  reduced  to  a  loaded  half-space. 

In  order  to  effect  a  solution  by  numerical  means,  however,  the  extent 
of  the  region  affected  by  the  load  must  be  restricted  to  be  finite.  Figure  44 
shows  two  possible  approaches.  In  Figure  44a,  a  rectangular  area  has  been 
designated  as  the  region  of  solution  while  a  semicircular  area  is  indicated 
in  Figure  44b.  In  each  case,  an  attempt  is  made  to  select  the  dimensions 
of  the  region  so  that  the  load  has  negligible  influence  on  the  displacements 
at  the  extremities  of  the  region. 

Boundary  Conditions 

-  Under  the  load  the  normal  stress  is  equal  to  the  uniform 

vertical  pressure  (p  )  and  the  shear  stress  is  equal  to  the  uniform  shear 

n 

distribution  (p^)  (Figure  44). 
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SOIL 


Xj-  Tire  Footprint  Length 

p  =  Uniform  Vertical  Pressure 
n 

p  =  Uniform  Shear  Distribution 
rs 


Figure  43.  Simulated  Loading  During  Braking 
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soil  surface 


Figure  44.  Region  of  Solution  for  Finite  Element 
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several  tire 
diameters 


-  The  shear  and  normal  stresses  are  zero  on  the  remainder 
of  the  soil  surface. 

-  The  displacements  are  zero  on  the  artificial  boundaries 
which  limit  the  extent  of  the  soil  medium  (Figure  44). 

PART  2  -  Mathematical  Model 

The  technical  literature  pertaining  to  the  elastic-plastic  analysis  of 

structures  by  the  finite  element  method  has  been  reviewed  and  evaluated. 

A  comparative  study  of  the  most  pertinent  sources  appearing  in  the 

(24) 

literature  through  1967  is  given  by  Marcal  .  It  has  been  determined 
that  the  best  existing  approach  to  the  prediction  of  elastic-plastic  response 
is  the  technique  presented  in  Reference  24  and  explained  in  greater  detail 
in  Reference  25  .  Essentially,  the  method  involves  applying  the  loading 
in  small  increments  and  performing  a  linearized  direct  stiffness  analysis 
at  each  stage  of  the  load  incrementation.  The  mathematical  mode)  is  the 
incremental  potential  energy  of  an  assembled  set  of  finite  elements.  The 
discretized  stiffness  equations  are  based  on  assumed  finite  element  dis¬ 
placement  approximations  since  a  potential  energy  formulation  is  basically 
a  displacement  formulation.  The  primary  disadvantage  of  Marcal' s  tech¬ 
nique  is  that  the  predicted  stress  state  is  discontinuous  between  the 
individual  elements;  this  situation  is  common  to  all  displacement  approaches 
however.  One  valuable  attribute  of  the  method  is  that  it  is  equally  appli¬ 
cable  to  both  perfectly  plastic  and  work  hardening  materials. 

It  is  believed  that  an  improved  method  for  solving  elastic-plastic 
problems  can  be  developed  by  using  an  incremental  Reissner  energy^^ 
expression  rather  than  an  incremental  potential  energy  formulation  rs  the 
mathematical  model.  The  Reissner  energy  treats  both  displacements 
and  stresses  as  primary  variables.  Therefore,  the  discretized  stiffness 
equations  will  be  based  on  assumed  finite  element  stress  and  displacement 
states  rather  than  just  assumed  element  displacement  states.  Conseq  lently, 
it  will  be  possible  to  obtain  continuous  stresses  as  well  a3  continuous 
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displacements  with  a  Reissner  energy  formulation.  Some  of  the  principal 
characteristics  of  the  method  to  be  used  for  predicting  the  sinkage  of 
braked  aircraft  tires  are  outlined  below. 


Finite  Element 

Two  different  finite  elements  are  considered  as  possible  candidates 
for  the  modeling  of  the  region  of  solution:  a  rectangular  element  (Figure  45) 
and  a  triangular  element  (Figure  46).  The  displacement  and  stress  approx¬ 
imations  within  each  of  the  two  types  of  finite  elements  are  taken  as  bilinear 
functions  of  the  spatial  variables  in  the  form: 


( r) 

(x,  y)  =  Cj  +  c^x  +  c3y  +  c^y 

(32) 

(x,  y)  =  dj  +  d£x  +  d3y 

(33) 

where  the  superscripts  r  and  t  refer  to  the  rectangular  and  triangular 
elements,  respectively.  The  function  /(x,  y)  represents  any  of  the  dis¬ 
placement  and  stress  variables  u,  v,  o  ,  a  ,  c  .  Both  forms  of  the 

x  y  xy 

approximating  functions  given  in  Equations  (32)  and  (33)  contain  sufficient 
generality  to  insure  continuous  displacements  and  stresses  across  the 
common  boundaries  of  interfacing  elements. 


As  in  the  normal  application  of  the  finite  element  method,  the 
displacement  and  stress  approximations  can  be  expressed  in  terms  of  the 
nodal  values  by  evaluating  Equations  (32)  and  (33)  at  the  nodes  (Figure  4? 
or  46)  and  solving  for  the  constants  c^,  i=l,  2,  3,  4  or  d^,  i=l,  2,  3  depending 
on  which  element  type  is  utilized.  For  example,  the  displacement  and 
stress  states  of  a  rectangular  finite  element  can  be  put  in  the  form 


=  1  P(x,  y)  1(6} 


(34) 


90 


where  f  6  1  is  a  vector  which  contains  the  nodal  values  of  the  stresses 

and  displacements  and  (P(x,  y) }  is  a  vector  of  bi-linear  functions  of  x  and 

y.  The  vector  of  nodal  variables,  (6)  contains  20-elements  since  there 

are  4-nodal  values  (in  the  case  of  the  rectangular  element)  for  each  of  the 

6-variables  u,  v,  a  ,  a  ,  and  a 

x  y  xy 

Stiffness  Matrix 


The  Reissner  energy  for  a  typical  finite  element  in  plane  strain 
which  is  situated  within  the  interior  of  an  assemblage  of  finite  elements 
can  be  written  as 

r  r  1 

Tf_  “II  [3u  +av  +  a  (u  +v  )  -  -7r[(l-v2)(c2  +c 2) 

R  J  -'a  X  ,  X  v  .  v  xv  .v  .  v'  2P  '  ' '  v  v  ' 


xy 


,  x 


-  2v  ( 1  +  v)  a  a  +  2(  1  +  v)  T  2  ] }  dx  dy  (35) 

x  y  xy 

Substituting  the  assumed  displacement  and  stress  states  given  by 
Equations  (32)  and  (33)  into  the  Reissner  energy  expression  and  performing 
the  indicated  integrations  over  the  element  area  results  in  a  quadratic  form 
in  the  nodal  variables 

irR  =  j  JsJ  [ k ]  { 6  }  (36) 

where  [k]  is  the  stiffness  matrix  associated  with  the  Reissner  energy. 

The  total  Reissner  energy  of  an  assembled  set  of  finite  elements 
can  be  written  as 

7Tr  =  \  LAJ  -  [_aJ(p}  (37) 

where  ‘.A;  is  a  vector  containing  the  independent,  non-imposed  nodal 
variables,  K  :  is  the  assembled  stiffness  matrix,  and  {P}  is  a  vector  of 
imposed  nodal  displacements  or  stresses.  The  first  term  in  Equation  (37) 
is  the  total  Reissner  energy  of  an  assemblage  of  elements  and  the  second 
term  is  interpreted  as  the  external  work  associated  with  imposed  stresses 
and  displacements  (boundary  conditions). 


Invoking  the  stationary  Reissner  energy  principal,  Equation  (37) 
implies  that 

[K]{a}={P]  (38) 

from  which  the  unknown  nodal  displacements  and  stresses  can  be 
determined  by 

{A}  =  [K]"  1  {P }  (39) 

Typical  Finite  Element  Modeling 

Figure  47  shows  a  typical  finite  element  modeling  of  the  region 
of  solution  of  Figure  44b  using  triangular  elements.  In  this  case,  the 
material  in  the  vicinity  of  the  loading  has  been  modeled  with  small 
elements  while  increasingly  larger  elements  are  used  toward  the  extremity 
of  the  region.  Typical  modeling  of  the  region  in  Figure  44a  is  not  shown 
here  since  it  is  more  or  less  obvious.  Note  that  the  uniform  normal  and 
shear  distributions  have  been  converted  into  equivalent  concentrated  node 
forces.  The  eventual  solution  to  the  problem  may  or  may  not  require  a 
more  refined  modeling  than  the  one  shown. 

PART  3  -  Solution 

If  the  soil  material  were  completely  elastic  in  nature,  the  solution 
could  be  obtained  by  applying  Equation  (39)  a  single  time.  The  material, 
however,  is  assumed  to  be  elastic-plastic.  Therefore,  if  the  stresses  in 
an  element  arc  of  such  a  magnitude  that  plastic  deformation  occurs  the 
stiffness  matrix  for  that  element  depends  on  the  stress  state  and  is  no 
longer  constant.  Clearly,  this  makes  the  assembled  stiffness  matrix 
nonconstant  also.  Consequently,  the  solution  procedure  will  be  to  express 
the  Reissner  energy  in  terms  of  incremental  stresses  and  displacements 
and  derive  an  "incremental  stiffness  matrix"  which  includes  the  effects 
of  plastic  flow  as  well  as  initial  stresses  and  displacements.  The  loading 
applied  to  the  soil  surface  will  then  he  applied  in  small  increments  and  at 


0  3 


each  stage  of  the  loading  the  stiffness  matrices  will  be  updated  for  those 
elements  which  have  yielded.  The  total  deformation  and  stress  states 
caused  by  the  applied  loading  will  be  obtained  by  summing  the  incremental 
stresses  and  displacements  obtained  during  each  load  increment.  The 
residual  displacement  state  (sinkage)  and  stress  state  which  remains 
after  unloading  will  be  obtained  by  the  solution  described  above  minus 
a  completely  elastic  solution. 

PART  4  -  Progress -to-Date 

-  An  appropriate  incremental  Reissner  energy  expression  for 
two-dimensional,  plane  strain,  clastic-plastic  problems  has  been  derived. 

-  Finite  element  stiffness  relations  based  on  the  Reissner  energy 
principle  have  been  derived  fer  rectangular  elements  with  bilinear, 
assumed  stress  and  displacement  states.  The  bilinear  variation  is 
sufficient  to  assure  interelement  stress  and  displacement  continuity  and 
solution  accuracy  can  be  controlled  by  varying  the  number  and  size  of 
the  finite  elements  used  to  model  the  region  of  solution. 

-  A  test  computer  program  has  been  nearly  completed.  No 
example  cases  have  been  treated  as  yet. 

4.  Braked  Wheel  Verification  Tests 

The  simplifying  assumptions  made  in  the  development  of  the 
braking  analysis  equations  (Equations  30  and  31)  as  well  as  the  results 
of  the  comparative  study  certainly  points  to  the  need  for  additional 
braked  tire  on  soil  performance  data  from  both  controlled  laboratory 
testing  and  full  scale  field  tests.  Currently  scheduled  as  a  p«rt  of  this 
continuing  research  effort  are  a  series  of  Braked  Wheel  Verification  Tests. 
These  tests  will  begin  in  early  February,  1971  and  the  results  of  these 
tests  wili  very  likely  lead  to  some  modification  in  Equations  (30)  and  (31). 
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SECTION  IV 

ADDITIONAL  STUDIES  IN  TIR E/SOIL  INTERACTION 

1 .  Single  Wheel  Comparative  Study 

Under  Phase  II  of  the  research  program,  a  single  wheel  simulation 
computer  program  was  developed  for  predicting  sinkages  in  sand  and  clay 
soil  of  moving  aircraft  tires.  The  results  from  the  development  of  this 
program  which  used  a  lumped  parameter  approach  are  available  elsewhere^  ^ 
Subsequent  to  this  computer  program  development,  a  number  of  comparisons 
were  made  between  predicted  sinkages  as  determined  from  the  computer 
program  usage  and  the  experimentally  determined  sinkages  from  the 
Single  Wheel  Verification  Tests  of  Phase  11^  \  The  results  of  this  com¬ 
parative  study  are  given  below. 

The  results  of  the  comparisons  are  shown  in  Table  22  for  Tests 
No.  8,  0,  10,  12,  and  19  in  clay  soil  of  the  Single  Wheel  Verification 
Tests^l  Loading  and  soil  parameters  used  in  the  sinkage  prediction  were 
taken  from  data  developed  from  these  tests^  \  Reference  to  Table  22 
indicates  that  in  3  cases  the  predicted  sinkages  were  within  16%  of  the 
measured  sinkages.  Differences  in  the  other  two  cases  compared  ranged 
up  to  approximately  40%. 

In  order  to  study  the  influence  of  varying  such  parameters  as 
cohesion.  Young's  modulus,  and  duration  of  pulse  while  holding  other 
parameters  constant,  an  additional  series  of  runs  were  made  using  the 
single  wheel  simulation  computer  program.  A  summary  of  the  additional 
cases  is  shown  in  Table  23. 

The  results  of  varying  the  cohesion  and  Young's  modulus  are 
shown  in  Figure  48.  In  this  graph,  the  sinkage  characteristic  (defined  as 
the  instantaneous  sinkage  divided  by  the  footprint  length)  is  plotted  against 
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TABLE  22 


SINKAGE  PREDICTION  BY  STATIONARY  PULSE  LOAD  COMPUTER 
PROGRAM  COMPARED  WITH  EXPERIMENTAL  RESULTS 

(CLAY  SOILS) 


Test 

No. 

Vertical 

Load 

(lbs) 

P 

Peak 

Pressure 

p(psi) 

P 

rmax 

E 

(psi) 

c 

(psf) 

$ 

(deg) 

Experi¬ 
mental 
Sinkage 
Z,  (inch) 

Analytical 

Sinkage 

Z,  (inch) 

High  Strength  Clay 

8 

2019 

21.  3 

575 

560 

1 

0.  39 

0.  33 

9 

1000 

i7.  9 

575 

560 

1 

0.  30 

0.  17 

10 

1494 

25.  1 

575 

560 

1 

0.  52 

0.68 

Low  Strength  Clay 

12 

995 

10.  5 

178 

230 

1 

0.  91 

1.  03 

19 

503 

9.6 

178 

230 

1 

0.  60 

0.68 
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TABLE  23 


ADDITIONAL  SINKAGE  PREDICTION  CASES  -  CLAY  SOIL 


T  est 

No. 

Soil 

Type 

Verti¬ 

cal 

Load 

P(lbs) 

Con¬ 

tact 

Area 

Min2) 

Peak 

Pres¬ 
sure  p 
FVnax(psi) 

E 

(psi) 

c 

(psf) 

Time 

§  Duration 

(deg)  (sec) 

8 

Clay 

2019 

94.  8 

21.  3 

575 

450 

1 

0.  1 

1  1 

II 

II 

1  1 

It 

II 

500 

II 

II 

rt 

II 

If 

If 

It 

It 

560 

II 

1  1 

ti 

tr 

If 

IT 

1 1 

425 

II 

II 

/  II 

9 

ti 

1000 

55.  8 

17.  9 

575 

450 

II 

It 

tf 

if 

II 

1 1 

It 

1  1 

500 

II 

II 

n 

1 1 

M 

1 1 

1 1 

1  1 

560 

II 

1  1 

10 

1 1 

1494 

59.  5 

25.  1 

II 

450 

II 

II 

H 

1 1 

II 

II 

It 

1 1 

560 

It 

II 

11 

ti 

2016 

57.  9 

34.  8 

1  1 

450 

It 

II 

12 

ti 

995 

95.  0 

10.  5 

178 

230 

II 

II 

H 

n 

II 

II 

1 1 

II 

II 

1  1 

0.05 

II 

1 1 

It 

1 1 

1 1 

II 

It 

II 

0.01 

II 

1 1 

II 

1 1 

It 

1! 

255 

It 

0.  1 

1 1 

II 

It 

It 

II 

280 

II 

M 

II 

1 1 

1 1 

1  1 

II 

1  * 

350 

II 

M 

19 

1 1 

503 

52.  3 

9.6 

1 1 

230 

1  1 

M 

II 

1 1 

II 

I  1 

II 

M 

255 

1  1 

II 

II 

1 1 

1  1 

1  1 

1 1 

1 1 

280 

II 

1 1 

II 

1 1 

II 

1  1 

II 

1 1 

350 

II 

It 
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the  load- st r ength  ratio  (defined  as  the  average  tire  contact  pressure,  a, 

divided  by  the  average  cone  index  of  the  soil).  Shown  in  the  same  graph 

are  experimental  data  previously  compiled  from  the  literature.  The 

analytically  predicted  results  show  a  trend  of  dependence  on  the  dimensionless 

parameter,  o/E,  where  E  is  Young's  modulus  of  elasticity.  For  lower 

values  of  load  strength  ratio  (<0.  5),  the  analytically  predicted  sinkage 

ratios  are  quite  close  to  the  experimental  results.  This  indicates  that 

the  computer  program  simulation  is  satisfactory.  But  for  higher  load 

strength  ratios,  the  predicted  sinkages  are  larger  than  their  corresponding 

experimental  results,  and  the  trends  are  indicated  by  the  solid  line  curves. 

This  trend  of  higher  sinkage  occurs  because  a  non-strain  hardening  soil 

model  was  used  and  for  a/Cl  >0.  5,  the  pressures  were  many  times 

avg 

greater  than  the  limiting  yield  strength. 

The  results  of  varying  the  time  duration  of  the  load  pulse  are 
shown  in  Figure  49.  Ir,  this  graph,  the  time  duration  of  the  load  pulse  was 
converted  to  horizontal  ground  velocity  of  the  aircraft  based  on  the  relation¬ 
ship  given  in  Reference  13.  The  sinkage  decreases  as  the  velocity  in¬ 
creases.  The  known  increase  of  sinkage  in  the  Region  III  velocity  range 
does  not  occur  in  this  instance  because  the  stationary  load  sinkage  pre¬ 
diction  computer  program  cannot  account  for  the  plowing  and  hydroplaning 
effects  of  the  high  speed  forward  motion. 


SECTION  V 


MULTIWHEEL  FLOTATION  CRITERIA  AND  RESEARCH  PLANNING 

I.  Multiwheel  Flotation  Criteria 

The  multiwheel  flotation  criteria  described  herein  is  an  extension 
of  the  previously  developed  single  wheel  flotation  criteria^.  The 
multiwheel  criteria  permits  the  evaluation  of  aircraft  flotation  performance 
rather  than  single  tire  performance.  It  also  will  permit  aircraft  designers 
to  determine  optimum  landing  gear  configuration  (twin,  tandem,  spacing, 
etc.)  for  aircraft  leading  to  drag  minimization. 

The  basis  of  the  multiple  wheel  drag  criteria  is  the  previously 
defined  multiple  wheel  dr*.  «  modifier  (see  Section  II,  Part  4)  which  leads 
to  the  definition  of  the  multi-wheel  drag  ratio  as 

<fi>m  ■  ,40> 

where 

(— )  =  single  wheel  drag  ratio 

P  u 

M  =  multiple  wheel  drag  modifier 
M 

By  using  the  following  definition  of  terms, 

N  =  total  number  of  wheels  per  landing  gear 

N  =  number  of  wheels  (of  N)  in  a  twin  situation 
n 

N  =  number  of  wheels  (of  N)  in  a  tandem-tracking  situation 
m 

=  number  of  wheels  (of  N)  in  a  tandem-nontracking  situation 

K  =  drag  modifier  for  twin  wheel  situation  (see  Figure  50) 
n 

K  =  drag  modifier  for  tandem-tracking  situation  (see  Figure  51) 
m 

=  drag  modifier  for  tandem-nontracking  situation  (see 
Figure  52) 

the  drag  variations  due  to  multiple  wheel  configurations  are  given  by 
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N  (R-K  R)  =  drag  increase  or  decrease  due  to  twin  wheel  effects 
n  n 


N  (R-K  R)  =  drag  increase  or  decrease  due  to  tandem-tracking 
mm 

wheel  effects 


N^R-K^R)  =  drag  increase  or  decrease  due  to  tandem- nontracking 
wheel  effects. 


where 

R  =  single  wheel  drag. 

The  resulting  multiwheel  drag  per  tire  is  then  given  as 


^multiwheel  per  tire  ^ 


N  N  N 


(41) 


and  the  multiwbe'el  drag  ratio  becomes 


(— )  =  (— ) 
'p'm  'p's 


N  N  N 

l-j*rr~(l-K  )+-~— (1-K  )+-f(l-K). 
1  N  n  N  m  N  r  J 


The  multiple  wheel  drag  modified  becomes  then 


(— ) 
'p'm 

(-) 

'p's 


=  mm  = 


N 


N 


1HTf<UKn>+Tf(UK«> 


N  1 

+Tf<UKr>J 


(42) 


(43) 


The  results  of  the  multiwheel  drag  analysis  presented  in  Section  II  were 

used  to  develop  the  drag  modifier  relationships,  (1-K  ),  (1-K  ),  and 

n  m 

(1-K^)  shown  in  Figures  50,  51,  and  52.  The  curves  shown  in  Figures  50, 
51,  and  52,  reflect  the  results  of  all  the  work  conducted  to  date  on 
multiple  wheel  interaction  effects  and  represent  an  attempt  to  bracket  these 
results  while  still  trying  to  indicate  h;gh  and  low  sinkage  magnitudes.  High 
sinkages  are  defined  as  greater  than  approximately  one  inch  while  low 
sinkages  generally  are  less  than  one-half  inch.  In  order  to  facilitate  the 
determination  of  M  ,  Equation  43  was  put  into  nomographic  form  as  shown 
in  Figure  53. 


► 


The  current  multiple  wheel  drag  criteria  is  defined  for  the  Region  II 
velocity  range.  Extensions  of  this  criteria  to  Regions  I  and  III  will  be 
made  as  analytical  and  experimental  results  become  available.  Based  on 
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Twin  Drag  Modifier,  (1- 


c 
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-0.4 
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1.5  2.0  3.0  4.0 

Twin  Spacing  Parameter,  n 
(a) 


1.  5  2.0  3.0  4.0  5.0 

Twin  Spacing  Parameter,  n 
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Figure  50.  Drag  Modifiers  for  Twin  Wheel  Situations 


Tandem-Tracking  Drag  Modifier,  (1- 
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Kf-H I 


Nont  ra 


Figure  52.  Drag  Modifiers  for  Tandem-Nontracking  Situations 
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Multiple  Wheel  Modifier  Nomogra 


available  evidence,  it  is  unlikely  that  the  relative  flotation  capacity  as 
defined  for  Region  II  will  change  when  considering  Regions  I  and  III 
(although  the  absolute  drag  ratios  will  change). 

2 .  Multiple  Wheel  Flotation  Criteria  Example 

In  order  to  illustrate  the  use  of  the  multiple  wheel  flotation  index, 
(FI)^,  an  example  was  prepared  using  the  C-5A  aircraft.  In  order  to 
arrive  at  (FI)  ,  the  single  wheel  flotation  index,  (FI)  ,  as  determined  from 

A/I  o 

Figure  54  is  required. 

C-5A  FLOTATION  RATING 


Given  the  C-5A  aircraft,  calculate  the  Multiwheel  Drag  Ratio, 
(R/P)^,  for  the  normal  landing  weight  and  tire  deflection  on  a  CBR  8 
(250  Cl)  clay  soil. 


Known  Information 

-  Clay  Soil,  Cl  =250 

avg 

-  C-5A  landing  gear  configuration  (Figure  55) 

-  Tire  size  is  49  x  17-26  PR  Type  VII 

-  Load  distribution  to  Main  Gear  (max)  is  94.2%, 

.  641,200  x  0.942  - - - 

P(per  tire)  = - ~ -  =  25,  200  lbs. 

-  Normal  tire  deflection  is  35% 

-  Normal  landing  weight  is  641,200  lbs. 


Solution 

A.  Using  Reference 
be  calculated  for  the  main  gear. 

Tire  diameter,  D 
Tire  width,  b 
Tire  deflection,  d 
Tire  print  length,  l 
Tire  contact  area,  A 
Tire  load,  P 


27,  all  tire  data  and  parameters  can 

=  48. 4  in. 

=  16.7  in. 

=  4. 36  in. 

=  13.86  in. 

=  240  sq.  in. 

=  25.  2  kips 
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Figure  55.  C-5A  Landing  Gear  Configuration 


B.  Using  the  nomograph  in  Figure  54,  calculate  the  single 
wheel  performance: 


Calculate  P/CI 


avg 


25200 

250 


=  101 


b. 


Intersect  240  sq.  in.  (A  )  with  P/CI  =  101,  which 
^  1  avg 

gives  a/ Cl  =  0.  42 

avg 


c.  Move  horizontally  to  intersect  the  cohesive  (clay)  line 

-2 

and  down  vertically  to  get  Z/ 1  =  5.  Ox  10 

-2 

d.  Intersect  Z/l  =  5  x  10  with  l  -  13.  9  to  get  Z  =  0.  7" 

e.  Intersect  Z  =  0.  7"  with  D  =  48.  4"  to  get 
(FI)S  =  (R/P)s  =  0.065 

C.  Using  Figures  50,  51,  52  and  53,  calculate  the  Multiwheel 


Modifier,  M 


M* 


a.  Calculate  the  spacing  parameters  m,  n,  and  r. 

As  the  multiwheel  criteria  is  limited  in  its  scope, 
some  simplifying  assumptions  must  be  made  to  rate 
all  aircraft  landing  gear  configurations.  The  C-5A 
main  landing  gear  consists  of  two  identical  units  of 
12  tires  each  as  indicated  in  Figure  55.  The  M 

M 

for  either  unit  is  identical  and  the  calculations  will 
be  based  on  the  following  assumptions. 

ASSUMPTIONS 


-  Assume  that  all  eight  center  tires  in  the  indicated 

unit  (shaded  in  Figure  55)  act  as  tandem-tracking 

tires  at  S  =  65". 
m 

-  As  sume  that  all  12  tires  in  the  indicated  unit  act 

as  twin  tires  at  S  =  42"  even  though  the  actual 
n 

S^  ranges  from  34"  to  53". 

-  Interaction  effects  due  to  tandem-nontracking  tires 
are  negligible  (see  Figure  52). 


Then: 
S 


m  = 


m 


65- 


48.  4- 


*  1.4 


42' 


n  = 


16.  7- 


'  2.  5 
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b.  Determine  the  Drag  Modifiers: 

-  From  Figure  50,  for  clay  soil,  move  vertically 
from  n=2.  5  to  a  point  between  the  two  lines.  The 
upper  line  represents  the  relation  for  sinkages 

of  1"  or  greater,  the  lower  line  represents  sinkages 

of  0.  5"  or  less.  Therefore,  for  Z  =  0.  7"  (see 

above)  the  point  should  be  »l/2  the  distance  between 

the  lines.  Now  move  horizontally  to  read 

(1-K  )  =  0.0. 
n 

-  From  Figure  51,  for  clay  soil,  move  vertically 

from  m=l.  4  to  a  point  between  the  two  lines  keeping 

in  mind  the  relative  values  of  the  sinkages.  Move 

horizontally  from  the  extrapolated  value  to  read 

(1-K  )  =  -0.  10. 
m 

c.  Determine  the  wheel  configuration  parameters 


N  12 


1.0 


N 

~~( 1  -K  )  =  -0.067 
N  m  - 

Total  -0.067 

M  =  1  -  (-0.067)  =  1.06V 
M 

Determine  the  Multiwheel  Drag  Ratio 

(R  /P)  =M  (R  /P)  =  1.  067x0.  065=0.  069 

MM  o 
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Note  that  this  value  is  slightly  different  from  that  which 
would  be  calculated  if  the  exact  equations  were  used  versus 
the  general  equations  used  to  construct  the  single  wheel 
nomograph'”. 

3.  Current  Aircraft  Flotation  Ratings 

The  procedure  as  detailed  above  was  used  to  evaluate  the  currently 
used  cargo  type  military  aircraft  as  to  their  relative  flotation  capacity. 

The  ratings  are  based  on  the  Multiwheel  Drag  Ratio,  (R/P)  ,  using  the 

specific  conditions  listed  in  Figure  56  which  also  shows  the  resulting  ratings 
in  a  chart  format.  The  most  commonly  used  aircraft  type  tires  have  been 
previously  rated  for  flotation  capacity  and  their  ratings  can  be  found  in 
Reference  1. 

4.  Summary  of  Flotation  Criteria 

A  summary  of  the  current  state  of  the  art  in  aircraft  flotation 
criteria  is  presented  in  Table  24.  This  table  is  an  attempt  to  present  in 
a  summary  fashion  all  of  the  current  flotation  criteria  as  developed  by 
both  UDRI  and  others. 

5.  Research  Planning 

As  a  part  of  the  continuing  research  program  by  the  Air  Force, 
research  planning  in  knowledge  deficient  areas  has  been  conducted  as  a 
part  of  the  current  effort.  The  planning  efforts  to  date  have  been  in  the 
areas  of  braking,  multipass,  and  speed.  The  results  of  these  planning 
efforts  are  currently  under  review  by  AFFDL/FEM. 

b.  Current  Research  Results  and  Reports 

In  order  to  provide  an  up  to  date  overview  of  the  current  research 
work  being  conducted  in  the  areas  of  aircraft  flotation  and  vehicle 
mobility,  the  University  has  established  a  research  library.  Efforts 
arc  continually  made  to  obtain  the  results  of  past  and  on-going  research 
programs  from  governmental  and  private  agencies.  The  library  system 
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The  aircraft  listed  are 
rated  based  on  <  ata 
taken  from  References 
27,  28,  and  29. 
using  the  conditions 
listed  below: 

Soil  Strengths^  CBR 

(25°  CI)  1281 
Normal  Landing  Wt.'  ^7) 

Normal  Tire  Deflection 

Clay  Type  Soil 


Figure  56.  Cargo  Type  Aircraft  Flotation  Ratings 
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currently  has  over  1000  documents  on  file.  Due  to  the  increasing  volume 
of  technical  reports  and  literature  reviewed  and  stored  in  the  library 
system,  it  became  necessary  to  expedite  the  indexing  and  retrieval 
process  through  an  information  system.  The  system  selected  is  referred 
to  as  the  KWIC  index  system  and  its  development  is  available  through 
services, provided  by  W right-Patter son  Air  Force  Base.  The  KWIC  index 
is  a  machine  generated  printed  listing  of  documents  which  uses  the 
keywords  in  the  document's  title  as  the  retrieval  key.  The  KWIC  index 
displays  on  continuous  computer  printout  sheets  (which  are  bound  in  book 
form)  each  significant  word  of  the  title  alphabetically  in  the  center  of  the 
page.  The  KWIC  index  also  identifies  the  author  on  the  same  line  as  the 
title.  All  flotation  related  research  articles  are  then  stored  alphabetically 
by  the  author  names.  Those  articles  having  no  identified  author  are 
numbered  consecutively  and  located  by  this  number  since  it  also  appears 
in  the  KWIC  index. 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS  FOR  RESEARCH 


1.  Conclusions 

The  results  of  the  landing  gear/soil  interaction  research  effort 
have  shown  that: 

1.  The  effects  of  geometric  configuration  (twin,  tandem,  etc.) 
are  a  significant  factor  in  defining  the  multiple  wheel  drag 
ratio,  (R/P)m. 

2.  Results  from  the  plate  tests,  analytical  solution,  and  the 
multiple  wheel  tests  are  in  general  agreement  as  to  load  and 
sinkage  interaction  effects. 

3.  Proper  twin  spacing  is  more  critical  than  tandem  spacing  in 
minimizing  rolling  drag. 

4.  Large  increases  in  sinkage  occur  for  braked  tires  on  sand 

while  only  slight  increases  in  sinkage  occur  for  braked  tires 
on  clay.  1 

5.  The  braking  prediction  equations  were  partially  verified  and 
are  suitable  to  use  on  a  preliminary  basis  for  predicting  braked 
tire  on  soil  drag  ratios. 

6.  The  analytical  sinkage  prediction  equation  is  a  reliable  pre¬ 
diction  techn.que  f  r  clays  in  the  lower  sinkage  range.  High 
sinkage  sit  rations  cannot  be  predicted  accurately  because  of 
the  lack  of  strain  hardening  in  the  current  soils  model. 

7.  The  single  wheel  drag  ratio,  (FI)<j,  has  been  modified  to  a 
multiple  wheel  drag  ratio,  (FI)j^.  This  multiple  wheel  drag 
ratio  can  now  be  used  to  determine  the  flotation  performance 
of  aircraft  during  takeoff  operations. 

2.  Recommendations  for  Research 

Immediate  efforts  for  landing  gear/ soil  interaction  and  flotation 
criteria  research  should  he  directed  at  the  following  areas. 


1.  High  speed  (Region  III)  sinkage  and  drag  performance  for 
rolling  tires. 

Z.  Additional  verification  of  the  braking  sinkage  and  drag  analysis 
equations . 

3.  Incorporation  of  strain  hardening  in  the  currently  used  soils 
model. 

Long  range  flotation  research  efforts  should  include: 

1.  Multiple  pass  (operations)  effects  as  related  to  runway  deteriora¬ 
tion  and  operational  capability. 

Z.  Landing  gear  loads  in  turning  operations. 

3.  Effects  of  roughness  on  the  flotation  performance  of  aircraft 
on  soil  runways. 
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TWIN  PLATE- VERTICAL  LOAD  TEST  PROGRAM, 
SOIL  TESTS  PREPARATION, 
PROCEDURES,  AND  UNIFORMITY  TESTS 
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Soil  Classification 


The  two  soils  selected  for  testing  were  sand  (mixture  of  mortar 
sand  and  Yuma  sand)  and  buckshot  clay.  The  grain  size  distribution  for 
the  sand  is  given  in  Figure  1  ,  .  m’  ..he  grain  size  distribution  and  limits 
properties  for  the  buckshot  clay  are  given  in  Figure  58.  The  intent  of  the 
soil  placement  was  to  provide  a  uniform  soil  test  section.  The  cone 
penetrometer  could  then  evaluate  the  soil  strength.  The  results  of  the 
moisture-density  determinations  for  the  test  sections  are  summarized 
in  Table  25. 

TABLE  25 

MOISTURE-DENSITY  DATA 


Average  Conditions 

Soil  Type 

Dry  Density 
Vd(pcf) 

Moisture  Content 
W(7o) 

Sand  ( IS  to  1 2S) 

=  95.  4 

1.  5 

Clay  1C 

87.  7 

27.  5 

2C 

82.  3 

29.0 

*C 

85.  9 

27.0 

4C 

75.  5 

29.  1 

5C 

82.  1 

27.  3 

6C 

97.  4 

27.  5 

7C 

82.  b 

28.  3 

8C 

85.  2 

29.  3 

92.0 

00 

r- 

'Nl 

Soil  Preparation 

a.  Sand 

A  standard  sand  soi 

l  placement  and  di 

•nsity  preparation 

technique  was  used  in  an  effort 

to  develop  an  acce 

ptable  uniformity  of 

!  24 


tvlc«  stand  aro  Steve  numbers  hydrometer 

6  S  to  >4  20  2ft  35  46  6S  IOO  ISO  200 


soil  section  for  testing.  After  placement  of  the  sand  in  the  test  box, 
a  1/4"  thick  steel  plate  was  placed  on  the  soil  surface  as  a  uniform  sur¬ 
charge  completely  covering  the  soil  surface. 

The  box  with  the  surcharge  load  was  then  placed  on  the  Load 
platform  in  the  MTS  load  frame.  An  aluminum  block  was  centered  on  the 
steel  plate  and  a  small  pressure  was  applied  by  actuating  the  MTS  hydraulic 
loading  system.  Once  this  small  pressure  was  applied,  the  MTS  Load 
system  was  programmed  to  vibrate  the  sand  under  a  sawtooth  type  load 
with  a  peak  intensity  of  150  lbs.  at  a  frequency  of  7  cps  until  0.  2  of  an 
inch  downward  movement  of  the  soil  surface  occurred  (determined  by 
observing  the  sinkage  trace  on  the  Visicorder).  The  application  of  an 
initial  pressure  insured  that  contact  was  maintained  by  the  load  piston 
throughout  the  vibration  process.  A  number  of  trial  and  error  procedures 
were  used  in  arriving  at  the  peak  intensity  and  frequency  selected  for 
compacting  the  soil.  Each  time  the  sand  was  to  be  vibrated,  the  sand 
was  thoroughly  mixed  to  bring  it  back  to  a  loose  condition. 

b.  Clay 

Buckshot  clay,  a  fine  grained  soil  found  in  the  Vicksburg, 
Mississippi  region,  was  shipped  to  the  University  in  a  premixed  condition 
at  a  moisture  content  of  approximately  277*  to  287o  in  sealed  containers 
by  the  U.S.  Army  Corps  of  Engineers  (WES)  located  in  Vicksburg.  The 
se^ction  of  buckshot  clay  which  exhibits  a  consistent  water  contentT 
strength  relationship  vas  made  to  permit  comparisons  with  previous  and 
future  tests  conducted  by  the  University  on  the  same  soil. 

Before  compaction,  the  test  box  was  lined  with  a  plastic  sheet 
to  reduce  moisture  loss.  An  approximate  2-1/2  inch  layer  of  loose 
soil  was  placed  in  the  box  and  compacted  by  100  blows  from  a  standard 
Proctor  hammer.  A  total  of  4  layers  were  so  placed,  which  constituted 
the  test  section.  Each  layer  was  scarified  before  the  placement  oi  the 
next  layer.  The  soil  surface  was  leveled  and  suitable  measurements  made 
to  define  the  wet  density  of  the  placed  soil. 
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The  top  of  the  test  box  was  then  covered  with  a  plastic  sheet 
and  stored  in  a  moist  room  for  3  days  prior  to  testing.  Moisture  content 
determinations  were  made  on  each  test  specimen  after  the  loading  tests. 
Two  moisture  content  samples  were  taken  from  each  test  specimen. 

Tost  Procedure 

a.  Sand 

Once  the  sand  soil  test  section  had  been  prepared,  either  a 
cone  penetration  test  or  a  plate  test  was  conducted.  For  those  test 
sections  for  which  only  cone  penetration  tests  were  performed,  three 
penetrations  were  made:  one  at  the  center  of  the  section,  and  one  toward 
each  end  of  the  section  approximately  5"  from  the  center.  Both  sinkage 
and  load  were  recorded  as  a  function  of  time  to  a  penetration  of  2-1/2". 

For  the  single  plate  tests,  the  plate  loading  was  applied  at 
the  center  of  the  test  section,  with  the  sinkage  and  load  recorded  as  a 
function  of  time.  The  plate  was  penetrated  at  a  constant  rate  of  12.  5"/sec. 
to  a  depth  of  1-1/4".  After  this  plate  test,  three  penetration  tests  were 
conducted:  one  toward  each  end  of  the  section  approximately  5-1/2"  from 
the  center,  and  one  where  the  single  plate  had  penetrated.  The  purpose 
of  the  penetration  test  in  the  previously  loaded  plate  region  was  to  obtain 
information  on  the  modified  soil  properties. 

In  the  twin  plate  testing,  three  different  spacings  between 
plates  were  used,  as  described  previously.  In  each  twin  plate  test,  both 
plates  were  penetrated  simultaneously,  with  sinkage  and  load  (total  load 
on  both  plates)  recorded  as  a  function  of  time  on  the  Visicorder.  After 
the  twin  plate  test,  one  plate  was  removed  from  the  test  frame  and  the 
remaining  plate  was  used  to  conduct  a  single  plate  test  in  one  of  the 
previous  plate  load  regions.  A  cone  penetration  test  was  then  conducted 
in  the  other  previously  loaded  plate  region.  Two  cone  penetration  tests 
were  also  performed  toward  each  end  of  the  test  section  approximately 
5-1/2"  from  the  center. 
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Calibration  checks  on  load  and  sinkage  were  made  prior  to 
and  after  testing  on  each  test  section.  This  calibration  was 
accomplished  by  substitution  of  a  known  resistance  into  the  circuitry 
of  the  MTS  system. 

b.  Clay 

After  the  test  box  was  placed  on  the  MTS  platform,  cone 
penetration  tests  were  first  run  in  each  of  the  four  corners  of  the  box 
(cone  test  Nos.  1,  2,  3,  and  4),  approximately  2"  from  the  edge.  The 
plate  loading  test  (either  single  or  twin)  was  then  conducted  at  the  center 
of  the  test  section.  After  each  plate  test  (single  or  twin),  a  cone 
penetration  test  (cone  test  No.  5)  was  conducted  in  the  previously-loaded 
plate  region.  For  all  the  twin  plate  tests,  one  piate  was  removed  and  a 
subsequent  plate  test  was  conducted  within  the  remaining  previous  plate 
test  region.  Finally,  an  additional  cone  penetration  test  (cone  test  No.  6) 
was  made  in  a  relatively  undisturbed  region. 

Both  sinkage  and  force  were  recorded  as  a  function  of  time 
on  the  Visicorder.  For  twin  plate  tests,  the  recorded  force  was  the 
total  force  on  both  plates. 

Uniformity  Tests 

a.  Sand 


The  results  of  the  cone  penetration  tests  are  given  in 
Table  26,  with  the  penetration  resistances  given  at  depths  of  1-1/4”, 
1-3/4”,  and  2-1/4”  as  determined  from  the  Visicorder  traces.  The 
average  soil  cone  penetration  resistance  for  each  penetration  tevt  was 
defined  as 


CPR 


N 


,cp\ 


1-1/4" 


+  CPR 


ft  1-3/4" 


+  CPR 


(6  2 


(A  -  1) 
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where  N  =  1,  2,  3  . 


TABLE  26 


Sand  Tests  -  Cone  Penetration  Test  Results 
(Initial  Loading) 


Tost  Tost  Resistance  (CPR)  in 


Setjucnc  v 
No. 

Cimimcnl  s 

Cone 

No, 

1-1/4 

lbs.  at 

"  1-3/4"  2 

-1/4" 

cpRN 

(tbs) 

CPRavg 
(lbs)  B 

Date 

IS 

Cone  T ost 

1 

5.  8 

8.  1 

1  1.  3 

8.4 

Only 

2 

5,  1 

7.6 

10.6 

7.  8 

(8.0) 

4-3-70 

3 

5.4 

7.  4 

10.  3 

7.  6 

2S 

1  st  Single 

1 

_ 

in  plate  area- 

Plato  T ost 

2 

5.  8 

8.  8 

12.  8 

9.  1 

(8.2) 

4-3-70 

3 

4.  8 

6.  9 

9.  9 

7.  2 

3S 

2nd  Single 

1 

„ _ 

in  plate  area- 

Plato  T ost 

2 

6.  5 

10.  2 

15.  4 

10.  7 

(11.  1) 

4-3-70 

3 

7.  8 

11.  2 

15.  5 

11.  5 

4S 

Cone  T osts 

1 

7.  5 

1 1.  4 

17.  5 

12.  1 

Only,  *Box 

2 

9.  1 

13.  3 

19.  2 

13.9 

(13.7) 

4-3-70 

Jiggled 

3 

10.  3 

15.0 

20.  3 

15.  2 

5S 

3rd  Single 

1 

« - 

in  plate  arca- 

Plato  Tost 

2 

6.  3 

9.0 

12.  3 

9.  2 

(9.  7) 

4-3-70 

3 

6.  4 

9.  8 

34.  3 

10.  2 

6S 

Cone  T ost 

1 

6.  3 

8.6 

11.  3 

8.7 

Only 

2 

5.0 

7.  1 

10.0 

7.4 

(8.  2) 

4-3-70 

3 

5.4 

CO 

00 

11.8 

8.  5 

7S 

1st  Twin 

1 

„ _ 

in  plate  area- 

at  1-1/2D 

2 

6.2 

9.  5 

13.  5 

9.  7 

(9.  0) 

4-3-70 

3 

5.  3 

8.4 

11.  3 

8.  3 

8S 

2nd  Twin 

1 

m - 

in  plate  area 

at  1-1/2D 

2 

5.  8 

8.  1 

10.7 

8.  2 

(9.4) 

4-3-70 

3 

7.4 

10.  5 

14.0 

10.6 

9S 

Cone  Tost 

1 

4.7 

6.6 

9.8 

7.0 

Only 

2 

4.6 

6.  3 

8.  7 

6.  5 

(7.6) 

4-3-70 

3 

6.  6 

9.  1 

11.  8 

9.  2 

10S 

1st  Twin 

1 

—  in  plate  area 

at  2-1/2D 

2 

4.  4 

7.  3 

10.  8 

7.  5 

(8.4) 

4-3-70 

3 

6.7 

9.4 

11.7 

9.  3 

1  IS 

2nd  Twin 

1 

- - 

in  plate  area 

at  2-1/2D 

2 

3.6 

5.  1 

7.  1 

5.  3 

(5.4) 

4-3-70 

3 

3.  9 

5.  4 

7.4 

5.6 

1  zs 

1st  Twin 

1 

*  ~  in  plate  area 

— •* 

at  4D 

2 

5.  5 

7.  8 

10.6 

8.0 

(7.0) 

4-3-70 

3 

4.  2 

6.  0 

8.  1 

6.  1 

®  ©  (D 


Typical  Cone  Test  Locations 
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The  average  soil  strength  (or  resistance)  for  each  test  section  in  which 
plate  tests  were  conducted,  the  results  of  which  were  used  in  the  analysis, 
was  given  by 


CPR 


cpr2  +  cpr3 


avg 


(A -2) 


In  general,  the  results  of  the  penetration  tests  in  sand  show 
that  each  test  section  axftibited  a  relatively  uniform  soil  strength,  but 
that  significant  variations  did  occur  in  soil  strength  between  test  sections. 
To  allow  for  these  f trength  variations  between  test  sections,  a  dimension¬ 
less  quantity  (F  /CPR^^  was  used  in  comparing  the  results  of  the  plate 
tests. 


b.  Clay 

The  results  of  the  cone  penetration  tests  are  given  in  Table  27, 
with  the  cone  penetration  resistances  (CPR)  given  at  depths  of  1/2"  and  1", 
in  addition  to  the  average  cone  penetration  resistance  reached  subsequent 
to  the  full  penetration  of  the  cone  tip  (CPRp).  The  average  cone  penetra¬ 
tion  resistance  (CPR^.^.)  for  each  test  section,  the  results  of  which  were 
used  in  the  analysis,  was  determined  by  averaging  the  uniform  penetration 
resistance  reached  in  cone  tests  nos.  1,  2,  3,  and  4. 

Reference  .o  Table  27  shows  the  variation  in  strength  between 
test  sections.  Note  that  two  different  stiffness  ranges  were  selected 
to  insure  chat  the  test  results  would  be  representative  over  a  strength 
range. 

Four  of  the  nine  prepared  clay  test  sections  exhibited  a  slight 
(  c*  i07o)  loss  in  strength  with  depth  (Test  Sequence  Nos.  3C,  5C,  7C, 
and  9C),  and  two  of  the  nine  test  sections  exhibited  a  slight  increase  in 
strength  with  depth  (Test  Sequence  Nos.  1C  and  8C).  Dimensionless 

quantities  (F  /CPR  )  were  used  in  all  comparative  analyses  between 

avg  unif 

tests. 
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TABLE  27 


Clay  Tests  -  Cone  Penetration  Test  Results 
(Initial  Loading) 


Test  Test  CPR  in 

sequence  Cone  lb.  at.  CPRp  CPRur^-£  Date 

No,  Comments  No.  I  /  Z"  *  '  £ 


1C 

lat  Single 

1 

3.  6 

13.  5 

33.  4 

Plate  Test 

3 

4,  6 

16.  9 

37.  5 

3 

5.  1 

17.  5 

35.  5 

4 

4.  1 

14.  5 

33.  6 

5 

« —  in 

plate 

area  m 

6 

... 

33.  9 

2C 

1st  Twin 

1 

3.  0 

10.  7 

22.  2 

at  1-1/ 2D 

2 

3.  0 

8.  7 

15.  7 

3 

— 

---. 

25.0 

4 

3.  0 

10.  2 

22.  3 

5 

m —  in 

plate 

a  rea  — ■» 

6 

21.  3 

3C 

1st  Twin 

1 

... 

— 

.... 

at  2-1/2D 

2 

6.  2 

1 9.  8 

35.  6 

3 

4.  5 

19.  3 

34.  2 

4 

7.  1 

21.  2 

37.4 

9 

« —  in  plate 

area - — 

6 

... 

.... 

34.  4 

4C 

2nd  Single 

1 

3.4 

10.  2 

19.  5 

Plate 

2 

3.  8 

10.  7 

20.  1 

3 

3.  7 

11.0 

21.6 

4 

3.8 

12.  1 

23.  4 

5 

m —  in  plate 

area  — » 

6 

... 

.... 

21.  7 

5C 

1st  Twin 

1 

5.4 

17.  5 

30.6 

at  4D 

2 

1 

3.9 

15.0 

26.  8 

J 

4 

5.0 

15.  2 

27.  1 

5 

4 —  in 

plate 

area  — •» 

6 

... 

.... 

29.0 

6C 

2nd  Twin 

1 

3.  6 

13.9 

29.0 

at  2-1/2D 

2 

4.  5 

14.  8 

28.  8 

3 

3.6 

13.6 

36.  9 

4 

4.  3 

14.  3 

27.0 

5  m  -  in  plate  area  «■ 

6  .  25.0 

7C  3rd  Single  1  5.4  17.3  31.8 

Plate  2  6.  2  18.9  32.3 

3  6.  2  19.  3  30.1 

4  5.4  15.7  30.7 

5  m  in  plate  area  m 

6  .  30.  3 


8C  2nd  Twin  1  3.7  10.5  21.0 

at  1-1/2D  2  6.  2  19.  4  - 


*  Average  of  Cone  Teat  Noa.  1,  2,  3,  and  4 
**  Cone  Teat  No.  6 


35.0  33.4  4-20-70 


21.  3  21.  3  4-20-70 


35.  7  34.  4  4-20-70 


21.  2  21.7  4-22-70 


28.2  29.0  4-22-70 


27.9  25.0  4-22-70 


31.  2  30.3  4-24-70 


20.0  19.3 


26.  1  26. 7 
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APPENDIX  II 

GOVERNING  EQUATIONS  AND  NUMERICAL  PROCEDURES 
FOR  THE  TWO  DIMENSIONAL  PLANE  STRAIN 
MULTIWHEEL  VERTICAL  PULSE  LOAD  PROBLEM 
-  TWIN  WHEEL  SIMULATION 


13  3 


Governing  Equations 


The  equations  of  continuum  elasticity  and  plasticity  used  in  the 
lumped  parameter  iteration  method  are  listed  in  this  section  in  the  form 
applicable  to  the  lumped  parameter  model. 

a.  Dynamic  Equations  of  Motion 


PU(i.j) 


Q^i+1,  j  +  1) 


h/VT 


TnC(itl'i) 


h/ 


VT 


(A -3) 


pV(i,  j) 


Qc(i.j-H)  -  j)  T^Ci+I,  j  +  1)  -  TT|C(i»  j) 


h/VT 


h/VT 


(A -4) 


where 

U  and  V  are  the  displacements  in  the  T|  and  C  directions,  respectively; 
a  and  a,  are  the  noimal  stresses  and  t  ,  is  the  shear  stress; 

n  C  tic 

p  is  the  density  of  the  soil; 

h  is  the  grid  size;  and 

the  dots  indicate  time  derivatives. 


b.  Quadrature  Equations 

t  t-At  -  t-At  (At)3  r  "t-At  ••t., 

U  =  U  +  (At)  U  +  [2U  +  U  1 


t  t-At  •  t-At  (At)3  r  ,■  -t-At  ..t- 

V  =  V  +  (At)  V  +  [2V  +  V  ] 

D 

t-At  .  At  r  lit  -  At  .  •  it  . 


[V^  +  V1] 


(A-5) 
(A -6) 
(A -7) 
(A- 8) 


where 

At  is  the  time  increment,  and  superscript  (t -At)  indicates  the 
variables  of  the  previous  load  increment. 
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c.  Drucker-Prager  Yield  Criterion 

The  criterion  states  that  if  the  yield  function,  f,  as  defined 
below  is  less  than  zero,  the  stress  point  is  elastic,  and  if  f  is  equal  to 
or  greater  than  zero,  the  stress  point  has  yielded. 


Yield  function  =  f  =  +  -  k 


(A-9) 


where 


1=0  +  O.  +  0. 

T1  C  ? 


J=  i  [(0n  •  V8  +  (a'  ■  °*>*  +  (°?  •  °")E  +  6T-r  1 


?  rf 


(A- 10) 


(A- 11) 


+  °C> 


(A- 12) 


2  sin< 


VT  (3  -  sintj>) 


(A- 13) 


k  =  — z_r-C — C(~ =  yield  stress  in  shear 
V3  (3-sin4i) 


(A- 14) 


c  =  cohesion 

<)>  =  friction  angle 

d.  Incremental  Strain-Displacement  Relations 


Ac  AUji.iL: 

n  h/VT 


(A- 15) 


c  h/VT 


(A-  16) 


Ae?(i.  j)  =  0 


(A-  1 7) 


%c(i,j)  = 


2h/VT 


[Ur.lL  + 


2h/VT 


(A- 18) 
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where 


and 


where 


and 


Ae  ,  Ac,,  and  Ae„  are  the  normal  strain  increments, 

T]  C  ? 

Ay  „  is  the  shear  strain  increment, 

■nC 


AU  =  Ul 

- 

(A- 19) 

t 

t-At 

AV  =  V 

-  V 

(A- 20) 

e.  Incremental  Stress-Strain  Relations 


A o  =  \Ae  +  2GAe  -0Q 

T1  T)  r 


zVJ 


+  B 


A =  XAe  +  2GAe^.  -  pQ 


At  =  2GAv  -  pQ 


AW 

2VT 

AW 


+  BAe 


(4-+b 

\2V J 

Maw+bJ 


+  BAe 


,C  nC  K  \2V7/\2V7 


(A-21) 


(A- 22) 


(A- 23) 


Ae  =  Ae  +  Ae, 


n  ltv  I 

^  V6V? 


Q  = 


4G 


!  4  10, 

l-2v 


X  =  Lame's  constant  in  Hooke's  law  = 


Ev 


(l+v)(l-2v) 


G  a  Modulus  of  rigidity  = 


2(1 +v) 

AW  =  Increment  work  done 

=  0  Ae  +  3 -A tr  +  r  Av 
r»  ^  C  C  r'- 


p  = 


1  if  f*0  and  AW>0  (loading) 

0  if  f*0  and  AW<0  (unloading)  or  f<0  (elastic) 


(A- 24) 
(A- 25) 


(A -26) 


(A-27) 
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Then  the  stresses  at  time  t  are: 


t  t-At  t 

o  =  a  +  Aa 

•■nr)  n 

(A-28) 

t  t-At  A  t 

V°c  +  *°e 

(A-29) 

T  ! 

TK  T)C  TlC 

(A- 30) 

f.  Stress  Correction  Equations  for  Perfectly  Plastic  Yielding 

o'  =  (MJO  +  T)  C  (1  +  6o^)  \  -  2ojk] 

•p  C  T1  c  3 

(A-31) 

°'c  =  +  \  +  6“8)  3;  -  2aikl 

(A -3  2) 

r'  =  (1-n  )  t 

nC  c 

(A  -  33) 

where 

J  -  (K-aI)a 

"c  =  2J  +  12af  (k-QlI)s 


(A -34) 


and 

o'  ,  a'  ,  and  t'  are  the  corrected  stresses. 

m  C  vC 

Numerical  Procedure  for  the  Development  of  the  Computer  Program 

For  convenience  in  the  numerical  calculations,  the  above  governing 

equations  are  first  expressed  in  terms  of  dimensionless  variables.  The 

dimensionless  variables  are  formed  in  the  following  manner:  the  variables 

having  a  dimension  of  stress  arc  divided  by  the  yield  stress  in  shear, 

k  (sec  Equation  A- 14);  variables  having  a  dimension  of  length  are  divided 

by  the  width,  a,  of  the  applied  surface  pressure  strip;  and  the  time  is 

divided  by  the  time  duration  of  the  pulse,  t  . 

d 

The  pulse  curve  of  the  applied  pressure  on  the  surface  on  the  soil 
is  approximated  by  linear  segments.  For  each  segment,  the  change  in 
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pressure  through  an  increment  of  time  is  determined.  For  each  incre¬ 
ment  of  time  the  stresses  at  the  fictitious  stress  points  in  the  plane  (j=l) 
are  changed  based  on  the  pressure  increments  of  the  particular  linear 
segment  of  the  approximate  pulse  curve.  Then  the  following  steps  are 
performed  starting  with  the  mass  point  at  (i,  j)  =  (1,  1):  ** 

I  • 

(1)  The  accelerations  U  and  W  at  time  t  are  obtained  by  means 
of  the  dynamic  equations  of  motion,  Equation  (A-3),  using  the  most 
current  stresses  that  are  known  at  time  t  (if  not  known,  those  at  time 
t-At  are  used). 

(2)  The  accelerations  at  time  t  and  the  accelerations,  velocities, 
and  displacements  at  time  t-At  are  then  substituted  into  the  quadrature 
equations.  Equations  (A-5  to  A-8),  to  give  U,  V,  U  and  V  at  time  t.  Then  the 
incremental  displacements,  AU*  and  AV*  are  obtained  from  Equations 

(A- 19)  and  (A- 20). 

(3)  The  stresses  at  the  two  stress  points  immediately  below  the 
mass  point  (i,  j)  are  recalculated  according  to  the  following  steps  for 
each  stress  point: 

(i)  The  yield  indicating  table  is  checked  to  determine  if  the 
stress  point  had  yielded.  (Initially,  the  table  would  indicate  all  stress 
points  to  be  elastic). 

(ii)  The  incremental  strains  at  time  t  are  calculated  using 
the  incremental  strain-displacement  relations.  Equations  (A-15  to  A-18). 

(iii)  The  stress  increments  are  then  calculated  from  the 
incremental  stress  strain  relations,  Equations  (A-21  to  A-23).  The 
stresses  at  the  time  t  arc  then  calculated  by  Equations  (A-28  to  A-30). 

(iv)  The  newly  obtained  stresses  are  then  substituted  into 
the  yield  criterion.  Equation  (A-9).  to  check  if  the  stress  point  had 
yielded.  Th-'  result  is  then  recorded  in  the  yield  indicating  table.  Stress 
correction  is  performed  using  Equations  (A-31  to  A-33),  if  it  was  re¬ 
quired. 

% 

(4)  Steps  1  through  3  are  repeated  for  the  rest  of  the  mass  points, 
proceeding  from  the  axis  of  symmetry  toward  the  right  and  then  row  by 
row  downward. 
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(5)  Using  the  new  stresses  obtained  for  all  the  stress  points, 
steps  1  through  4  are  repeated,  thus  starting  the  iteration  cycle.  This 
is  done  until  the  desired  accuracy  is  reached. 

(6)  Steps  1  through  5  are  repeated  for  the  subsequent  time 
increments,  in  which  the  applied  pressure  on  the  boundary  is  incremented 
according  to  the  vertical  pulse  load  curve. 


APPENDIX  III 

MULTIWHEEL  VERTICAL  PULSE  LOAD  ANALYTICAL 
SINKAGE  PREDICTION  COMPUTER  PROGRAM 
TWIN  WHEEL  SIMULATION 


with 


A.  Sonic  Preliminary  Remarks  A^>out  the  Computer  Program 

B.  Fortran  Source  Listing  of  the  Computer  Program 

C.  List  of  Symbols 
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Some  Preliminary  Remarks  About  the  Computer  Program 


a.  The  load  curve  was  assumed  to  be  symmetrical,  that  is,  the 
peak  occurs  at  one-half  the  time  duration  and  the  shape  of  the  loading 
portion  of  the  curve  is  the  same  as  that  of  the  unloading  portion.  Both 
portions  were  approximated  by  20  equal  linear  segments.  The  'basic 
load  curve"  in  the  program  has  a  peak  load  of  10,  000  psf,  as  shown  in 
Table  28,  and  the  curve  is  scaled  down  proportionately  in  the  program 
for  lower  peak  loads. 


TABLE  28 

BASIC  LOAD  CURVE 


Dimensionless  Time 

0.0 
0.  025 
0.050 
0.075 
0.  100 
0.  125 
0.  150 
0.  175 
0.  200 
0.  225 
0.  250 
0.  275 
0.  300 
0.  325 
0.  350 
0.  375 
0.  400 
0.425 
0.  450 
0.  475 
0.  500 


Applied  Pressure  in  psf 

0.  0 
-20.0 
-110.0 
-270. 0 
-580.  0 
-1050.0 
-1730.0 
-26  30.0 
-3750.  0 
-4890.0 
-58'>0.  0 
-6770. 0 
-7540.0 
-8190. 0 
-8730. 0 
- C)  1 7  0 .  0 
-9500. 0 
-  *740. 0 
-9'*oo.O 
-9  *80.  0 
-  10000. 0 


b.  One  of  the  input  data  items  to  the  program  is  the  time  incre¬ 
ment,  DT.  It  is  calculated  previous  to  u*e  of  the  program  using  the 
stability  criterion  discussed  in  the  Phase  II  Final  Report^.  The  following 
procedure  should  be  followed. 


14! 


-  An  approximate  time  increment  is  obtained  by  the  formula, 

(At)  =  -y— 

approx.  ZCj 

where  h  and  are  the  grid  size  and  the  dilatational  wave  velocity, 
r  es  pectively. 

-  With  the  above  approximate  time  increment  as  a  guide,  a 
smaller  time  increment,  At,  is  chosen  such  that  the  number  of  increments 

‘d 

in  the  entire  load  curve  —  is  a  multiple  of  40.  This  ensures  that  each 

At 

linear  segment  of  the  load  curve  will  have  an  integral  number  of  time 
inc  r  ements . 

c.  Since  the  computer  time  required  to  run  through  a  load  curve 
is  quite  long,  the  computer  program  is  written  such  that  portions  of  the 
load  curve  can  be  run  at  separate  times  by  using  magnetic  tape  input  and 
output.  This  can  be  done  by  specifying  the  starting  load  increment 
number,  LB,  and  the  ending  load  increment  number,  LEN,  and  setting 
up  the  appropriate  tapes. 

d.  The  computer  run  is  monitored  by  printing  out  the  vertical 
normal  stresses,  the  vertical  displacements,  and  the  yield  indicating 
table  of  the  region  under  the  load  intermittently;  the  number  of  load 
increments  skipped  is  given  by  the  Index  ILL  The  vertical  normal 
stresses,  vertical  displacements,  and  yield  indicating  table  of  the  rest 
of  the  region  are  also  printed  out  at  a  less  frequent  rate,  and  the  number 

of  load  increments  skipped  in  this  case  is  given  by  IEI.  The  other  stresses 
and  displacements  are  not  printed  out  because  the  volume  of  print-out 
would  be  prohibitively  large;  however,  at  a  much  less  frequent  rate  all 
results  of  a  load  increment  are  saved  on  the  output  tape.  The  number 
of  load  increments  skipped  for  this  case  is  given  by  JLI.  The  indices 
ILI,  IEI,  and  JLI  are  all  input  data  specified  in  the  last  data  card. 
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e.  The  numerical  value  of  each  element  of  the  yield  indicating 
table  supplies  the  following  information: 

(i)  If  -1.0  <  YIT  <  0.0,  the  stress  point  is  elastic. 

(ii)  If  -10.  015  <  YIT  <  -10.  0,  the  yield  function  is  greater 
than  zero  but  has  not  exceeded  the  tolerance  for  yielding  (which  is  0.  015 
in  this  case),  thus  the  stress  point  is  still  considered  elastic. 

(iii)  If  0.  0  <  YIT  <  10,  000,  the  stress  point  has  yielded  and  is 
loading,  and  no  stress  correction  was  applied.  The  digits  to  the  right 

of  the  decimal  point  give  the  stress  correction  factor  (l-1"  ),  which  is  a 
number  between  0.0  and  1.0.  The  four  digits  left  of  the  decimal  point 
gives  the  value  of  the  yield  function  which  is  a  number  between  0.  0  and  1.  0. 

(iv)  If  30000.0  <  YIT  <  40000.0,  the  yield  function  has  exceeded 
the  tolerances  for  stress  correction,  and  stress  correction  has  been 
applied.  The  digits  other  than  the  ten  thousand  place  digit  gives  the  same 
information  as  (iii). 

(v)  If  20000.0  <  YIT  <  30000.0,  the  yield  function  is  negative 
but  has  not  gone  below  the  tolerance  (-0.015)  for  becoming  elastic  again, 
thus  the  stress  point  is  still  considered  plastic.  The  digits  other  than 
the  ten  thousand  place  digit  give  the  same  information  as  (iii). 

(vi)  If  40000.0  <  YIT  <  70000.0,  the  stress  point  is  plastic  and 
unloading.  The  digits  other  than  the  ten  thousand  place  digit  give  the 
same  information  as  (iii). 

f.  Before  making  a  continuation  run  of  the  computer  program, 
with  the  soil  medium  being  still  all-elastic,  the  value  of  the  cohesion 
may  be  changed  without  affecting  the  results.  However,  since  the 
stresses  are  normalized  with  respect  to  the  yield  stress  in  shear  which 
is  proportional  to  the  cohesion,  the  values  of  the  stresses  must  be  con¬ 
verted  by  the  conversion  factor,  CONV,  during  read-in  of  the  tape  data. 
This  is  done  by  specifying  the  control  index  ICV;  if  conversion  is  desired, 
ICV=1;  if  conversion  is  not  desired,  ICV=-0.  The  value  of  the  cohesion  of 
the  saved  data  on  tape  must  also  be  specified  as  an  input  data. 
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g.  The  procedure  for  running  the  computer  program  is  as  follows; 

(1)  Specify  on  the  first  data  card  a  title  of  length  less  than 
24  characters  including  blank  spaces. 

(2)  Specify  the  next  three  data  cards: 

Second  card  -  specify  five  soil  parameters:  weight  density 
(pcf),  Poisson's  ratio,  Young's  modulus  (psi), 
cohesion  (psf),  and  friction  angle  (degree): 

Third  card  -  specify  six  load  parameters:  time  duration  of 
load  pulse  (seconds),  width  of  loaded  area 
(inches),  peak  pressure  of  load  curve  (vertical 
load/contact  area,  psf),  tolerance  for  un¬ 
loading  (dimensionless),  the  i-index  of  the 
mass  point  at  which  the  loaded  area  starts, 
the  i-index  of  the  mass  point  at  which  the 
loaded  area  ends. 

Fourth  card  -  specify  computational  parameters:  space  grid 
size  (inches),  time  increments  (seconds,  as 
explained  in  Item  b),  number  of  grids  for 
depth  (use  30),  number  of  grids  for  width 
extent. 

(3)  The  next  five  data  cards  specify  the  basic  load  curve  (see 
Sample  Data). 

(4)  Specify  on  the  tenth  data  card  the  value  of  the  cohesion  of 
the  saved  data  on  tape  if  conversion  is  desired.  If  conversion  is  not  de¬ 
sired,  a  blank  card  must  be  supplied  (see  Item  f). 

(5)  Prepare  two  magnetic  tapes;  supply  tape  numbers  for  the 
"setup  cards".  Specify  on  the  second  to  the  last  data  card  these  two  tape 
numbers  arranged  with  the  one  for  Unit  9  first. 

(6)  Specify  the  last  data  card:  Twelve  integer  numbers  are 
required.  Eight  of  them  must  be  the  following,  the  other  four  may  be 
specified  accordingly. 
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1  (LEN)  (ILI)  (I El)  (JLI)  0  0  0  1  9  10  0 


LEN  is  the  ending  load  increment  number.  It  is  desirable 
to  have  it  equal  to  multiples  of  JLI.  It  is  also  desirable  to  have  JLI  being 
multiples  of  IEI  and  IEI  being  multiples  of  ILI.  It  is  suggested  to  use 
LEN=50,  JLI=50,  IEI=10,  and  ILI=5. 

(7)  To  make  a  continuation  run,  the  last  data  card  is  the 
only  card  needed  to  be  changed.  The  information  needed  to  change  this 
card  is  always  printed  at  the  end  of  the  preceding  run.  Only  LEN  needs 
to  be  supplied  by  the  operator. 
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FORTRAN  SOURCE  LISTING  OF  THE  COMPUTER  PROGRAM 
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ISETUP  9 
SSFTUP  10 
«IFi JOH 


(TAPE  JO.) 
(TAPF  -JO.) 
NAP,  ALT  I  U 


tIBFTC  MlJVL  ;  M94  ,v«  J 

C  HENRY  LUW!  JO  UP  PESEAPCM  INST. 

C  PRHGPV'  FCR  "1ULT I-WHEEL  VERTICAL  LOAD  SINKrGE  PREDICTION 
C  PAIN  CALLING  THE  TWJ  OVERLAID  SUOR  JUTI  iES. 

CPV  -i  r./.»UMMY/  VD<38»?8),  Uf>(  38,28),  V<38,2RI,  U(38,28), 

1  :.X{  38, 28  )  ,  SY{  38, 28  I  ,  SXY(  38,2  8)  , 

2  IJT(  38,28  ), UDTI38, 28), JODT(38,2e),  Ul(38,28), 

3  VT(38,2«),VOT(38,28),VODT( 38,28),  VI(3b,28), 

4  SXT<  38,28),SYT(38,?H),SXYT(38,28>,YIT(38,2h) , 

0  $7(38,28),  W{  38,  20) ,PCU(21 ) ,PIN(20) ,NLO(30) 

CO^MU-t/SOM/  WR.KRM,  T",DT,N,M,Nl,ril,  18,  IE  1,  SI  ,  SI  K  ,001  ,HH  ,  OTT ,  AL  , 
i  GG,P'),C03,P,  AP,44P,LPP,SQ*2,lS,N7D,NlN.KT,KU,KN,WaT 

COP-i'J  N/CNUTR/LB.LEN,  IL  I,  ILP,  IE  I,  I E  ,  JL I ,  JLL ,  NTI ,!  J  T2  ,NT  3 ,  NTO ,  NT  I  , 
l  :iI7 ,  NOT 

CALL  S iJ -t T i.  («>99) 

CALL  SCALC 
99  STOP 
Fr.n 


THE  FOLLOWING  IS  A  ‘‘mP  SUR-PRNGRAM  tq  DFFINE  THE  FILES  FOR  THE  TAPE 
UNITS  u,  9,  10. 


$  I  R  w  A  P 

N V L AP 

F  ITRY 

.tJNOR. 

. UNO 8 . 

p;r 

U  Jl'Orf 

UNI  TOti 

FILE 

,  >3  (  1), READY,  imUT,HIN,«3LK»?56 

FN'PV 

.O’  09, 

• UN09 . 

P?F 

ON ! T  09 

UNI r09 

FILF 

,3(2  )  ,RCAOY,  I  J3UT,  aiN,P-LK*2S6 

f  !■  RY 

,U  110. 

•  U'liO. 

p/F 

oo:t ic 

U  M I  T 1 0 

FILF 

,8(3)  .READY,  INOU' ,  B  I N ,  RLK  *  2<?6 

FN~ 

tQRiGIN  SEGl 

MttFTC  PVLDiT  ••94,XK7 

C  THIS  SHARHUT I OF  P  i:  AOS  IN  PARAMETER  NATA,  RFAC'S  DATA  FROM  TAPF,  AND 
C  DOFS  O’HFk  PRFL! ’UNARY  CALCULATIONS. 

SU3nnO’-I  If  <  DATA  (  *) 

C2VVIU  </OUMMY/  VD( 38 »  28 ) ,  UD(  38,  28),  VC  38.28),  U(38,28), 

;  SX(38, 28),  SY( 38,28) ,  $XY<38,28), 

2  UT ( 38,28), UDT (38,28), UDDT (38,28),  UI(38,?8), 

3  V’(38,?8),VDT( 38, 28 ) , VDOT ( 38 , 2 8 ) ,  VI(38,28>, 


o  o  o  o 


4  SXT(38,2rt),SYT(la,?8),SXYT(38,2H),YIT(38,28), 

8  S7(3A,28),  W( 38,28), PCU<21) ,PIN<20) ,NU)(30) 

CPV'1CN/S  0  \T/  W  R , WRM  *  T  M » VT  ,N,M,N1,M1,  Id,  If: ’I,  SI ,  SI K  ,C01 ,  HH  »OTT*AL* 

1  -  GG»  P .’)»  C  J 3,  P ,  AP,  A4P,LPP,  SWR2,LS,N’0,N:  1 » KT » KU » KN » WOT 

CU*M0  ;/CO\|TR/ia,LCN,  ILI,  IIP,  IF  I ,  I F  ,  JL  I ,  JLL  ,  NT1  ,NT2  ,  NT  3 ,  NTO ,  NTI  , 

l  MT.'JJT 

DIMENSION  HTLf(4) 

C  READ  H  ViO  V. R  I  T  F  T  I T L E  .IF  -rHF  RUN. 

RE  A"*  (6,12V)  TITLE 
WRI’F  (6  ,  144  )  TITLE 

READ  IN  UU4  -  FIRST  REAP  CONSISTS  OF  SOIL  P^RmMFTEKS 
SECOND  RFAn  CONSISTS  OF  LOAD  PARAMETERS 
TMIki)  RFAD  CONSISTS  OF  COMPUTATIONAL  PARAMETERS 
FOURTH  RF AO  CONSISTS  OF  LOAD  CURVE 
REAM  (5,100)  RHU, PJ, F  *  C »  aHl 
RF  \  '  (8,  101 )  T0,FPL,PYP,WOT, IR, IP  ) 

RE  A^  (5,102)  H, Dr , M, N 
R  F  A''  (8,140)  (PCU<  I  ),  1*1,21  ) 

N1  *’■!-! 

Ml  *  v._ i 

C  CALCUL  AT  CTHCR  COIL  PARVMTRS  AND  P31  |T  THEM  OUT  FOR  REFERENCE 
G* 144 . »F / (  I.4Pj)/2. 

C2»SCRT(G*3:.2/RHO) 

Cl*:2*(2.*(I.-PJ)/(  l.-2.*P’J)  I **0.5 
AL*?.*PO«G/( l.-2.*P0> 

WRIT  (6,103) 

wi<rr  (6,104)  kho, po, f,g,c2,ci 

WPI*F  (6,105)  C.PHI 
PHI*PHI*3, 1415*27/180. 

CCM3.-SI  I(PHI) )*3.**0.5 
AP*2.*SlN(PHl )/CC 
AAP= :p*ap 

YS*6.*C*C JS(PHI)/CC 
WRIT  (6,106)  AP ,YS 
WR  l  r  c  (6,111) 

WRIT  (6,113)  TO,FPL,PKP, 13, IFM 
WRIT  (6,107) 

WRIT  ( o  ,  108  )  HjO^M,  N.WQT 
WRIT  (6,10V) 

WR  IT  (6 , 1  1U  )  rPL,TO,YS 

C  NCN-GI MCNS 1 0‘lALl  7  InG  ML  PARAMETERS  AND  CMCULATF  SOME  CONSTANTS  THAT 
C  WILL  OF  USE  >  In  THC  LATFK  L^  IPS 
AL* 'L/YS 
G*0/YS 
G  G  *  2  »  •  .» 

H*H/FPL 
$CR?»$CR7 ( ? . ) 

HH»H • SCR  ? 

FPL*FPL/ 12. 

OTT*D’/TO 

COl«YS*Tl)*TD*32.2/(RHO.Ft>L*FPL  ) 

Cf  3«  AP*|  1.4PD/I  I.-2.*P0» 

P«GS/(0.543.*C1)J*AP) 
dtg*ot"»  ;g 
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DTP  =  f;T~*P 

C  MODIFY  BASIC  LOAD  CURVE  FOR  th{:  SPECIFIED  PEAK  LOAD,  CALCULATE  LOAD 
C  INCREMFTG  vn  * ;» I *1 6*4 S  ION  AL  IZ  IWG  DOTH,  CALCULATE  the  number  OF 
C  LOAD  INC^F  ir *JTS  I*i  A  LINE  SFGMFMT.  PRI  NT  OUT  ALL  FOR  REFERENCE 
WPIrr  (6, lib) 

ND»O.02r»/DTT*o.OOI 

fm:*nm 

*jjn  =  Ao«  j i'i 

tt.)=d. 

;>iq»o 

FPsPKP/IUOOO. 

PC*PCU(1 ) 

on  ift?  j*i,?o 

TT0  =  TTTH0.02e. 

nc=nc+  ii  ki 

JJ=.H1 

PI=(°CU( JJ)-PC)*FP/F  1IN 
PC*PCU( J  J) 

PCm=PC*FP 

p i m ( j)-p:/ys 

PCU( JJ)=PCM/Y< 

163  WHI'F  (6,116)  J , TTD, PC, PCM, P I , PCU( JJ ) ,  P  I  N( J) ,  NO 
C  READ  IN  VALUT  OF  THE  COHr:SlDM  OF  THE  TAPE  DaTA  WHICH  IS  USED  FUR 

C  CONVERSION. 

R  F  AD  ( 5 ,  1 A  0  )  r  T 

C  RF AD  IN  TA^E  JUMPERS  OF  THE  TAPE  SETUP  VI  UNIT  9  AND  1 0 , RE SPECT I VEL Y 
READ  (5,  1 A  1 )  MIT, MOT 

C  READ  IN  STARTING  LOAD  INCREMENT  JUMPER  AND  FiOING  LOAD  INCREMENT 
C  NUMRCR,  AND  OTHER  CONTROL  INDICES.  CALCULATE  ThE  PR!  JT  CONTROL 
C  INDICES,  THE  Tf;.jr  OF  7HF  PULSF,  AND  THG  INITIAL  APPLIE)  PRESSURE  SI, 
C  AND  PR  I*,T  OUT  FOR  RFFEREMCE « 

READ  (6,141)  LB.LEN, ILI, Ifcl, JLI*NTI,NT2,NT3,LPP,  1TI , JTO.ICV 

ILP  =  LPMLI-1 

IF=LB« IEI-1 

JLL=LE+JLI-1 

HL=LH-1 

T  P-RL  * DT 

Te  =  TM/T:)/0. 025*1.0 
KT='rn*O.OOl 

IF  (Kt.lT. 41)  SO  TO  170 
S I K  *0 . 

S  I  =T>  . 

GL  TQ  17  1 
170  TK  =KT 

KM  =  KT 

L  S  =  K  T  *  N I  N 

KU=1 

IF  (KT.LE.20)  GO  TO  166 
KT=42-KT 
KU  =  -l 
Kf,  =  KT— I 

165  SIK=FLOAT(KU)*PIM(KN> 

ktu=kt*kij 

Sl  =  (T J-TK)*(PCU(KTU)-PCU(KT) ) ♦PCU( K  T ) ♦ S IK 
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o  o  n  o 


171  wei’r  (6*11/) 

K  P  I  ”  F  L3.LEN,SI,SIK,TM 

WR  rf  (f.  ,UM  LK.lEN, ILT, !EI, JLI.NT1,  IT2,NT3,LPP,  IT  I  ,NT3*ICV 
C  !  F  CnwFRSI  1  i  Ii»  NEEDED.  CALCULATE  the  CONVERSION  FACTJR  FOP  THE 
C  STRESS^. 

IF  ( 1C7.E0.0)  GJ  ~3  159 
CC(/=CT/C 
169  CCNTIUC 

IF  (L!i.lr.l)  GO  T)  164 

C  IF  THIS  IS  THE  V  Fr'Y  FIRST  KIJN,  ALL  STRESSES  AND.  DI  SPLACF.MENTS  ARE 
C  FIRST  SFT  FUUAL  TO  /  F  R  L) ,  AND  THE  YIELD  INDICATING  MATRIX  IS  SET  TO  BE 
C  ELASTIC 

DO  7  J*1*N 
nc  7  1  =  1  f  -1 
ll(I»J)-0. 
v(  :,.i )  =o. 

»D I  I,J)*0. 

V0(  I,J)-0. 

in n=o. 

VT(  I,J  )=0. 

UD’I  !»J  »«0, 

V0~(I,J>=0. 

uc;n c :  * j ) =0. 
vmT  ( i  *  j  )  *o. 
hi  < i . j i =o . 
v!  ( :*j )=o. 

S  X  (|  ,  J  )  =  0  • 

SY(  I  »  J  >  =  0  * 

SXY(T»J)=0. 

SX>(  I ,  J  )  =0 . 

SYT( Tf J)3C. 

S  XYT  (  I  *  J  )  *0  . 

SMI,J)=0. 

7  YIU  !  ,  J  I  «  — 1  • 

RETURN 

C  RE  AkR  AN  SF  DESIOmJN  JF  UNIT  OF  TAPFS  DEPENDING  ON  NT  I  AND  NTO 
164  IF  CJTI.ro. 7)  GO  TO  173 
I LL  =  N  IT 
NIT-NCT 
NO~=ILL 

173  Wkl-F  (6, 136)  SIT 

IF  THIS  IS  *  CO.!*  rilJATITl  RUM,  THE  STRESSES*  D I SPLACE ME 'IT? ,  AND 

YIELD  I  *niC  .TT'JG  TAPLF  nF  THE  PRECEDING  RUN  ARE  READ  IN  FROM  INPU 
’APF.  j  \  T  PFAU  II,  IF  M3  REDUNDANCY  OCCURRED,  ARE  SAVED  IN  DISK 
UNI*  I  F'N  TRANSFER  TO  0'JTPUT  TAPF. 

DC  162  1=1,  *ri 

RFAn  (  ;T  M  UTtu;)T,UO!)T,UI  ,VT,  VDT,VOI.»T,  VI,  SXT,SYT,SXYT,YIT,SZ,W 
I LL=S YT( I . 1 ) 

WRI’C  (6,139)  ILL 
IF  (  I  .  IE.  IT2)  GO  TO  162 

WPITF  (d  )  l  JT  ,urjT,UOOT,UI,VT,  VDT.VDDT,  VI  ,  SXT  ,  SYT,  SXYT  ,YIT,SZ,W 
162  CONTINUE 

C  READ  FROM  DISK  IlMIT  8  THE  DFSIRE  STARTING  STRESSES  AND  DISPLACEMENTS 
C  FOR  THE  CO  IT  riUATJON  RUN. 
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ir  (NT2.E0.0)  GO  n  1 75 
REWIND  8 

174  RE4P  (81  UT,UOT,UOOT,UI,VT,VOT,VJOT,VI ,  SXT , SY7 , SX YT , Yl T , S 7 , W 

1 LL’SYT ( 1 » 1 ) 

REWIND  NT! 

175  IF  ULL.'lE.lLR-i))  GO  99 
NTI  =  0 

C  ECUATE  rHE  CURRENT  DISPLACEMENTS  AMO  STRESSES  WITH  THE  PREVIOUSLY 
C  SAVED  DISPLACEMENTS  AND  STRESSES.  MAKE  STRESS  CO'JVFRSIDM  IF  THF  TAPE 
C  OAT  A  IS  NORMAL  I?  FO  DIFFERENTLY  • 

DC  8  J*l ,N 
DO  H  1-1,1 
UII.JI*  UT  t I » J  ) 
v ( i »  j  >*v  r ( :,j) 

UP  1 1  #  j>*ur>Tt  I  *  J  I 
VDITt JI»VDT( It J) 

IF  ( ICV.EO.O)  on  T0  176 
$XY7( I, J)*SXYT( I,J)«C3NV 
S Y T ( 1 1  J ) x  SYT( I,J)*CnNV 
SXT(t,J)=  SXT  t  I  tsJ ) *CONV 

176  SX(!,J»=  S  XT (  I ,  J  ) 

SY( !,J»x  SYT(  I, J) 

s<y( :,j)=sxyt(  i, ji 
8  CLNMNlJF 
RETURN 
99  RETURN  . 

100  F0RMAT(fR»l,r6.2»3F8.1) 

101  FC'RNfi-  (?F1Q.3,F10.2,E10.2»2I5I 

102  FGR-U-  (  r8.3,F12.6,2I5) 

103  FORMA-MlMli  19X,  15HSDIL  PROPFkTICS) 

104  Ff)R^A-(23Xt7MnENSITY,  17Xt5HPH0  =,F10.1,10H  LBS/CU-FT/ 

1  2  3X , 14HP0I SSOMS  PATIO, 11X.4HPJ  «,F10.2/ 

2  2  3X,  !4HYfiU'lGS  MOr'JL JS,  12X ,  3HE  *,FiO. 1,4H  PSI/ 

3  2  3X  ,  13KSHr-AR  MODULUS  *  1 IX,  3HG  »,F10.1,10H  L8S/SJ-FT/ 

4  23X,?9HSHEAR  WAVE  VELOCITY  C2  «,F10.1,7H  FT/SEC/ 

5  2  JX  »  29HDI L  AT  AT IQN4L  WAVE  VEL.  Cl  «,F10.1,7H  FT/SEC//) 

105  FCR',A~(23X,8HC0HE$Ir>N,18X,3HC  «,F10.1,ICH  LBS/SU-F  T/23X,  14HFRICU0 
IN  ANGLE, 10X,5HPMI  x,F10.l,4H  OES//J 

106  EORv  AT  (2  3X,29HFilR  YIELD  CRITERIA  ALPHA  *  *F  16*8/49 X ,3 HK  «,El6.b, 
1 1 OH  LRS/SC-rT//) 

107  F  C'RU  \  ~  { l  HO  ,  19  X  ,  24HC0MPU’’  A  r  I  ON  AL  PARAMFTERS) 

108  FOR»’-’-(2  3X,  10HSPACE  MF$H* 

116X.3FH  * , F  1C .4 , 3H  IN/23X , 29H0ASI C  t|ve  INCREMENT  OT  », 

2F10.7.4H  SEC//23X, 11HUUM9FR  JF  I,15X,3MM  » , 1 4/2 3X , 1 1HNUMPER  OF  J, 
315X.3HN  »,  I<./2  J<.29HUNLDAniNG  TOLERANCE  WOT  «,lPfcl0.2//l 

109  Ft  AMA'UHO.IDX.SOHCHARAC'FRISTIC  Pj.»AMFtERS  F'V  NO  4-  )  I  ME  USI  COAL  1  l  i 
l  MG  I 

110  FfiRM4^lt3X,  14HLENGTM  »  FPL  »,F10.2,3M  IN/2  3* » 1 4HT I  M*"  *  TO  «,F10 

1.2, AN  iFC/23X,  l<«HSTRf  «  «  »,F10.2,10H  US/SD-F 7 / // ) 

111  F(R'*A’  t  lHO,  19X,  I5HL JAO  PARAMETERS) 

112  FfN'UT  ( 1M/// ) 

113  FORM!  <?3X,29h'IME  DURATION  OF  PULSE  TD  ».ri9.).«H  SEC/ 

1  ?»X,2'RH'IVE  F  nuTFR  I  NT  LENGTH  FPL  »,HC.3,3H  lu/ 

2  23X,29Mi»r.’.K  PnFSSURE  OF  PULSf  P*p  «,F10.1,10m  L^S/SC-FT/ 
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3  23X,29HL0-\0  BORDER  INDEX-BFGIN  IB  *,110/ 

4  23X,?9HLJA0  BORDER  INOEX-EMD  IEN  *  *  1 1 0/ / » 

115  FORMAT  (6X.99H  AT  DIMENSIONLESS  BASIC  LOAD  DIMENSIONED 

IDIMPNSI.riFO  OIMFNS IONLESS  0 1  MENS  IONLESS/ 105M  J  TIMF 

2  T /▼:)  CURVE  LOAD  CURVE  LOAD  INCREMENT  LOAD 

3  CURVE  LH/'D  INCREMENT//) 

116  FORMAT  (lX,I3tF14.3,E18.1,Fl6.3,Fl6.5,lPC19.7,E18.7,I8) 

117  FORMAT  (  tHl,  lrtX,  3AHPARAMETERS  FOR  THIS  PARTICULAR  RUN )  ..  7 

llfl  FORMAT  (  ?3X,37HSTARTI‘iG  LOAD  INCREMENT  NUMBER  LB  =  ,  I6/23X,37HfcND 

Uf.S  LOAD  INCREMENT  NUMBER  LF.N  =  ,  1 6/2  3X ,  37HSTART I NG  SURFACE)1  PRES 
2 S U R c  SI  =, 1PE17.7. 16H  {  DI MFNSIONLF  SS  ) /23X ,  37MPRESSURF.-  (OR  L 

3OA0)  INCRFMFNT  SlK  * ,  1 PF 1 7 . 7 , 1 6H  (  DIMENSiONLE  SS  ) /23X  .UHSTARTI  N 
AG  T i MF ,20X ,  4M"M  *,OPF12.7,4H  SET / /  ) 

129  rORM-.T  (4A6J 

135  FORMA"  (  1HO,  1HX,  B2HLAST  DATA  CARD  QP  THIS  RUN  IS—// 

22‘>X,b7HlB  IFN  ILI  ICI  JL I  Nri  NT  2  NT  3  LPP  NT  I  NTO  ICV// 
322X.12IS) 

138  FORMA"  (  1H0,  18X,  UHlAPf  NUMBER, 1 6 , 54H  CONTAINS  THfc  RESULTS  OF  LOA 
ID  INCREMENT  NUMBER  ILL  *//> 

139  F  °  R  y  "  (  23X,  IS) 

140  E uRVaT  C  SE 1 1 . ?  ) 

141  FCRM\-  ( IX, 1415) 

144  FOR  < A “  ( IH1////1H0.47X, 44HMULTI-HHCEL  VERTICAL  LOAD  SlNXAGE  PREOIC 
1TI0N//S9X.4A6) 

END 


tCRIGlN  SEG1 

$1 BFTC  MVLCAL  M94.XK7 

C  THIS  SUBROUTINE  DJFS  THE  MAIN  CALCULATIONS  OF  THE  ITERATIONS.  PRINTS 
C  OUT  N«Tnro  RPSULTS.  AMP  SAVES  RESULTS  ON  TAPES. 

SUBROUTINE  SCALC 

COyMQN/ 3UMMY/  V‘)(  38.  28).  UO(38,2B),  V'36,28),  U(38,’b), 

1  SX<38,?8),  S  Y  ( 38.28).  >XY(3R,28). 

2  U"(  38,  ?8).U’0T(  3tt.28),Ur>0T(38,2a).  UI(J«,28), 

3  VT( 38,28), VDT( 38,28) ,V0JT(38, 28) ,  VK36.28), 

4  SXT(38, 28), SYT<  38,  28 ) , SXYT ( 38 ,28 ) , YI T ( 38 ,28 >  , 

5  S7  (  38,  2b  ) ,  W(38,?8),PCU(2l),PIN(2O),NLO(.30) 

COMMON/ 1, nr/  AR,WRMtT«tDT,N,.M,Nl,Ml,in,  IEN, SI  .SlK.COl  ,HH,OTT  ,  AL  , 

l  Gu.PD.C )3,P, AP,AAP,LPP.S0R2,LS,NTD,NI 1,KT,KU,KN,W0T 

CUMv-IN/CJ|TO/Ld,LfN,  ILI,  ILP,  I E I ,  IF ,  .IL  I  ,  JLL  ,NT1  ,NT2  ,NT3  ,NT0,NTI  , 
l  MIT.NU" 

C  ThF  FOLNO^nO  U.*rA  STATFMEN?S  SUPPLY  THE  VARIABLE  FORMATS. 

L  CGI C  *  L  TRN*,LASt 

DI  Me*,S  /  IN  FMS(  3),FMw(3),FMV(  3) 

PAT*.  F'*SU,/1BH(1X,»2,12  F10.6)/, 

i  r»A  m/ipnux,  12,  i?  fic.m/, 

:  F‘*ri;)/>m»*t!2,;2  F  tC .  6 )  / , 

3  FXB  .FX*.,  r  X  ) ,  F  X2  /  6HF  10.5),  1>HF  10. 4), 6  HF  10.3),  6HF  10.21/ 

4  . F«6/4Hf 10.6 )/ 

SJf (Vl.VT.V  ),V4)»S  5RT|  (  (  V  1-V 2 1  ••  2*  (  V2- V  3  ) ••  /♦  I  V 3*  V l )  •  •  2 ) /6.  ♦  V4 * v«. ) 
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STARTIN'1,  POINT  nr  MOST  CVJTFR  LOOP,  FOR  CALCULATION  OF  EACH  LOAD 
INCREMENTS. 

UPL*-1 . 

169  DO  250  ILLSLH,L£N 
LPT=0 
TM=T*«CT 

START  I  NO  POIl"  Jf  THF  1T£°ATI0N  LOOP* 

JA*1 
I  T*0 

6  IT*»T*1 

LAST -  M.F0.2 
DF(J*0. 

START lS*,° POINT  FOR  THE"  LOOP  INCREMENTING  EACH  RUw  G3MC  DOWNWARD 

on  yo  j»i,Ni 
jp=j+ l 

or:  «o  i*i, mi 

CALCULA-r!  AT  E AC-  MASS  POIM",  tHF  ACCELERATIONS  FROM  THE  U W&MIC 
ECU  AT  I  ON  5  CF  M0TI04.  THF  PARTICULAR  EQUATION  TO  USE  DEP.l.'S  0 
LOCATIN'!  OF  THE  MASS  PUINT. 

IF  (J.GT.ll  GO  TO  9 

IT  (I.r.C.m.ANO.I.LF.IEN)  GO  TO  4 

SIX=0. 

SIY=0. 

on  r> 

A  IF  I ( I.MC.IE.OR.IB.EO.ll.ANO.I.NF.IFNl  GO  TO  11 
SIX*«I/2. 

S!Y=SI/2. 

GO  "0  5 
ll  S I x» 5  I 
SI  Y-sS! 

5  UCD*  2 «  *C  01*( -  X ( I P » JP ) *SVY ( I »  JP ) ”S I X ) /HH 
VD0=2.*C:U*(SY( I, JP1*SXY( IP, JP)-STY)/HH 
GC  rr  10 

9  Ut;Q*CMl.(SXUP,JP)-SXCl,J)*SXY(I,JPl-SXY(IP,jn/HH 
VDD*C01*ISY{I,JP)-SY(IP,J)*SXY(IP,JPl-SXY(I,jn/HH 

CALCULATE  THF.  DISPLACEMENTS  AND  VELOCITIES  AT  TiMF  i  FkO*  ,Hfc 
^  ^"sijt'/iIj  HOTT*ur)T«  I ,  J  )  ♦OTT*QTT*<  UDDTII ,  J)  *2.  ♦UrM>) /f». 

VB  =VT(  l,JMJTr.vDT(  T,J1-K)TT*0TT*(  VDDTU  ,J)*?.*V0j)/6. 

ULR  s unT (  !.J»*DT-MUp*ot(  I,  J  I *11001/2. 

VI  R  *VOT(  !,J)«DTT*(vnuTC  I ,  J  )  ♦VDi) )  /2» 
ir  (T-.iiF.l)  GO  TO  12 
UII.JI*  UR 
V  C  T  ,  0 )  *  VB 

ucu ,  j>  =  m0 

VI  (T  ,  ! )  =  V  08 

;  AVERAGE  WITH  THF  DISPLACEMENTS  from  the  preceding  iteration 
12  UII,J)*(UK«UI I,J) >/2. 

V  ( ! ,  J  l  -(  VB*V  (  »,JI  )/?.. 

UC(I ,  M=  (UDr,<UD{  I,  J  )  )  /2. 
v i; {  t  ,  ji=  (  vdmvoi  i ,  j  n  /?. 


THE  PRECFDING  ITERATION 
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C  CALCi-I \r''  T'»f  >• ! SPL4CEr,ENT  INCREMENTS 
ia  ui  ( : ,  ji*  / { 1 1  j i -II" 1 1,  j  i 
)-j  rc  i ,  j ) 

:r  (  ji.  ir. 2)  r.  ’  “"i  17 

urn*  c  i  f  j  i  =ur.-.' 

V**N'l  I  ,J  l*V.-i'’ 

17  lr  ( i " . l r . 2 »  'j  ri  to 

C  CALCULI"  -nr  pror.r-r  CONVERGENCE  OF  TMF  VERTICAL  DISPLACEMENT 
C  HFTwf.l  l  7“ I S  VJU  'HL  PRECEDING  ITERA'IQN.  LOCATE  THE  LARGEST  PERCENT 
C  AND  S-Wr  If  1  1R  L  ' T 1  REFERENCE. 

(-•ii=  *.^r  ( i ,  j  i-i.  > 

:,ji-i.i 

IF  n'.LF.lHi)  i,l  TO  111 
D( u- ru 
!  UJ  =  I  LL 

I7IMT7 
Ilia  r 
J(ja  J 

is  ir  (ov.Lf.crvi  so  rn  ?o 
of  v*c •/ 
l lv* I LL 
l Tv* 
iv*: 

J  V  =  J 

CAL CilL  •’’r  7hF  STF  FSSrS  OF  THE  TWO  STRESS  POINTS  BFLtJW  THF  MASS  POINT. 
INOEx  fC  ijr.ICAT^S  WHICH  STRESS  POINT  IS  BEING  CONSIDERED. 

IC*1  IS  THr  inc  ON  tmf  RIGHT 
I  C  =  7  IS  THE  ONE  ON  tHF  LEFT 
iO  If*l 
K*!  J 
L  -  IP 

2?  kj.  =k-i 
Lr'*L-! 

Fk  TC.  Ml.  1  .  ANO.K.NE.M 
C  SAVC  TH^  STiFSSFS  Jf  THF.  PRECEDING  ITERATION. 

SX".=  S  <  (  K  *  L  1 
5:  VS*  SY(K,l.l 
S/Yrs' XYK.Ll 
YITCaY:r(K,Ll 

c  calculat  THr  s-pain  increments. 
r  x  *  ( 1 1 :  (  k  ,  l  » I  (  k lm  ))  /hh 
f  y*  (  ;:ck',,li-vi(k,lm»i/hh 
I  XY.CIK  v  l,L>-u:(K,LM)*VI(K,L»-Vl(KMtLMn/HH/2. 
r  F  a  A  L  •  (  E  X*l7Y} 

C  CHFCF  'wr  Ylf-Ll)  roiCATir.r.  ’■APLC  TCI  DETERMINE  WHICH  STRES-STRAIN 
C  *<FL-T  !  .’■<  7  j  usr. 

IF  (  VI'S  .GI-.C  ..  A  .O.YI-S.LE.ACCOO.  )  GO  TO  35 
C  STRESS  p-j:  r  IS  ELASTIC. 

<■XlK.Ua  SXT  (K  ,L  MGG*tX4EE 
SYlK.Lla  SY’IK.L  MGS*r.Y4FE 
SXY1K.L1  =r.  XW  (K,l)*GG*rXY 
r.n  r>C 

C  15  si  ^CS  O'U-IT  IS  PLASTIC. 

15  l!  ••  L  =  —  1  • 
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L'o=sxs*r:x4  sys  *ey4Sxys*cxy 
SZS  =  Pf>(  SXS  +  fYS  ) 

SS*rXf,4SYS  +  5?S 

SJ*r.  |E<<  XS,SYS  »j/SfTXYS) 

tU  *CTl-S</SJ/f>. 

S  J2=S.I*2  . 

wr  =  (  ■*  ;/s  J2  +  nr  * r- 1  / al  i  *p 
$X(K,U=  SXT)K,L )  ♦  GOEX4EE-(  SXS/SJ.?4BE)*WE 
SYCK  ,L)  =  SY~  (K  ,L  )+GG*fcY4EE-<  SY  S  /  5J  2*0E  )  *  WE 
SxY(K.L) =SXYT (K,L )+GG*tXY-SXYS*WF/SJ2 

c  so  check  :r-  'hf  stress  point  has  yielded.  this  is  done  only  for  the 

C  FTHL  iTESa'liM.  JA  INDEX  CONTROLS  THE  ENTkY. 

SO  IF  t  .»  0T  »L  AST  )  GO  TO  6'. 

C  CALCUL  \TC  "HE  YITLO  FU'lCT  TO'I 
Sr*SX(K,L>  +  SY(K,U 
SZ/=PC*S  $ 

SS=C$4SZZ 

SJ»rjF(SX(h,L  ‘.SY(K,L ).SZZ,SXY(K,L)  I 
FC=  '.P*SS  *SJ-1 . 

C  CHECK  IT  "itt  v  t  EL  ^  FUNCTION  IS  GREATER  THAN  THE  TOLERANCE  ABOVE  ZERO, 

C  AND  MAKr  THi;  APPROPRIATE  CHANGE  IN  THE  YIELD  INDICATING  TABLE 
IF  (EC. G“. 0.0  IS)  GO  TO  SS 
IT  ( VI'S  »  L  T  •  C  •  )  GO  TO  SI 
IF  ( rr.GT.-0.nis)  GO  to  5* 

S3  IF  (  Fu  ;•  )  G  l  7  1  66 
y:=»-c 

IF  (FC.G’.O.)  YI*-10.-FC 
Gt.  T :  m 

SS  ST«r$S  >5.1  J  NT  HAS  YIELDED.  CHANGE  CONTROL  INDEX  TO  SAVE  DATA  OF  THIS 
LOAD  INCREMENT  ON  TAPE.  CHANGE  THE  FORMAT  OF  'HE  YIELD  INDICATING 
TMJLF  PR  I  NT  OUT. 

55  ir  (  r  ^  47  >  GO  T  O  58 
I F  ( L°P.NE.  I  )  GO  to  GS 
LP3  =  ? 

FFY( 3)=FXS 
C  CHECK  FDR  U  ) LOAD TNG 

56  V.C  =  SX(K,L  )  *rX4SY(K,L  )*EY4SXY{K,L)*F^Y 
IF  (HC.LT.-WQY )  GO  '0  S7 

IF  (.VT.  iT. WOT. jr, YITS.lt. AOCOC.  >  GO  TO  58 

57  UPL=i. 

C  C  ALGOL  A'F  THE  STRESS  CORRECTIO.N  FACTOR. 

58  SSP= tl ,-AP*SS )«*? 

S JS-S  l*S J 

PF=(SJS-S;P)/(2.*SJS4l2.*AAP*SSP) 

PZK-PO*(  (  1.4f,,*AAP)»SS/3.-2.*AP) 

R  IT* l ,-PH 

FCJ=AI NT( IOOOO.*FC) 

IF  (FC.LT.O.)  FCJ=20000.-FCJ 

YI  «rc.i 

IF  (UPL.f.T.O.)  Y!*Y  I  44,0000. 

C  S'RFSi  CORA EC'IOf.  IS  MADE  BY  THE  FOLLOWING  STATEMENTS  ONLY  IF  THE 
C  YIELD  FONCTION  Ic  GREATER  THAN  A  TOLERANCE. 

:r  ( r c , l T • o ♦ 0 ^ 0 )  go  to  6a 
yi  =vy  «  joooo. 
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F  NY { 3 ) =  FX3 

SX(K,L)=  SX(K,L )*«4T4P*K 
SY(K,L>=  SY(K,L)*RAT*PAK 
SXYfK.L) =SXY(K ,L )*R4T 
64  IF  (FRNT)  GO  TO  66 
IF  (L^ST)  Y  !  T  {  K  »  L  ) =Y I 
IF  <  IC.EO.2)  GO  TO  67 

66  K=I 
IC  =  ’ 

Go  rn  n 

67  IF  (K.NF.2)  CO  TO  RO 

C  SET  SYMMP’RY  CONDITIONS  AT  AXIS  OF  SYMMETRY. 

S  Y ( 1 . L )*  SX ( 2 » L ) 

SX(i,L)=$Y<2,L> 

S X Y (  I  , L) =SXY( 2, L  ) 

IF  (Ll’ST)  YIT(l,L)-YI 
80  CONTINUE 

IF  ( IT.lf.2)  GO  TO  6 

WRITE  (6,143)  OFU,  I LU , I TU, I U, JU , DF V , I L V , I T v , ! V , J V 
IF  (LAST)  GO  TO  35 

C  CHECK  IF  THF  CONVERGENCE  IS  GOOD  ENOUGH,  RETURN  TO  CALCULATE  ANOTHER 
C  ITERATION  IF  NOT  ACCURATE  ENOUGH. 

IF  ( (DFJ.GT. 0.002. OR. DFV.GT. 0.002). 4ND.IT.LT. 7)  GO  TO  6 
C  IF  ACCURATE  FNJUGH,  AOJUSt  JA  INDEX  AND  CALCULATE  FINAL  ITERATION. 
JA  =  2 
GfJ  T0  6 

C  THE  FOLLOWING  LOOP  SAVES  ALL  THE  DISPLACEMENTS  AND  STRESSES  FOR  THE 
C  CALCULA'ION  OF  THE  (EXT  TIME  INCREMENT. 

C  WITHIN  THIS  LOOP,  THE  VERTICAL  DISPLACEMENTS  AND  STRESSES  ARE  ALSO 
C  CALCULATED  FROM  THE  DIAGONAL  01 $PLAEMENTS  AND  STRESSES. 

85  DO  70  J*I,  I 
DO  RO  1=  l,M 

UT ( I , J  )  »  U (  I , J  ) 

V’-(I,J)=  V(I,J) 

UDT ( I , J ) *  UD( I , J ) 

VnT ( ! , J  )  *  VD ( I , J ) 

W( I , J ) s (U ( I » J  )  +  V  ( I,J) )* 1 .0E06/S0R2 
IT  (J.GT.l)  GO  TO  89 

IF  ( I ,GT . IEM.OR. ( I.LC. I B . AND . 1 8 .NE . I ) )  GO  TO  87 
SZ ( T ,1)»SI*100. 

GO  ’(j  90 
87  s z  <  r , l >  =  0 . 

GL  TU  90 

89  SXT(i,j)=  SX  C  I , J  » 

SYMI,J)  =  SY(I,J) 

SXYT ( I , J ) *  SXY (  I , J) 

SM  I  ,J  )  =  (  (SX(  I,  J)+SY(  I,  J)  )/2.+SXY(  I,  J  H*1C0. 

90  CUNT.iJF 

C  IF  THE  ENG  OF  A  l INEAR  SEGMFNT  OF  THE  LOAD  CURVE  IS  REACHED,  SIK  IS 
C  CHANGED  T1  THE  APPLIED  RRFSSURE  INCREMENT  OF  THE  NEXT  LlNFAR  SEGMENT. 
IF  (  ILL.LT.LS  )  GO  TO  182 
IF  ( ILL.GE.NTD)  GO  TO  179 
L  S  *  L  S  ♦  N I  N 
KT=KT.KJ 
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IF  (KT  *  N  E  «  2 1 )  GO  TO  178 

KN»?  1 

KU=-1 

178  KN*KN4KLI 
S1K*FL0A'(KU)*P1N(KN» 

SI*PCU(KTJ 

GO  T0  1X2 

C  179  IF  F'-O  (-)F  THF  LUAO  CURVE  OCCURS,  THE  SURFACE  PRESSURE  IS  NO 
C  LUNGER  INCREMENTED. 

179  ST  =0. 

S 1 K“0  . 

182  CONTINUE 

IF  (LRP.LT.2)  GO  TO  185 

C  THE  STRESSES  AN1)  Cl  SPLACEMENTS  OF  THE  LOAD  I NCREMFNT  1 :4  WHICH  THE 
C  FIRST  PLASTIC  POINT  OCCURS  ARE  SAVED  ON  THE  UUTPUT  TAPE. 

L  P°-0 

wEI’-r  (6,125)  ILL, IT 
WRI’F  (6,119) 

LPT=1 
GO  'rG  245 
185  CONTINUE 

C  PRINT  OUT  rn r  VFR’ICU  STRESSES  AND  OISPLACEMENTS  IN  THE  REGION  UNDER 
C  THE  LfiAnpr  AFFA  FUR  MONITORING  THE  COMPUTER  RUN. 

193  IF  (ILL.LT.RP)  GO  ’U  240 
C  SELFCTI'N  OF  F  OP  NAT  FOR  THE  PRINT  OUT. 

IF  ( AHSISZ( !FM,1) l.LT.10.  }  GO  TO  209 
IF  ( ARM S?( IFN, 1) J.LT.100.)  GO  TO  2 UB 
Ff’S(  3)=FX4 
GO  TO  212 

208  F F  S  (  3  )  =F  X6 
GC  Tfj  212 

209  F  F  S  (  3  I  ®  F  X  6 

212  IT  (ARSlWdB4ltin.LT.  100.)  GO  TO  218 
IF  ( AR$(W( IB+1, l> ).LT. 10000. )  GO  TO  216 
FFW(  3)=F.X2 
GG  *0  238 
216  F  E  W  (  3  )  -  F  X  4 
GC  TC  218 
218  FFW(3)=FX6 

238  CON” ! NUG 
J  F  =  1 

J  F  =  1 2 

I L°= I LP4 I L I 

239  WRITR  (6,125)  ILL, IT , TM , ( J , J * JB , JE >  „• 

WFITF  (6 , FMS )  (J,(  S7( I» J), I=JB,JE>, J*l,N>  . 

WRdF  (6,126)  ILL,  IT,TM,  (  J,  J»JB,  JE) 

WRITr  (6 , FMW )  (J,(  W(I,J),I*ja,JE),J=l,N) 

WF I"E  (6,136)  ILL, IT,TM. ( j.jsjR, JE) 

W  R  I T  E  ( 6  ,  FMY  )  (J,(YIT(  I,  J),  [»JR,  JE)  ,  J*l,rj) 

IF  ( ILL. ME. IF  )  GO  TO  240 
J  B  *  1 3 
JE  =24 
IE  =  !FdC  I 
GC  T  239 
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r»  o  o  o 


C 

C 


the  ir;r)Fx  jll  controls  the  saving  of  the  results  of  a  particular 

LOAn  r  ICRCNENT  on  THE  OUTPUT  TAPE  FOR  LATFR  REFERENCE.  THE 
INTERVAL  IS  GIVEN  BY  JLI.  THE  RESULTS  OF  THE  LAST  LUAO  INCREMENT 
IS  ALSO  SAVFO  FUR  THE  CONTINUATION  PUN. 

240  IF  ( ILL.NE.JLL.aND.ITIMEINT) .GT.IO)  GO  TO  250 
245  NTl«\*M«  I 

NLO(NTl) =  I LL 
SXT { 1,1) =  TM 
SYT  U , 1) « ILL 
JTM- I T IMF ( NT  I 

WRlTr  (NT 0 )  UT,unT,UOnT,UI»VT,VDT, VOOT.VI , SXT, SYT , SXYT , YI Tf SZ,W 
IF  (LPT.FO.I)  GO  TO  185 
IF  (  JTM. LT.  10)  GO  TO  99 
JLL=JLL+ JLI 
250  SI=St+SIK 
NT.?«MT1 

WP l~E  <6,138)  NOT 
Of)  205  1  =  1, NT1 
205  UP  I ”E  <6,  l  39)  NLO(I) 

THE  PFSULTS  FROM  THE  INPUT  TAPE  SAVED  IN  UNIT  8  IS  TRANSFFERRED  TO  THE 
OUTPUT  T  A  PE  BEHIND  THF  OUPUT  OF  THIS  RUN. 

IF  INT3.rO. 0)  GO  TO  199 
OP  210  1=1, N~3 

REA  )  ( NT  I)  UT.UDT.UODT.UI, VT,VOT,VOOT, VI.SXT, SYT, SXYT, YIT,SZ,W 
NT1*NT1* 1 
ILL=SYT< 1,1) 

WRITE  (NTO)  UT,unT,U00T,U!,VT,VDT,V9DT, V I , SXT, SYT , SXYT , YI T, SZ , W 
WRI~F  (6,139)  ILL 
IF  ( I T IMF  I  NT ) , LT • 3 )  GO  TO  199 
210  CONT  I NUE 
199  NT 3  =  NT  1 
LEN=LEN+ i 
ICV=0 

MMrE  (6,137)  LEN,  IL l , I F I , JL I , NT  1 ,NT2 , NT  3 , LPP , NTU , NT! , I C V 
99  RETURN 

119  FORMAT  ( 3X,23HSTRFSS  POINT  IS  PLASTIC//) 

125  FORMA”  ( 49H1 VERTICAL  STRESS  DISTRIBUTION  (SZ*100.)  AT  ILL  =,I5, 
16H,  I’  = , I  3 , 8  H,  TIME  =,F10.6,5H  SEC // 1 X , 19 , 1 1 1 10/ / ) 

126  FORMAT  ( 5RHIVFR7 ICAL  0ISpL4CFMENT  DISTRIBUTION  <W*10.0F  06)  AT  IL 
1L  =  ,  I  r> ,  6  H ,  IT  =,n,RH,  TIME  =,  F10.6, 5H  SEC //I  X ,  I  9 , 1 1 1 1 0// ) 

136  FORM*-  ( 40H1YJELP  INDICATING  TABLE  C  YIT)  AT  ILL  =,I5,6H,  IT  =, 
113, UH,  TIME  =,FI0.6,5H  S EC/ / IX , I  9, 1 1 1 1 0/ / ) 

137  FOR** A ”  (  1H0,3X,71HF0R  CONTINUING  RUN,  ONLY  NEFD  TO  CHANGE  THE  LAST 
l  Pft“4  CAPO  AS  FOLLOWS  —  // 

211X.57HLB  LEN  TL I  IE!  JLI  NT l  NT2  NT3  LPP  NT  I  NTO  ICV// 
38X»!5,5X,10I5) 

138  FORM  i  ”  (  1H0,18X,11HTAPE  NUMBER , 1 6, 54H  CONTAINS  THE  RESULTS  OF  LOA 
H)  INCREMENT  NUMBER  ILL  *//) 

139  FORM*.,r  (  23X,  15) 

142  F0R"4-  (ru.i.FU.S) 

143  FORMA”  ( 2  X , 2 ( F20 .8,418)  ) 

ENG 
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THE  FOLLOW  t  •-,'3  ! 3  A 


SAXPLF.  SFT  OF  INPUT  3ATA  FOR  THE  PROGRAM 


SOATA 

SPACING  FACTJR.  N=4.65 

130.0  0.45  8950.0  2000,0  15.0 

0.05  12.0  24600 .0  5.0CE-05  11  15 

3.00  6.25F-05  36  28 

0.0  -20.0  -110.0  -270.0  -580.0 

-1050.0  -1730.0  -2630.0  -3750.0  -4890.0 

-5890.0  -6770.0  -7540.0  -8190.0  -8730.0 

-9170.0  -9500.0  -9740.0  -9900.0  -9980.0 

-10000.0 

400.00  500.00 

930  2445 

1  50  5  10  50  0  0  0  1  9  10 

tEOF 


0 
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LIST  OF  SYMBOLS 


AL 

AP 

BE 

C 

CT 


Lame  constant  (i  )  in  psf;  later  becomes  dimensionless 
(X/k) 

Soil  parameter  (o') 

A  variable  in  the  plastic  relation  related  to  the  stress 
invariants  (B) 

Cohesion  (c)  in  psf 

The  value  of  the  cohesion  of  the  saved  data  on  tape;  used 
in  conversion 


Cl 

Cl 

DFU.DFV 


DT 
DTT 
DU,  DV 


Dilatational  wave  velocity  in  fps 
Shear  wave  velocity  in  fps 

The  value  of  the  maximum  percent  convergence  of  U  and 
V  displacements  between  successive  iterations  among 
all  the  mass  points 


Time  increment  (At)  in  seconds 
Dimensionless  time  increment 

Percent  convergence  of  the  U  and  V  displacements 
between  successive  iterations 


EX.EY.EXY  Strain  increments  (A€  ,  Af  ,  and  Ay  ,)  at  time  t 

n  C  T)  C 

FC  Yield  function 

FP  Load  curve  adjustment  factor 

FPL  Width  of  loaded  area  (a)  in  inches;  later  becomes  in  feet 

G  Shear  modulus  (G)  in  psf;  later  becomes  dimensionless 

(G/k) 

H  Grid  size,  the  distance  between  mass  points  (h)  ir 

inches;  later  becomes  dimensionless  (h/a) 

I  Index  in  the  horizontal  direction  (i) 

IB  I-index  of  the  mass  point  at  the  beginning  of  the  loaded 

area 


IC  Index  for  controlling  the  particular  surrounding  stress 

point  to  be  calculated 

IEI  The  number  of  load  increments  skipped  in  the  printout 

for  the  whole  region 

1EN  I-index  of  the  mass  point  at  the  end  of  the  loaded  area 


160 


I  LI 

ILP 

IT 

J 

JA 

KT 

LB 

LEN 

M 

N 

NIN 

P 

PCU 

PH 

PHI 

PIN 

PKP 

PO 

RAT 

RHO 

SI 

SIK 

SJ 

SS 

SX.SY.SXY 


The  number  of  load  increments  skipped  in  the  printout 
for  a  limited  region 

Index  for  controlling  which  load  increment  is  to  be  printed 
out 

Iteration  index 

Index  in  the  downward  direction  (j) 

Index  for  controlling,  during  the  final  iteration,  the 
entry  to  program  for  checking  if  the  stress  point  has 
yielded 

Index  for  indicating  the  particular  linear  segment  of  the 
load  curve  that  is  being  considered 

Starting  load  increment  number  for  the  particular 
computer  run 

Ending  load  increment  number  for  the  particular 
computer  run 

The  number  of  grid  points  in  the  horizontal  direction 

The  number  of  grid  points  in  the  downward  direction 

The  number  of  load  increments  in  a  linear  segment  of 
the  load  curve 

A  constant  in  the  plastic  relation  (Q) 

Magnitude  of  the  applied  pressure  after  each  linear 
segment  of  the  load  curve 

Stress  correction  factor  (T[  ) 

Frictional  angle  (<f>)  in  degrees;  later  becomes  in  radians 

Applied  pressure  increment  within  each  linear  segment 

Peak  load  of  load  pulse  in  psf 

Poisson's  ratio  (v  ) 

Stress  correction  factor  (]-•*"  ) 

3  C 

Weight  density  (p)  in  lb/ft 

The  current  applied  surface  pressure  that  is  prescribed 
at  the  fictitious  stress  points 

The  particular  applied  pressure  increment  being  considered 
Second  stress  invariant  of  the  stress  tensor  (J) 

First  stress  invariant  of  the  stress  tensor  (I) 

Normal  and  shear  stresses  in  the  (e.  ")  directions 

(0  ,  0,,  t  ,)  at  time  t 
-n  -v 


lol 


SXS,  SYS, 
SXYS 

SXT,  SYT, 
SXYT 

SZ 


SZZ 

TD 

TM 

TTD 

U, l'D,  UDD 

UI 

UPL 

UT,  UDT, 
UDDT 

V,  VD,  VDD 

VI 

VT.VDT, 

VDDT 

W 

WO 

WR 

YI 

YIT 

YS 


The  stresses  of  the  previous  iteration  at  the  surrounding 
stress  point  being  considered 

Normal  and  shear  stresses  in  the  (p,£)  directions 
(n  ,  n  ,  t  )  at  time  t-At 

n  c  nc  2 

Vertical  normal  stress  (0  )  x  10  ;  use  for  printout 
purpose 

Normal  stress  in  the  direction  normal  to  (r|,£) 

Time  duration  of  the  load  pulse  (t^)  in  seconds 
Time  in  seconds  (t) 

Dimensionless  time  (t/t^)  at  the  end  of  each  linear 
segment  of  the  load  curve 

Displacement,  velocity,  and  acceleration  in  the  n-direction 
(u,  u,  u)  at  time  t 

Displacement  increment  (Au)  in  the  p -direction  at  time  t 

A  variable  for  indicating  unloading;  UPL  =  -1.  0  is  not 
loading  and  UPL  =  1.0  is  unloading 

Displacement,  velocity,  and  acceleration  in  the  n-direction 
(u,  u,  u)  at  time  t-At 

Displacement,  velocity,  and  acceleration  in  the  ^-direction 
(v,  v,  v)  at  time  t 

Displacement  increment  (Av)  in  the  f-direction  at  time  t 

Displacement,  velocity,  and  acceleration  (v,  v,  v)  at 
time  t-At 

Vertical  displacement  (w)  x  10  ,  use  for  printout  purpose 
Incremental  plastic  work  done  (AW) 

A  reference  plastic  work  done 
Yield  indicating  table  at  time  t 
Yield  indicating  table  at  time  t-At 
Yield  stress  in  shear  (d)  in  psf 
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APPENDIX  IV 


MULTIPLE  WHEEL  VERIFICATION  TESTS 


Soil  Tests  and  Preparation 


Classification 


Both  of  the  soils  selected  for  ‘he  test  program  have  been  used 
extensively  by  WES  in  previous  mobility  studies  and  the  classification 
properties  have  been  reported  previously.  For  comparison  purposes, 
the  grain  size  distribution  and  limits  properties  are  given  for  the  buck¬ 
shot  clay  in  Figure  59,  and  the  grain  size  distribution  for  the  mortar 
sand  is  shown  in  Figure  60. 

California  Bearing  Ratio 

The  CBR  is  a  plate  bearing  test  using  a  three-square  inch  piston 
which  is  penetrated  continuously  into  the  soil  to  a  depth  of  one-half  inch 
while  continuously  recording  the  load  resistance  with  depth.  Annular 
surcharge  weights  are  placed  around  the  piston  prior  to  its  penetration. 
The  ratio  of  the  load  at  0.  1  inch  penetration  to  that  load  supported  bv  a 
standard  well  graded  crushed  gravel  multiplied  by  100  is  defined  as  the 
CBR  of  the  soil. 

Cone  Penetrometer  Resistance 

The  mobility  cone  penetrometer  is  a  rod  device  having  a  thirty 
degree  cone  tip  and  has  a  cross  section  base  area  of  0.  5  square  inch. 
The  shaft  is  narrowed  above  the  cone  to  minimize  the  friction  between 
the  side  of  the  shaft  and  the  hole.  The  cone  penetrometer,  which  is 
pushed  into  the  soil  at  a  standard  rate,  measures  the  resistance  to 
penetration  (Cone  Index)  in  pounds  per  square  inch.  The  Cone  Index 
is  a  measure  of  soil  shear  strength  and  its  variation  with  depth.  The 
Cl  value  is  usually  given  as  the  average  resistance  over  the  first  b" 
of  depth. 

Table  dn  gives  the  cone  index  values  taken  in  the  tire  ruts  after 
one  pass  ot  the  test  carriage.  For  tandem-tracking  tests,  the  rut  was 
was  formed  bv  both  tire*  and  the  strength  value  given  represents  the 
after  test  soil  strength  caused  by  t\i<>  passes  of  a  tire.  For  the  twin 


1  b  4 


GRADATION  CURVES 


NCHES  U  S.  STANDARD  SIEVE  NUMBERS  HYDROMETER 

i  4  6  8  10  14  16  20  30  40  50  70  100  140  200 


Soil  Strength  Consistency  and  After  Test  Strength 
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and  tandem-nont  racking  tests,  the  value  given  represents  only  one  pass 
soil  strength  alteration.  Also  included  in  this  table  is  the  before  traffic 
Cl  values  and  the  high  and  low  value  of  all  the  tests  to  indicate  the  soil 
bed  uniformity.  As  noted  previously,  the  CIa  indicates  that  the  values 
given  are  an  average  of  at  least  five  separate  tests. 

T  riaxial  T  ests 

Undrained,  unconfined  triaxial  tests  were  run  on  undisturbed  clay 
specimens  taken  from  the  test  soil  carts  for  the  soil  conditions  of  0=20. 
The  triaxial  specimens  with  zero  confining  pressure  were  approximately 
1.7  5  inches  in  diameter  by  4.  5  inches  high.  Each  test  specimen  was 
extruded  from  a  piston  type  soil  sample  similar  to  shelby  tube  samples. 
Only  the  specimen  ends  needed  trimming  before  testing.  The  results  of 
the  unconfined  tests  are  plotted  in  Figure  61. 

Triaxial  tests  were  also  conducted  on  air  dried  samples  of  the 
mortar  sand.  The  specimen  size  was  2.  8  inches  in  diameter  by6.  0  inches 
high  and  the  samples  were  prepared  to  a  unit  weight  of  98.0  pcf  ±1.0  pcf 
which  was  the  average  soil  unit  weight  in  the  soil  carts  for  the  CI=40 
test  condition.  The  results  of  the  triaxial  tests  on  the  sand  are  shown  in 
Figure  62. 

Soil  Preparation 

The  buckshot  clay  was  processed,  placed,  and  compacted  at 
predetermined  moistures  and  densities  corresponding  to  the  desired  soil 
strength.  The  soil  is  first  passed  through  a  roller  crusher  and  then 
placed  in  a  pug  mill  where  the  soil-water  mixing  takes  place  at  a  selected 
moisture  content.  The  soil  is  then  placed  in  the  soil  carts  in  6  inch 
layers  with  each  layer  compacted  by  pneumatic  tired  rollers.  Previously 
developed  empirical  relationships  between  buckshot  clay  moisture  content 
and  eompactive  effort  permitted  the  soil  to  be  placed  near  the  design  soil 
conditions. 
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Buckshot  Clay 
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Figure  61.  Triaxial  Results,  Buckshot  Clay 


Figure  62.  Triaxial  Results,  Mortar  Sand 


The  mortar  sand  was  prepared  in  an  air  dried  condition.  Different 
soil  strengths  were  achieved  by  varying  the  density  of  placement  of  the 
soil.  The  sand  was  placed  in  uniform  layers  which  were  screened  and 
vibrated  on  the  surface  as  the  filling  progressed.  Empirical  relationships 
between  the  thickness  of  layer,  vibratory  effect,  and  soil  density  permitted 
the  sand  to  be  placed  near  the  design  soil  conditions. 

Test  Procedures 

The  test  procedures  for  running  a  single  test  in  sand  are: 

(1)  Loosen  sand  section  by  plowing 

(2)  Check  weakened  state  by  Cone  Index  tests 

(3)  Vibrate  sand  to  desired  soil  strength 

(4)  Check  strength  by  Cone  Index  tests 

(5)  Take  final  soil  strength  and  surface  elevation  profile 

(6)  Calibration  checks  between  recording  station  and  computer 

(7)  Make  a  test  run  without  load  on  test  tires  (in-air  run) 

(8)  Load  tire  and  check  inflation  pressure  and  deflection 

(9)  Run  test 

(10)  Take  post  test  soil  strength  and  rut  depth  profile. 

The  test  procedure  for  running  a  test  in  clay  differs  from  the  sand 
only  in  soil  preparation  techniques.  The  clay  test  bed  has  to  be  prepared 
at  least  ten  days  in  advance  of  a  test  so  that  the  water  content  of  the  clay 
will  stabilize  as  described  above.  The  only  preparation  needed  to  run  a 
test  after  the  soil  bed  has  been  constructed  involves  smoothing  the  soil 
surface  and  measuring  the  soil  strength  and  surface  profiles.  Steps  (7) 
through  (10)  are  then  completed  as  in  sand.  Preparation  for  succeeding 
tests  in  the  same  clay  soil  bed  can  involve  re-rolling  and  smoothing  the 
soil  surface.  In  both  the  sand  and  clay  test  procedure,  a  test  is  not  run 
unless  the  desired  soil  strength  and  uniformity  is  attained  during  the 
preparation  stages. 
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Sinkage  Measurement 


The  sinkage  measurement  determined  for  each  clay  test  is  the 
result  of  the  use  of  deflection  pins  which  were  developed  in  previous 
testing  programs  such  as  the  C-5A  program.  The  devices  consist  of  a 
flat  circular  plate  (2"  diameter)  and  a  short  1/4"  diameter  rod,  pointed 
on  one  end.  Before  each  test,  a  rod  is  penetrated  into  the  soil  and  then 
the  flat  plate  is  placed  on  top  of  the  pin  flush  to  the  soil  surface  in  the 
tire  path.  After  the  test,  the  differential  movement  of  the  pin  and  plate 
is  the  value  of  the  elastic  sinkage  (rebound).  Normally,  the  pin  does 
not  rebound  with  the  soil,  but  the  flat  plate  normally  does. 


APPENDIX  V 

WES  BRAKING  TEST  DATA 
SAND  AND  CLAY 
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UNIVERSITY  OF  DAYTON  TIRE/SOIL  INTERACTION 
RESEARCH  REPORTS  AND 
COMPUTER  PROGRAMS 
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A  summary  listing  is  given  on  the  following  pages  of  each  report 
and  each  computer  program  ceveloped  by  the  University  of  Dayton  under 
Air  Force  sponsorship  (Air  Force  Flight  Dynamics  Laboratory,  Vehicle 
Equipment  Division)  in  the  research  program,  "Landing  Gear/Soils 
Interaction  and  Flotation  Criteria," 

The  computer  programs  are  available  for  use  by  other  organizations 
with  Air  Force  permission.  Additional  information  may  be  secured  by 
contacting: 

Dr.  David  C.  Kraft 

Dept,  of  Civil  Engineering  and  Research  Institute 
University  of  Dayton 
Dayton,  Ohio  45409 

Mr.  George  J.  Sperry 

Project  Engineer,  Vehicle  Equipment  Division  (AFFDL/FEM) 
Air  Force  Flight  Dynamics  Laboratory 
Wright-Patterson  Air  Force  Base 

LANDING  GEAR/SOILS  INTERACTION  AND  FLOTATION  CRITERIA 
PUBLICATIONS  AND  REPORTS 

Kraft,  David  C.  ,  Analytical  Landing  Gear-Soils  Interaction,  Pha^e  I, 
AFFDL-TR-68-88,  Air  Force  Flight  Dynamics  Laboratory,  Wright- 
Pattcrson  Air  Force  Base,  Ohio,  August  1968. 

Kraft,  David  C.  ,  Preliminary  -  ingle  Wheel  Relative  Merti  Index, 

UDRI-TR -69- 16,  University  of  Dayton,  Dayton,  Ohio,  May  1°69. 

Kraft,  David  C.  ,  Flotation  Performance  of  Aircraft  Tiros  on  Soil  Rur.wayt  . 
Journal  of  Terramechaniv  s,  Vol.  6,  No.  1,  1969. 

Kraft,  D.  C.  ;  Luming,  H.  ;  Hoppcnjans,  J.  R.  ,  Aircraft  Landing  Giv.r- 
Soils  Interaction  and  F’otation  Criteria,  Phase  II,  AFFDL-TU -6‘>-?h. 

Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson  AFB,  Ohio, 
November  1969. 

Luming,  Henry,  F  rite  Element  Approach  to  Axisymmet, Dynamic  S  uil 
Deformations,  Symposium  on  Application  of  Finite  Element  Methods  in 
Civil  Engineering,  *  andcrbilt  University,  Nashville,  Tennessee,  November, 

1969. 
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Kraft,  David  C.  ,  Liguori,  Albert  E. ,  Hoppenjans,  J.  Richard,  Twin- 
Vertical  Plate  Verification  Tests,  Test  Report,  UDRI-TR-70-27,  University 
of  Dayton  Research  Institute,  Dayton,  Ohio,  May,  1970. 

Luming,  Henry,  Multiwheel  Vertical  Pulse  Load  Analytical.ginkage 
Prediction  Technique  and  Computer  Program,  UDRI-TR-70-22,  University 
of  Dayton  Research  Institute,  Dayton.  Ohio,  May,  1970. 

Luming,  Henry,  Analytical  Aircraft  Landing  Gear-Soil  Interaction  - 
Phase  III,  Rolling  Single  Wheel  Analytical  Sinkage  Prediction  Technique 
and  Computer  Program,  AFFDL-TR-70- 142,  Air  Force  Flight  Dynamics 
Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio,  September,  1970. 

Kraft,  David  C.  ,  and  Luming,  Henry,  Multiple  Rolling  Tire  Sinkage  and 
Drag  Interaction  Effects,  accepted  for  presentation  at  Joint  ISTVS-SAE 
Meeting,  Detroit,  Michigan,  January,  1971. 

LANDING  GEAR/SOILS  INTERACTION  AND  FLOTATION  CRITERIA 
COMPUTER  PROGRAMS 

Title  of  Computer  Program:  Transient  Loading -Sinkage  Analysis,  Com¬ 
puter  Program  -  1 

Brief  Description:  The  computer  program  calculates  the  instantaneous 
(time-dependent)  sinkage  into  a  soil  medium  of  a  rolling  aircraft 
tire.  The  rolling  aircraft  tire  loading  is  simulated  as  a  dynamic 
pulse  loading  applied  in  a  vertical  direction  through  a  mass  at  the 
interface.  The  duration  of  the  pulse  is  varied  to  simulate  different 
forward  velocities  of  the  tire.  The  soil  medium  is  assumed  to  be 
elastic  and  the  load  is  applied  as  a  uniform  pressure  over  a 
circularly  loaded  area.  The  input  parameters  are  the  magnitude  of 
the  mass,  shape  of  load  pulse,  duration  of  pulse,  radius  of  loaded 
area,  intensity  of  pressure,  soil  density,  and  soil  shear  modulus. 

Computer  Language:  Fortran  IV  -  (IBM) 

Equipment:  The  computer  program  was  originally  written  for  use  on  the 
IBM  7094  at  Wright-Patter son  AFB,  Dayton,  Ohio. 
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Reference:  Kraft,  David  C.  ,  "Analytical  Landing  Gear-Soils  Interaction  ~ 
Phase  I,  "  AFFDL-TR-68-88,  Air  Force  Flight  Dynamics  Labora¬ 
tory,  Wright-Patterson  AFB,  Ohio,  August  1968. 

Title  of  Computer  Program:  Flotation  Index-Operations  Inuex,  Computer 
Program  -2 

Brief  Description:  The  computer  program  calculates  the  flotation  capacity 

of  single  and  multiple  wheel  landing  gear  configurations  for  operation 
on  unprepared  (soil)  runways.  The  flotation  capacity  is  expressed 
by  the  Flotation  Index  (FI)  and  Operations  Index  (OI)  which  are 
calculated  based  on  sinkage  and  drag.  The  FI  is  the  drag  ratio  of 
a  given  aircraft  based  on  specified  operating  conditions.  The  OI  is 
the  ratio  of  sinkage  to  load  at  the  same  operating  conditions. 

Current  program  results  include  flotation  ratings  of  all  currently 
used  aircraft  tires  on  cargo,  bomber,  and  fighter  type  aircraft. 

These  results  permit  aircraft  designers  to  select  tires  and  landing 
gear  configurations  for  optimum  flotation  (minimum  drag).  Program 
was  revised  6-70. 

Computer  Language:  Fortran  IV  -  (IBM) 

Equipment:  The  computer  program  was  originally  written  for  use  on  the 
IBM  7094  at  Wright-Patterson  AFB,  Dayton,  Ohio. 

Reference:  Kraft,  David  C.  ,  Luming,  Henry,  and  Hoppenjans,  J.  R.  , 

"Aircraft  Landing  Gear-Soils  Interaction  and  Flotation  Criteria  - 
Phase  II,  "  AFFDL-TR-69-76,  Air  Force  Flight  Dynamics 
Laboratory,  Wright-Patterson  AFB,  Ohio,  November  1969. 

Title  of  Computer  Program:  Sinkge  Wheel  Stationary  Pulse  Load-Sinkagc 
Prediction,  Computer  Program  -  3 

Brief  Description:  This  computer  program  calculates  the  instantaneous 
sinkage  into  a  soil  medium  caused  by  a  rolling  aircraft  lire.  The 
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interface  contact  of  the  rolling  tire  is  simulated  by  a  stationary- 
circular  surface  contact  pressure  which  is  uniform  over  the  con¬ 
tact  area  and  varies  with  time  in  the  form  of  a  pulse.  The  mag¬ 
nitude  of  the  pressure  changes  in  the  same  manner  as  the  pressure 
experienced  by  a  soil  particle  near  the  surface  of  the  soil  as  the 
tire  rolls  over  it.  The  soil  is  assumed  to  be  an  clastic-plastic 
material  with  elastic  deformation  governed  by  Hooke's  law,  the 
plastic  deformations  governed  by  an  incremental  stress-strain 
relation  which  is  based  on  the  normality  flow  rule,  and  the  plastic 
yielding  governed  by  the  Drucker-Prager  yield  criterion  with  no 
strain-hardening.  The  input  soil  parameters  are  the  density,  the 
Young's  modulus,  the  cohesion,  and  the  friction  angle.  The 
numerical  method  used  in  solving  the  boundary  value  problem  is 
the  lumped  parameter  iteration  method.  This  method  uses  an 
axisymmetric  lumped  parameter  model  to  approximate  the  con¬ 
tinuous  medium  and  an  iterative  procedure  to  calculate  the  dis¬ 
placements  and  the  stresses  at  the  discrete  points  of  the  model. 

Computer  Language:  Fortran  IV  -  (IBM) 

Equipment:  This  computer  program  was  originally  written  for  use  on  the 
IBM  7094-DCS  at  Wright-Patter son  AFB,  Dayton,  Ohio.  It  has  a 
32K-word  core  capacity. 

Reference:  Kraft,  David  C.  ,  Luming,  Henry,  and  Hoppenjans,  J.  R.  , 

"Aircraft  Landing  Gear-Soils  Interaction  and  Flotation  Criteria  - 
Phase  II",  AFFDL-TR -69-76,  Air  Force  Flight  Dynamics 
Laboratory,  Wright-Patter  son  AFB,  Ohio,  November  1969. 

Title  of  Computer  Program:  Rolling  Single  Wheel  Sinkage  Prediction, 
Computer  Program  -  4 

Brief  Description:  This  computer  program  calculates  the  instantaneous 

sinkage  into  a  soil  medium  caused  by  a  rolling  aircraft  tire.  The 
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interface  contact  of  the  rolling  tire  is  simulated  by  a  surface  contact 
pressure  which  is  applied  uniformly  over  an  area  equivalent  to  the 
tire  footprint  area  and  is  moved  across  the  surface  at  the  aircraft 
horizontal  ground  velocity.  The  magnitude  of  the  uniform  pressure 
increases  over  a  finite  rise  time  from  zero  to  a  pressure  equal  to 
the  vertical  tire  load  divided  by  the  contact  area.  The  soil  medium 
is  assumed  to  be  an  elastic-plastic  material  with  elastic  deformations 
governed  by  Hooke's  law,  the  plastic  deformations  governed  by  an 
incremental  stress-strain  relation  which  is  based  on  the  normality 
flow  rule,  and  the  plastic  yielding  governed  by  the  Drucker-Prager 
yield  criterion  with  no  strain-hardening.  The  input  soil  parameters 
are  the  density,  the  Young's  modulus,  the  c  ohesion,  and  the  friction 
angle.  The  numerical  method  used  in  solving  the  boundary  value 
problem  is  the  lumped  parameter  iteration  method.  This  method 
uses  a  three-dimensional  lumped  parameter  model  to  approximate 
the  continuous  medium  and  an  iterative  procedure  to  calculate  the 
displacements  and  stresses  at  the  discrete  points  of  the  model.  An 
input-output  scheme  is  also  used  for  utilizing  the  limited  core 
capacity  for  the  three  dimensional  problem. 

Computer  Language:  Fortran  IV-(IBM) 

Equipment:  This  computer  program  was  originally  written  for  use  on  the 
IBM  7094-DCS  at  Wright-Patterson  AFB,  Dayton,  Ohio.  It  has  a 
32K-word  core  capacity. 

Reference:  Luming,  H.  ,  "Analytical  Landing  Gear-Soil  Interaction  -  Phase 
III,  Rolling  Single  Wheel  Analytical  Sinkage  Prediction  Technique  and 
Computer  Program,"  AFFDL-TR-70-  142,  Air  Force  Flight 
Dynamics  Laboratory,  Wright-Patterson  AFB,  Ohio,  September  1970. 

Title  of  Computer  Program:  Multiwhocl  Stationary  Pulse  Load  Sinkage 
Prediction,  Computer  Program  -  5 
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Brief  Description:  This  computer  program  calculates  the  instantaneous 

sinkage  into  a  soil  medium  caused  by  a  pair  of  rolling  twin  aircraft 
tires.  If  the  sinkages  for  various  twin-wheel  spacings  were  cal¬ 
culated  and  compared  with  the  corresponding  single-wheel  sinkage, 

%  twin  wheel  effects  of  aircraft  landing  gear  configurations  could  be 
obtained  for  use  in  flotation  studies.  The  interface  contacts  of  the 
rolling  twin  tires  are  simulated  by  two  stationary  surface  pressure 
strips  which  are  uniformly  distributed  and  vary  with  time  in  the 
form  of  a  pulse.  The  magnitude  of  the  pressure  changes  in  the  same 
manner  as  the  pressure  experienced  by  a  soil  particle  near  the 
surface  of  the  soil  as  the  tire  rolls  over  it.  The  soil  is  assumed 
to  be  an  elastic-plastic  material  with  elast.c  deformations  governed 
by  Hooke's  law,  the  plastic  deformations  governed  by  an  incremental 
stress-strain  relation  which  i3  based  on  the  normality  flow  rule, 
and  the  plastic  yielding  governed  by  the  Drucker-Prager  yield 
criterion  with  no  strain-hardening.  The  input  soil  parameters  are 
the  density,  the  Young's  modc.los,  the  cohesion,  and  the  friction  angle. 
The  numerical  method  used  in  solving  the  boundary  value  problem  is 
the  lumped  parameter  iteration  method.  This  method  uses  a  plane- 
strain  two  dimensional  lumped  parameter  model  to  approximate  the 
continuous  medium  and  an  iterative  procedure  to  calculate  the  dis¬ 
placements  and  stresses  at  the  discrete  points  of  the  model. 

Computer  Language:  Fortran  I\  -(IBM) 

Equipment:  This  computer  program  was  originally  written  for  use  on  the 
IBM  7094-DCS  at  Wright-Patterson  AFB,  Dayton,  Ohio.  It  has  a 
32K-word  core  capacity. 

Reference:  Luming,  Henry,  "Multiwheel  Landing  Gear-Soil  Interaction  - 
Phase  III,  Multiwheel  Vertical  Pulse  Load  Analytical  Sinkage  Pre¬ 
diction  Technique  and  Computer  Program,"  RI«D  Computer  Program 
Report,  Report  No.  UDRI-TR-70-22,  University  of  Dayton  Research 
Institute,  Dayton,  Ohio,  May  1970. 
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